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HEAT  AND  MASS  TRANSFER  PARAMETERS 

THERMOPHYSICAL  CHARACTERISTICS  OF  HEAT- SENSITIVE  PAINTS 
B.  G.  Abramovich  and  L.  N.  Novichenok 

In  recent  years  heat-sensitive  paints  have  found  broader  use  in 
measuring  surface  temperatures  both  in  our  country  and  abroad.  This  method 
of  measuring  temperature  has  one  adherence  by  its  simplicity,  profitability, 
and  its  truly  endless  possibilities  of  application.  What  point  measurements 
achieve  by  protracted  laborious  investigations,  thermo-paiuts  arrive  at 
almost  without  any  labor  Involved.  They  make  it  possible  to  obtain  a 
graphic  picture  of  temperature  distribution  with  various  units  used  in 
measuring  temperature  on  moving  and  current-bearing  parts,  etc.  Balance 
study  of  the  properties  of  these  remarkable  temperature  indicators  throw 
open  new  vistas  of  their  application. 

To  evaluate  the  inertness  of  thermo-paints  and  their  effect  on  heat 
transfer,  we  determined  thermophysical  properties  of  a  set  of  heat-sensitive 
paints  which  are  now  being  readied  for  series  production. 

In  exploring  thermophysical  properties  of  heat- sensitive  paints,  we  used 
a  nonsteady-state  method  [2,  3]  the  basis  of  which  is  the  solution  of  the 
thermal  conductivity  problem  for  a  system  of  bodies  consisting  of  a  bounded 
and  a  semibounded  rod  with  boundary  conditions  of  the  first  and  fourth 
orders  [1], 

Thermophysical  characteristics  were  calculated  from  the  data  of  one 
short  experiment,  using  the  following  formula#: 
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for  Che  coefficient  of  thermal  diffusivity 


0  *  <2> 

for  specific  heat 

c-W 

(3) 


where  k  was  determined  from  previously  compiled  calculation  tables  [3]. 

A  schematic  diagram  of  the  apparatus  is  shown  in  Figure  1. 

The  specific  weight  of  a  thermo-paint  film  was  determined  by  the  method 

3 

of  titrating  a  potassium  Iodide  solution  CyAI  -  1500  kg/m  ),  in  which  the 
thermo-paint  film  was  placed,  using  water  [4].  The  method  was  based  on 
establishing  a  relationship  between  the  value  of  the  specific  weight  of  the 
film  and  the  ratio  of  the  volumes  of  titrated  and  titrating  liquids  at  the 
moment  the  film  falls  onto  the  bottom  of  the  vessel. 

The  calculation  made  use  of  the  formula 


v, 


to. 


All  the  results  of  the  study  are  shown  in  tabular  form. 


ThermophyBical  characteristics  of  thermo-paints  were  determined  with  an 
error  of  approximately  10-15  percent,  and  the  error  in  determination  of  their 
density  was  approximately  5  percent. 

As  we  can  see  from  the  table,  the  values  of  the  thermal  conductivity 

coefficient  lie  within  the  limits  0.19-0.99  w/m* deg,  the  coefficient  of 

-7  2 

thermal  diffusivity  ~  in  the  limits  (1. 1-29) ■ 10  m  /sec,  and  the  specific 
heat  —  <0. 22-1.45) *103  J/kg*deg. 
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Table  of  Thermophysical  Characteristics  of  a  Set  of  Heat-Sensitive  Paints 


i 

1 


j 

5* 


Critical 


No. 

No.  of 
:  thermo 
paint 

tempera* 

ture»  l« 
°C  at  cr  ' 

•^=120 

see  and 

V30 

d  sec 

Density 

v  kC 

»  m' 

Coeffi¬ 
cient  of 
thermal 
conduct¬ 
ivity 

A  nlf*o9 

Specific 

heat 

oc-^ 
kg* deg 

Coefficient 
of  thermal 
dlffuslvlty 

1.4 

I 

I  A 

45 

1445 

0,22 

1,12 

2 

17 

55 

1483 

0,22 

1,10 

1.4 

3 

31 

70 

1340 

0,19 

1,30 

I.I 

4 

40 

74 

1435 

0.24 

1,035 

1.6 

5 

32 

90 

1395 

0,21 

1,065 

1,4 

6 

50 

no 

1428 

0,23 

1,132 

1.4 

7 

4 

130 

1415 

0,22 

1,350  * 

I.I 

8 

30 

130 

1460 

0,23 

1,234 

1,3 

9 

19 

145 

1455 

0,26 

0,900 

.  2,1 

10 

5 

160 

1395 

0,38  • 

0.587 

4,6 

II 

29 

175 

1305 

0,83 

0,299 

19,9 

12 

6 

190 

1200 

0,27 

1,070 

2,1 

13 

18 

205 

1390 

0,43 

0,585 

5,3 

14 

230 

220 

1460 

0,62 

0,456 

7,8 

15 

240 

225 

1450 

0,46 

0,480 

6.7 

16 

8 

240 

I486 

0,36 

0,635 

3,9 

17 

320 

250 

1407 

0,31 

0,754 

2,9 

18 

35 

280 

1220 

!  0,26 

1,040 

2,1 

19 

33 

320 

1378 

0,24 

1,235 

1,4 

20 

26 

350 

1420 

0,24 

0,891 

1.9 

21 

12 

360 

1445 

0,84 

0,260 

20.6 

22 

230* 

400 

1455 

0,52 

0,455 

7.8 

23 

470 

440  ■ 

1420 

0,39 

0,524 

6,2 

24 

480 

480 

1445 

0.3Q 

0,838 

2,5 

25 

490 

490 

1430 

0,34 

0,700 

3,4 

26 

14 

550 

1455 

0,99 

0,229 

29.0 

27 

52 

555 

1450 

0,17 

0,614 

1.5 

28 

30 

570 

1460 

0,23 

1.234 

1.3 

29 

7 

£0 

1421 

0,53 

0.5*5 

7,0 

90 

!00 

830 

1978 

0,21 

1.450 

I.I 
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For  comparison,  we  present  the  vein**  of  thermephysical  characteristic,, 
inherent  In  the  following  materials: 

a)  lacquer  and  paint  coatings  and  resins  (organosllicon  lacquer  K-55, 
epoxide  resin  E-41,  etc.) 

X  -  0.11-0.17  w/m’deg, 
a  -  (0.43-0. 87)*10-7  m2/sec, 
c  »  (1. 4-1. 8) • 103  J/kg*deg; 

b)  heat-insulation  materials  (asbestos  cement,  asbestos  bakelite,  poly- 
chlorovinyl  plastic,  etc.) 

X  »  0.032-0.11  w/nwdeg, 
a  -  (0.2  -  0.4)  10-7m2/sec, 
c  -  (32  -  1.85)  103  J/kg*deg; 

c)  stainless  steel  lKhl8N9T 

X  ■*  9.5  -  12  w/m*deg, 

a  -  55*10-7m2/sec , 

c  -  (0.22-0.28)  103  J/kg-deg; 

From  comparison  of  this  data  it  is  clear  that  heat-sensitive  paints 
have  higher  coefficients  of  thermal  conductivity  and  thermal  diffusivity 
than  ordinary  lacquer-paint  coatings' and  heat-insulating  materials,  but 
lower  values  than  steel. 

As  a  result  of  studies  of  therraophys ical  properties  of  thermo-paints,  it 
has  become  possible  to  calculate  their  dynamic  error,  which  is  approximately 
0.1  percent,  and  to,  establish  that  thermo-paints  practically  speaking  do  not 
distort  the  temperature  field. 

Accordingly,  heat-sensitive  paints  can  be  regarded  as  low-inertia 
temperature  transducers  suitable  for  investigating  thermal  fields  of  various 
materials. 

Symbols 

9  -  relative  temperature;  X,  a,  c  and  y  *  coefficients  of  thermal 

I 

conductivity,  thermal  diffusivity,  heat  capacity  and  specific  weight  of  heat- 
sensitive  paint;  b  -  thermal  activity  of  etalon  material;  R  -  thickness  of 


9 


i 

I 

i 

i 

i 


I 

! 


X 


A 


thermo-paint  layer;  t  -  time;  ■  specific  weight  of  water; 
water  required  in  titration;  ■  volume  of  KT  solution. 


i 
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volume  of 


Figure  1.  Schematic  diagram  of  apparatuat  1,  Thermostat;  2,  Base 
3,  Vessel;  4,  Thermocouples;  5,  Specimen;  6,  Self-recording 
potentiometer;  7,  Holder;  8,  Thermocouple  switch 
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THEKMOPHYSICAL  CHARACTERISTICS  OF  FERRITES 
Ya.  M.  fiekker  and  A.  F.  Chudnovakiy 

The  special  physical  properties  of  ferrites ,  in  particular,  high 
electrical  resistance,  make  it  possible  to  directly  study  the  lattice 
component  of  thermal  conductivity,  and  also  effects  associated  with  magnetic 
ordering  and  crystallochemlcal  structure. 

The  present  study  deals  with  thermal  conductivity  and  thermal  expansion 
of  ferrites,  and  establishment  of  their  interrelationship  with  other  physical 
properties. 

Materials  and  Methods  of  Their  Investigation 

1.  Choice  of  chemical  compositions  of  the  materials  in  the  study  was 
dictated  by  the  necessity  to  vary  within  wide  limits  crystallochemlcal 
structure  and  magnetic  properties  of  solid  ferrite  solutions.  For  this 
purpose,  three  theories  of  ferrite  compositions  were  synthesized  containing 
cations  with  predominant  localization  in  octahedra  (Al,  Gr,  Ni)  or  in 
tetrahedra  (Zn,  Cd),  and  also  cations  (Mn,  Fe,  Mg),  the  localization  of  which 
was  determined  by  technological  factors. 

Series  1  consisted  of  solid  solutions  of  ferrite-chromites: 

Fe(t.l)1Crt(,.v)]04  , 

where  Q  <  y  (  1  is  the  concentration  of  solid  solution  (varied  by  0.1  from 
specimen  to  specimen); 

0  «  i,  <  1  is  the  extent  of  spinel  inversion.  Me  is  ions  of  the  metals 
Ml,  Mg,  Cd. 

Series  two  —  solid  solutions  of  ferrite-magnesium  aluminates: 


I': 
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Series  three  —  solid  solutions  of  ferrites: 


y  Ni  fea04  —  FejQ* 

where  Me  is  ions  Mg,  Zn,  Mn. 

To  exclude  the  effects  of  technological  factors  on  measurement  results, 
all  specimens  of  each  ferrite  system  were  prepared  under  the  aame  techno¬ 
logical  regime,  and  heat  treatment  was  conducted  simultaneously  and  jointly. 

2.  The  coefficient  of  thermal  conductivity  was  measured  by  a  method  in 
[1]  at  room  temperature,  by  a  method  in  [2]  in  a  wid»  temperature  range,  and 
by  a  method  in  [3]  under  exposure  to  an  external  magnetic  field.  The 
coefficient  of  thermal  expansion  was  measured  on  a  multi-purpose  vacuum 
dilatometer.  The  phase  composition  of  the  specimens  and  thair  crystalline 
structure  were  studied  on  the  MIM-7  microscope  (porosity,  granularity),  on 
the  URS-50I  X-ray  unit  (lattice  constant,  phase  composition),  and  pn  the 
YeM-6A  electronic  microscope  (structure  of  crystallites). 

Magnetic  parameters  in  static  fields  were  measured  by  the  ballistic 
method.  The  Curie  temperature  was  measured  on  the  NIM-4  bridge.  The 
electrical  conductivity  and  energy  of  activation  were  determined  by  the 
four-probe  method.  The  rate  of  propagation  of  sound  longitudinal  and 
transverse  waves  in  the  specimens  was  measured  on  the  UZIS-10  instrument. 
Microhardness  of  specimens  was  determined  on  the  PMT-3  instrument. 

Results  of  Measurements  and  Their  Discussion 

1.  Thermal  conductivity  as  a  function  of  mean  grain  size  is  shown 

in  Figure  1. 

Mlcrostructural  studies  revealed  that  initially  grain  size  growth  is 
accompanied  by  reduction  of  its  defect  content  and  prolongation  of  inter- 
crystalline  boundaries,  which  leadB  to  a  rise  in  K.  This  is  In  agreement 
with  magnetic  measurements.  A  structure-sensitive  magnetic  character¬ 
istic  —  the  coercive  force  (H  )  —  decreased  with  decrease  in  defect 

c 

concentration  [4]. 


Figure  1  (curve  J‘)  makes  clear  that  the  ljuik  ul  at*  a  function  of  l n 
corresponds  to  a  drop  In  K.  The  appearance  of  new  defects  accompanied 
further  grain  size  growth.  Electron-microscopic  studies  showed  that  for 
coarsely  crystalline  structure  (80-120  micromicronu)  the  specific  defect 
volume  rises  sharply,  which  reduces  the  value  of  K. 

A  relationship  exists  between  K  and  Ac,  as  predicted  by  theory. 

Results  of  experimental  study  of  this  relationship  are  Bhown  in  Figure  2, 
where  the  calculated  curve:  K  is  shown  for  purposes  of  comparison. 

It  is  clear  from  this  that  the  relationship  k**Aen  is  approximately  satisfied 
where  0,6  4  n  <  1,5  In  compositions  with  predominantly  ionic  bond  form  A£ 
proved  to  have  a  relatively  weak  (n  <  1)  effect  on  K.  Ac  more  markedly 
affects  (n  >  1)  K  in  those  compositions  where  locally  directed  covalent 
bonds  are  predominant  and  are  stronger  than  ionic  bonds. 

3.  The  crystal  lattice  constant  (sq)  reflects  [5]  the  energy  state  of 
the  given  crystalline  structure  and,  consequently,  affects  6  and  k  of 
ferrites.  The  mean  free  phonon  path  X c 


decreases  with  rise  in  aQ,  that  is,  k  is  reduced.  Additionally,  the 
coefficient  of  anhazmonicity  (8)  rises  with  rise  in  aQ,  which  Increases  the 
phonon  scattering  cross  section  (S ($)')  and  5.  Figure  3  presents  the  results 
of  study  of  6  and  the  reduced  thermal  resistance  (W^  -  W/Ac)  as  a  function  of 
a  .  It  is  clear  from  the  figure  that  experimental  results  are  in  agreement 
with  theory.  The  functions  derived  have  been  confirmed  by  calculations  of 
the  Debye  temperature  9a  *232  '-ras-  (here  the  minimum  spacing  between  ionB 
in  ferrites  /is  taken  into  account).  0D  is  reduced  with  rise  in  aQ,  that 
is,  the  number  of  phonons  in  the  state  of  statistical  equilibrium 
Ntp**  (exp-3j?--*i  )"*  rises  and  increases  their  scattering.  0^  for  the 

composition  studied  lies  within  the  limits  450  <  700°K. 


4.  Anal/ala  of  the  thermal  conductivity  components  of  ferriteB  In  a 
magnetic  field  showed  that  the  variation  in  K  in  the  temperature  range 
studied  can  be  caused  by  the  following: 

a)  the  Maggi-Righi-Leduc  effect  (MRL)  and  the  a  >  AE  effect. 

Estimation  of  the  contribution  made  by  the  MRL  effect  to  the  value  of  K  for 

iron  ferrite,  where  this  contribution  must  be  the  maximum,  shows  that  the 

presumed  effect  is  caused  by  [7]  a  supplementary  temperature  gradient 

_2 

approximately  equal  to  10  deg/cm,  which  is  three  orders  of  magnitude  less 
than  the  applied  gradient. 

Experimental  data  shown  in  Figure  4  demonstrates  that  the  thermal 
conductivity  of  ferrites  in  a  magnetic  field  initially  rises  .ftsZKg-—6  >0 
and  then  remains  unchanged.  The  quantity  AK^  —  the  effect  for  cobalt 
ferrite  is  greater  than  for  nickel  ferrite,  the  magnetization  of  which  is 
less  than  for  the  cobalt  ferrite.  These  results  can  be  accounted  for  by  the 
AE  effect  [8],  the  value  of  which  is  determined  from  the  relationship 


where  B  is  a  numerical  constant, 
o 

The  influence  of  the  AE  effect  lies  in  the  fact  that  when  ferrites  are 

magnetized,  their  modulus  of  elasticity  and  shear  modulus  determining  the 

rate  of  phonon  propagation  are  changed.  The  value  of  the  AE  effect 

Increases  with  rising  H  and  when  I  *  I  remains  unchanged,  AE  ■  ■  0.66  B  E  , 

8  0  0 

analogous  to  the  observed  function  of  K^. 

If  we  represent  the  AE  effect  [9]  via  microscopic  characteristics 
of  ferrites 

aE  2/U  Ef 
So  5©i 

then  it  becomes  understandable  that  the  value  of  K  of  ferrites  in  a  magnetic 
field  depends  materially  on  their  magnetostriction  (\g)  since  internal 
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hut  for 


stresses  [10]  6i  »  104  M/i»  .  Actually  [11]  for  CoFe^*,.  “  200*10 

NlFejO^,  Xg  “  -26*10  6,  respectively,  which  is  in  agreement  with  the 

AKl  of  these  farritest 
h 


5.  Thermal  conductivity  of  ferrites  follows  a  distinctive  temperature 
course  determined  by  the  state  of  magnetic  ordering.  Comparison  of  K 
demonstrates  a  certain  correlation  —  specimens  with  greater  I  are 
characterized  by  a  higher  value  of  K.  This  has  also  been  confirmed  by 
measurement  of  the  Curie  temperature: 


Q  -  IaNa>u> 


Specimens  marked  by  high  0R  correspond  to  a  higher  K  value. 

The  temperature  dependence  of  K  is  shown  in  Figure  5.  In  the  low- 

temperature  region  K  initially  rises  with  decreasing  T  owing  to  an  increase 

in  X  (T  180°K)  with  further  drop  in  T,  thermal  conductivity  falls 
c 

off  owing  to  effective  phonon  scattering  at  crystalline  defects,  since  they 
become  commensurable  with  the  size  of  defects.  In  the  positive  temperature 
region,  K  decreases  with  rising  T,  approximately  following  the  Eucken  law. 

In  the  0^  region  the  temperature  dependence  of  K  exhibits  the  maximum,  which 
is  caused  by  an  anomaly  of  C^.  It  can  be  concluded  from  measurements  of 
electroconductivity  and  6  that  CV*C, ■*<?„,  At  the  heat  capacity  exper¬ 
iences  a  jump  [9],  the  value  of  which 


A  C  «  1,5 


0K 

SI  o 


From  here  it  is  clear  that  the  maximum  on  the  temperature  dependence  of  K  is 

i 

proportional  to  1.  Thus,  curve  1  relating  to  zinc  ferrite,  whose  magnetic 
moment  equals  zero,  does  not  exhibit  a  K  maximum. 

The  maximum  on  the  curve  K  -  f(t)  for  the  ferrite  Ni^  2Fe2°4  is 

greater  than  for  the  ferrite  Nig  ^Zn^  gFe20^,  which  is  in  agreement  with 

their  I  and  0.  . 

k 


A 


I 
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The  smearing  out  of  maxima  on  the  experimental  curves  occurs  within  the 
Interval  50-80° K.  The  greater  the  1,  the  greater  the  smearing,  which  can  be 
accounted  for  by  fluctuations  of  spontaneous  magnetization.  Breakdown  of 
the  domain  structure  occurs  not  at  the  fixed  temperature  0^»  but  within  some 
temperature  range. 

When  ferrite  is  heated,  the  ordering  of  the  spinel  system  Is 
violated,  and  when  T  ■  0^  the  energy  of  thermal  agitation  becomes  comparable 
*•  with  the  energy  of  bulk  interaction  and  long-range  order  in  orientation  of 

spins  is  wholly  absent.  However,  owing  to  micros tructural  inhomogeneities 

« 

and  local  degaussing  fields  short-range  order  is  preserved.  Therefore,  when 
T  >  0^,  the  thermal  conductivity  after  going  through  the  maximum  does  not 
drop  off  sharply,  but  diminishes  smoothly  with  rise  in  temperature. 

Conclusions 

1.  Investigation  directly  of  the  phonon  component  of  thermal  conduct¬ 
ivity  and  thermal  expansion  of  ferrites  has  made  It  possible  to  establish 
their  interrelationship  with  other  physical  properties  of  ferrites  and,  iv 
particular,  with  magnetic  properties. 

2.  An  increase  in  the  mean  grain  size  brings  about  a  rise  in  thermal 
conductivity,  which  is  accounted  for  by  a  reduction  in  the  specific  volume  of 
defects  in  polycrystalline  ferrites. t 

3.  A  correlation  exists  between  thermal  conductivity  and  mean  atomic 
weight,  and  in  considering  this  correlation  the  nature  of  bonding  forces 

between  ions  must  be  borne  in  mind. 

A 

4.  The  values  of  coefficients  of  thermal  conductivity  and  of  thermal 

«  expansion  mutually  correlate  with  each  other. 

5.  Microhardness,  thermal  conductivity  and  thermal  expansion  of 
ferrites  exhibit  an  interrelationship  —  compositions  with  high  microhardness 
correspond  to  a  greater  thermal  conductivity,  and  those  with  low  microhard- 
ness  correspond  to  low  thermal  expansion. 

6.  Thermal  conductivity  of  ferrlteB  in  a  magnetic  field  rises,  which  is 
accounted  for  by  the  AE  effect. 
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7.  Compositions  with  high  magnetization  correspond  to  a  higher  thermal 
conductivity  coefficient.  Specimens  which  have  a  higher  Curie  temperature 
correspond  to  a  higher  thermal  conductivity  coefficient. 

Symbols 

K  -  coefficient  of  thermal  conductivity;  W  ■  coefficient  of  thermal 

expansion;  6  ■  coefficient  of  thermal  expansion;  ■  coefficient  of  thermal 

conductivity  in  the  magnetic  field  H;  H  -  Intensity  of  magnetic  field; 

I  -  magnetization;  I  -  saturation  magnetization;  I  -  magnetization  at 
8  0 

0°K;  Pg  ■  Bohr  magneton;  No  ■  number  of  atoms  per  gram-atom;  R  -  universal 
gas  constant;  h  »  Planck  constant;  k  -  Boltzmann's  constant;  A  »  mean 
atomic  weight;  D  -  density  of  mass;  C#  -  heat  capacity  of  crystal  lattice; 

-  heat  capacity  of  magnetic  ordering;  E  ■  Young's  modulus  in  the  field  H; 
Eq  «  Young's  modulus  in  the  absence  of  a  field. 
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Figure  1.  Coefficient  of  thermal  conductivity  of  ferrites  and 
coercive  force  as  functions  of  mean  grain  size:  Curves  1,  Ferrite 
2VT;  2,  MgFe204;  3,  NiFe^j  4,  H  of  ferrite  2VT 
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Figure  2.  Coefficient  of  thermal  conductivity  as  a  function  of 
mean  atomic  number  of  ferrite-chromites  and  ferrlte-aluminiteB 
with  equal  concentration  (y)  of  characterizing  ions: 

Curves  1,  y  ■  0. 9 j  2,  y  ■  0,6;  3,  y  ■  0.5;  4,  y  ■  0.2;  5,  Calculated 
values 
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Figure  3.  Reduced  thermal  resistance  and  thermal  expansion  as 
functions  of  crystal  lattice  constant 
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Figure  4.  Relative  change  of  thermal  conductivity  of  ferrites  In 
magnetic  field:  1 ,  MnFejO^S  2,  CoFe^O^ 
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Figure  5.  Temperature  dependence  of  thermal  conductivity  of 
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EXPERIMENTAL  STUDY  OF  THE  VISCOSITY  OF  THE  SYSTEM 

NjO^  £  ENOg  £  2NO  +  O2  in  the  Temperature  Range  300-780#K  and  the 

PRESSURE  RANGE  1  -SO  atm 

0.  V.  Belyayeva,  B.  G,  Maksimov,  V.  B.  Nesterenko, 

V.  N.  Pisarchik,  L.  A.  Pranovlch  and  B.  D.  Timofeyev 

Use  of  chemically  reacting  systems  of  the  type  NjO^  £  2N02  £  2N0  +  0^(*) 
to  intensify  heat  trannfer  [1],  cooling  of  nuclear  reactor  channels  [2],  and 
as  a  heat  carrier  in  a  gas-turbine  engine  [3]  requires  that  we  know  the 
thermophysical  properties  of  these  systems  in  a  wide  range  of  working 
temperatures  and  pressures.  The  present  paper  deals  with  an  experimental 
investigation  of  the  viscosity 'of  a  chemically  reacting  system  (*)  in  the 
temperature  range  300-780#K  and  the  pressure  range  1-50  atm  on  a  viscosimeter 
using  the  dropping  load  method,  the  design  of  which  and  methods  of 
viscosity  measurement  have  been  set  forth  in  detail  earlier  [4,  5].  Several 
load  dealgna  were  studied  for  use  of  the  viscosimeter  with  dropping  load  in 
inveetigating  the  viscosity  of  the  chemically  reacting  system  (*)  in  the 
indicated  range  of  temperatures  ahd  pressures.  These  designs  differ  from 
each  other  by  the  size  of  the  ring  gap  between  the  inner  diameter  of  the 
calibrated  channel  and  the  outer  diameter  of  the  load,  and  also  by  design 
elements  of  the  centering  devices  for  uniform  movement  of  the  load  over  the 
measured  section. 

Prior  to  and  after  the  experiments  on  the  system  (*),  the  viscosimeter 
was  weighed  in  nitrogen,  helium  and  carbon  dioxide  gas.  When  experimental 
data  on  viscosity  of  the  chemically  reacting  system  (*)  was  treated,  the 
following  expression  [5]  was  used: 
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where  the  coefficients  Aq  and  are  found  from  the  viacoelmeter  weighing 
data.  As  has  been  shown  earlier  [5],  the  coefficient  Bq  depends  on  the  gae 
flow  regime  in  the  ring  channel.  For  laminar  flow,  a  specific  function 
exists  for  each  load  with  specific  geometric  dimensions,  and  the  function  can 
be  represented  by  the  following  polynomial! 

*  Q°  4  ^  ’  (2) 


where 


Re  »  *jb(t  -to)]. 

The  coefficient  K  Is  quite  precisely  determined  based  on  geometric  dimensions 
of  the  load  and  the  calibrated  channel. 

Treatment  of  weighing  data  was  carried  out  on  the  Mlnsk-22  digital 
electronic  computer  at  the  Institute  of  Nuclear  Power  of  the  Academy  of 
Sciences  Belorussian  SSR. 

The  data  obtained  from  weighing  the  viscoslmoter  in  nitrogen,  helium  and 
carbon  dioxide  gas  were  used  to  determine  the  function  B 0  •  f(Re),  the 
coefficient  Aq,  and  it  was  also  established  that  in  experiments  on  load 
variant  IV  the  dimensionless  coefficient  of  load  friction  with  channel  wall 
f  «*  0  (Figure  1).  Here,  data  on  viscosity  and  the  thermodynamic  functions  of 
the  weighed  gases  necessary  for  the  calculation  were  taken  from  [6]  for 
nitrogen  and  helium,  and  from  [7]  for  carbon  dioxide  gas.  The  maximum  error 
of  individual  weighing  points  when  Bq  was  used  as  a  function  of  Re,  shown  in 
Figure  1,  and  the  coefficient  Aq  obtained  did  not  exceed  *5  percent,  while 
the  mean  deviation  of  112  weighing  points  in  absolute  value  was  1.9  percent. 
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Experimental  data  on  the  function  of  nitrogen  tetrcxlde  [8-10]  were  used 
in  treating  experimental  data  on  viscosity  of  the  system  (*),  and  in  the 
temperature  region  of  440®K  and  the  pressure  region  1-25  atm,  use  was  made 
of  data  calculated  according  to  the  method  of  the  Institute  of  Nuclear  Power 
of  the  Academy  of  Sciences  Belorussian  SSR  [11]. 

Figure  2  presents  graphs  of  experimental  data  on  the  viscosity  of  the 
chemically  reacting  system  (*)  obtained  in  the  Institute  of  Nuclear  Power  of 
the  Academy  of  Sciences  Belorussian  SSR,  on  the  viscosimeter  with  falling 
load  variant  IV,  Petker  and  Mason  [12]  used  a  viscosimeter  with  a  rolling 
sphere,  and  Beer  [13]  —  a  visccalmeter  with  a  spiral  capillary. 

It  is  clear  from  Figure  2  that  experimental  data  [12]  and  [13]  embrace 
the  temperature  range  only  from  3O0-443°K.  Satisfactory  agreement  of 
experimental  data  is  observed  only  at  room  temperature  and  atmospheric 
pressure.  At  pressures  of  1  and  5  atm,  experimental  data  of  the  Institute  of 
Nuclear  Power  of  the  Academy  of  Sciences  Belorussian  SSR  agreed  well  with 

[12]  throughout  the  entire  temperature  range  investigated.  Experimental  data 

[13]  at  the  maximum  temperature  of  396°K  attained  in  the  experiment  differed 
by  14.4  percent  from  the  experimental  data  of  the  Institute  of  Nuclear  Power 
of  the  Academy  of  Sciences  Belorussian  SSR  and  the  data  in  [12].  A 
reduction  of  the  experimental  data  obtained  in  [13]  obviously  must  be  made 
because  of  systematic  error  not  taken  into  account.  First  of  all,  the  error 
must  be  sought  in  imprecision  of  measuring  the  temperature  drop  in  the 
capillary  by  mercury  U-shaped  manometers,  the  mercury  which  was  in  contact 
with  nitrogen  tetroxlde  during  the  experiment.  Experiments  in  the  Institute 
of  Nuclear  Power  of  the  Academy  of  Sciences  Belorussian  SSR  have  shown 

that  during  the  period  of  contact  of  mercury  with  nitrogen  tetror.lde  when  the 
overall  impurity  content  (water  and  nitric  acid)  is  about  15  percent  salt 
formation  occurs,  which  for  all  practical  purposes  puts  the  mercury  U-shaped 
manometers  out  of  order. 

Conclusions 

1.  Experimental  data  obtained  at  the  Institute  of  Nuclear  Power  of  the 
Academy  of  Sciences  Belorussian  SSR  for  the  chemically  reacting  system 
NjO^  *  2N(>2  2N0  +  02  agree  well  with  data  in  [12]  at  pressures  of 


1  and  5  atm  In  the  temperature  range  investigated.  These  data  were  obtained 
for  the  load  variant  IV  with  gas-dynamic  stabilization  for  free  motion  in 
a  calibrated  channel. 

2.  The  mean  error  of  our  experimental  data  on  viscosity  of  the  system 

+■  2NC>2  2N0  +  0 2  was  estimated  at  2.5  percent;  the  maximum  error  of 

individual  points  did  not  exceed  5  percent  throughout  the  entire  range  of 
temperatures  and  pressures  investigated. 

In  conclusion  the  authors  express  their  gratitude  to  Academician  of  the 
Academy  of  Sciences  Belorussian  SSR  A.  K.  Krasin  for  constant  attention 
to  the  study  and  to  Professor  D.  L.  Tlmrot  for  his  kind  consultation  and 
interest  in  the  study. 

Symbols 

y  *  coefficient  of  dynamic  viscosity;  Ap  and  Bo  ■  coefficients  of 
viscosimeter;  t  -  dropping  time  of  load  for  uniform  motion  over  the  measured 
section;  pq  and  p  -  density  of  material  of  load  and  of  teat  medium,  , 
respectively;  f  •  dimensionless  coefficient  of  friction  of  load  against  wall 
of  calibrated  channel  when  the  channel  is  in  a  vertical  position; 
t  -  temperature;  g  -  coefficient  of  thermal  expansion  of  construction 
material  of  the  load  and  the  calibration  channel;  Re  ■  Reynolds  criterion 
characterizing  the  gas  flow  regime  in  the  circular  channel  between  the  load 
and  the  calibrated  tube. 
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EFFECT  OF  CRYSTALLINITY  ON  THERMOPHYSICAL  PROPERTIES  OF 
POLYTETRAFLUOROETHYLENE  (TEFLON-4) 

V.  S.  Bil' ,  Ye.  V.  Samardukov,  N.  Yu.  Gasteva 

and  R.  M.  Shchonova  , 

Polytetrafluoroethylene  (Teflon-4)  la  a  crystallizing  polymer,  the 
thermophysical  characteristics  of  which  have  not  been  closely  studied. 

Crystallizing  polymers  exhibit  the  same  kind  of  temperature  dependencies 
of  thermal  diffusivity,  heat  capacity  and  density,  and  different  kinds  of 
temperature  dependencies  of  thermal  conductivity.  Crystallization  induces  a 
rise  in  thermal  conductivity,  thermal  diffusivity,  density  and  the  reduction 
in  heat  capacity.  Thus,  for  example,  from  the  data  in  [1]  for  polyethylene 
the  thermal  conductivity,  like  the  low-molecular  characteristics  of 
corpounds,  reduces  with  temperature  rise,  but  for  polytetrafluoroethylene, 
like  low-molecular  amorphous  compounds,  thermal  conductivity  rises  with  rise 
in  temperature. 

The  effect  of  crystallinity  on  therraophysical  properties  of  polytetra¬ 
fluoroethylene  have  been  studied  with  an  instrument  described  in  [2]  using 
a  method  given  in  [3]  for  the  temperature  range  40-30Q°C.  The  density  was 
determined  by  the  hydrostatic  weighing  method  for  t  <  20“C. 

Figure  1  presents  experimental  data  on  thermal  conductivity  of 

polytetrafluoroethylene  as  a  function  of  temperature  for  two  density  values 

3  3 

(crystallinity)  of  2140  kg/m  and  2180  kg/m  .  (Scatter  of  experimental 
points  *4  percent,) 

According  to  the  data  in  [4]  for  polytetrafluoroethylene,  density  is 
unambiguously  related  to  crystallinity  by  the  relationship 


'  S  "2.0  , 


(i) 


-486- 


» 


This  affords  the  possibility  by  using  experimental  data  on  thermal  conduct¬ 
ivity  to  extrapolate  temperature  dependencies  of  thermal  conductivity  for 
density  (crystallinity)  values  corresponding  to  amorphous  and  crystalline 
phases  with  the  aid  of  the  two-phase  model  proposed  by  I.  Eiermann  [1]. 

Under  the  two-phase  model  it  is  assumed  that  a  crystalline  polymer  Is  a 
homogeneous  mixture  consisting  of  crystalline  and  amorphous  regions,  and 
therefore  thermal  conductivity  can  be  calculated  as  the  thermal  conductivity 
of  the  mixture  following  a  formula  presented  in  [1] 


The  temperature  dependencies  of  X,  and  X  were  calculated  from  formula  (2) 
using  temperature  dependencies  of  X(t)  that  we  obtained  for  two  dendity 
values.  The  results  sre  Indicated  in  Figure  1  by  solid  lines. 

The  results  of  measuring  thermal  diffusivlty  of  polytetraf luoroethylene 
3  3  3 

at  densities  of  2100  kg/m  ,  2118  kg/m  ,  and  2160  kg/m  shown  in  Figure  2 
demonstrate  that  within  the  limits  of  experimental  precision  the  coefficient 
of  thermal  diffusivlty  depends  linearly  on  temperature  within  the  range 
investigated.  Additionally,  the  higher  the  crystallinity  (density),  the 
higher  the  absolute  value  and  the  higher  the  rate  of  decrease  of  thermal 
diffusivlty  with  temperature  rise.  Experimental  points  can  be  described  by 
the  following  equation  with  precision  of  up  to  *2  percent: 

a  *  {  Q685  +  4, 2  ( o  -  2000)'10'*- 

k  :  2  (3) 

- [0.110  ♦  10,650(53  - 2000)'  10*‘](T#C  - 273)- 10*3}  ’ lO*’  ^—-1 

Equation  (3)  allows  us  to  calculate  the  temperature  dependencies  of 
coefficients  of  thermal  diffusivlty  of  crystalline  and  amorphous  phases,  and 
also  the  function  a(T)  for  intermediate  crystallinity  values  (Figure  2). 
(Calculated  curves  are  indicated  by  solid  lines.) 

Thus,  the  experimental  and  calculated  data  obtained  allow  us  by  using 
density  (crystallinity)  of  polytetrafluoroethylene  to  estimate  thermal 
conductivity  and  thermal  diffusivlty  as  function  of  temperature. 
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Symbols 

X  ■  thermal  conductivity  of  polymer,  w/m*deg;  X^  and  Xfl  -  thermal 

conductivity  of  Its  crystalline  and  amorphous  regions;  a  ■  coefficient  of 

2  3 

thermal  diffuslvlty,  m  /sec;  p  -  density,  kg/m  ;  W  -  bulk  extent  of  crystal¬ 
linity;  T  -  temperature  (°K);  t  ■  temperature  (#C). 


Figure  1.  Thermal  conductivity  of  polytetrafluoroethylene  as  a 

function  of  temperature  for  different  densities  (crystallinities) : 

3  1 

1.  P  ■  2300  kg/m  ,  wholly  crystalline;  2,  p  -  2250  kg/m  , 

83  percent  crystallinity;  3,  p  -  2200  kg/m3,  66  percent  crystal- 

Unity;  4,  p  -  2180  kg/m  ,  60  percent  crystallinity; 

5,  p  -  2140  kg/m3,  46.5  percent  crystallinity;  6,  p  -  2100  kg/m3, 

33  percent  crystallinity;  7,  p  -  2050  kg/m3,  16.5  percent  crystal- 

3 

Unity;  8,  p  ■  2000  kg/ra  ,  wholly  amorphous 
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Figure  2.  Thermal  dlffuslvlty  of  polytetrafluoroethylene  as  a 

function  of  temperature  at  different  densities  (crystallinities) : 

3  3 

1,  p  -  2300  kg/m  ,  wholly  crystalline;  2,  p  ■  2250  kg/m  , 

3 

83  percent  crystallinity;  3,  p  ■  2200  kg/m  ,  66  percent  crystal- 

3 

Unity;  4,  p  »  2160  kg/m  ,  53  percent  crystallinity; 

3  3 

5,  p  ■  2118  kg/m  ,  39.2  percent  crystallinity;  6,  p  ■  2100  kg/m  , 

3 

33  percent  cr-atallinity;  7,  p  ■  2050  kg/m  ,  16.5  percent 

3 

crystallinity;  8,  p  -  2000  kg/m  ,  wholly  amorphous 

I  I 
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THERMAL  CONDUCTIVITY  OF  THE  SYSTEM  NjO^  t  2N02  *  2N0  +  02  IN  THE 
TEMPERATURE  RANGE  300-780#K  AND  THE  PRESSURE  RANGE  1-50  atm  abs 

A.  A.  Bilyk,  F.  N.  Dresvyannikov,  Yu.  G.  Kotelevskly, 

L.  V.  Mishina,  V.  B.  Nesterenko  and  B.  D.  Timofeyev 

Several  studies  have  dealt  with  experimental  Investigation  of  thermal 
conductivity  of  nitrogen  tetroxide  [1-3].  However,  all  these  studies  used 
the  heated  filament  method  at  low  temperature  (up  to  1  atm  abs) . 

The  present  study  sets  forth  the  results  of  measurements  on  the  thermal 
conductivity  of  the  systemNjO^  +■  2NO,  +  2N0  +  02  at  pressures  from  1  to 
50  atm  abs  in  the  temperature  range  300-7 80®K. 

A  method  of  coaxial  cylinders  using  stainless  steel  of  the  grade 
KhlSNIOT  as  the  main  construction  material  was  selected  In  studying  thermal 
conductivity  at  the  Institute  of  Nuclear  Power  of  the  Academy  of  Sciences 
Belorussian  SSR.  This  material  satisfies  to  the  greatest  extent  require¬ 
ments  posed  In  investigating  chemically  reacting  gases. 

The  experimental  device  consisted  of  the  following  assemblies:  a)  a 
working  thermostat  serving  to  thermostat  the  experimental  section;  b)  the 
experimental  section;  c)  a  system  for  control  of  pressure  in  the  working 
cavity,  housing  a  model  manometer  with  a  membrane  divider  using  carbon 
tetrachloride  as  the  dividing  liquid;  d)  a  vacuum  system  including  a 
forevacuum  pump  and  a  vacuum  gauge;  e)  a  system  for  filling  the  working 
cavity  with  the  test  gas,  consisting  of  a  container  filled  with  condensed 
nitrogen  tetroxide  and  a  thermostat  serving  to  heat  It. 

The  working  thermostat  was  a  structure  thermally  insulated  from  all 
sides  containing  a  thick-walled  copper  cylinder  within,  over  whose  external 
surfaces  a  four-section  heater  made  of  nichrome  wire  was  wound.  Activation 
of  the  thermostat  to  the  desired  regime  and  maintenance  in  it  of  a  constant 
temperature  was  done  automatically  with  the  aid  of  an  electronic  temperature 
control. 


The  experimental  section  (Figure  1)  consisted  of  two  coaxially  arranged 
cylinders  made  of  stainless  steel.  The  outer  diameter  of  the  Innermost 
cylinder  was  8  mm  and  the  wall  thickness  was  1  mm.  An  approximately  0.8-mm 
gap  between  the  cylinders  was  maintained  with  a  ceramic  washer  (1)  placed  on 
the  side  of  the  cylinder  dead-ends  and  a  thin  metal  washer  (3)  on  the  side  of 
flange  joining  of  the  cylinders ,  which  together  with  the  Teflon  inserts  (4) 
ensured  hermetization  of  the  working  cavity  within  the  cylinders.  Temper¬ 
ature  was  measured  at  the  inner  surface  of  the  innermost  cylinder  with  the 
aid  of  three  spring-mounted  chromel-alumel  thermocouples  (7)  placed  50  mm 
from  each  other. 

A  rosette  of  five  thermocouples  (2)  was  soldered  onto  the.  outer  surface 
of  the  outermost  cylinder  using  condenser  walding  at  a  spacing  of  40  mm  from 
each  other  to  record  the  temperature  field  along  the  cylinder  length. 

A  regulated  heater  (5)  waa  provided  to  reduce  heat  losses  from  the 
flange  side  in  the  lower  section  of  the  experimental  area. 

The  heat  flux  was  generated  by  means  of  the  three-section  heater  (6) 
made  of  0.23-mm-dlameter  nichrome  enclosed  In  glass  tubing  of  about  5  mm 
diameter  and  placed  in  the  innermost  cylinder.  To  record  the  power  supplied 
over  100  mm  of  length,  potential  leads  of  copper  wire  were  connected  from  the 
central  section  of  the  innermost  heater. 

The  inner  heater  was  fed  from  a  storage  battery.  Measurements  of  the 
emf  of  thermocouples  and  voltages  over  the  working  section  and  the  model  coil 
in  the  inner  heater  circuit  were  made  on  the  R-307  high-ohmic  potentiometer. 

Choice  of  the  geometry  of  the  experimental  section  made  it  possible  to 
exclude  convection  over  the  entire  range  of  temperatures  and  pressures 
(Gr*Pr  <  1000). 

Use  of  the  relative  method  of  coaxial  cylinders  made  it  possible  to  take 
into  account  correction  for  noncoaxiality  of  cylinders  when  weighing  for  a 
known  matei-ial.  This  correction  was  included  in  the  apparatus  constant. 

Since  the  literature  contains  no  data  on  the  reduced  degree  of  blackness 
v  for  stainless  steel  passivated  by  nitrogen  tetroxide,  the  correction 
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for  radiation  was  provided  for  by  calculation.  It  was  determined  thusly. 

The  overall  quantity  of  heat  Q  transferred  from  the  inner  cylinder  to 
the  outer. 


Q 


X  aT 

— i —  » 


(l) 


where  A  is  the  apparatus  constant,  determined  by  the  geometry  of  the  section. 


The  amount  of  heat  transferred  by  radiation  was  determined  from  the 
formula 


(2) 


where 

8  *  grp 

(3) 

For  the  case  when  the  value  of  A  did  not  depend  on  temperature,  the 
coefficients  A  and  B  can  be  determined  from  the  following  scheme.  Since  at 
temperature  drops  up  to  5°K  at  low  temperatures  (up  to  373BK),  even  for 
oxidized  steel  which  has  the  highest  degree  of  surface  blackness,  the 
contribution  of  radiation  was  less  than  one  percent  of  the  total  amount  of 
heat.  Therefore  we  can  neglect  radiation  and  maintain  adequate  precision, 
and  determine  the  value  of  A  from  weighing  experiments  with  known  material. 

At  high  temperatures  we  can  calculate  B  ■  B(T)  for  weighed  gas  using  the 
value  of  A  obtained  at  low  temperatures. 

The  value  of  A  was  calculated  at  several  points  in  the  temperature  range 
296-356*K  in  the  weighing  experiments  with  air.  A  as  a  function  of  temper¬ 
ature  was  not  established  in  our  experiments. 

The  value  of  B  was  also  determined  by  weighing  experiments  in  air,  and  B 
as  a  well  defined  function  of  temperature  was  not  obtained.  A  scatter  about 
the  mean  value  was  observed. 
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5.27  were  used  in 


Thus,  the  data  obtained  for  A  =  0.236  and  B*10 
treating  measured  coefficients  of  the  thermal  conductivity  of  nitrogen 
tetroxide.  To  verify  the  correctness  of  the  determinations  of  A  and  B, 
control  experiments  with  helium  were  conducted.  Results  for  X  gave  satis-' 
factory  agreement  with  known  literature  data  on  helium  [4], 

The  measurements  made  of  thermal  conductivity  of  dissociating  nitrogen 
tetroxide  in  the  temperature  range  300-780*K  at  pressures  of  1,  10,  30  and 
50  atm  abs  are  shown  in  Figures  2  and  3. 

The  mean  error  of  experimental  data  was  estimated  at  6  percent,  and  the 
maximum  error  at  9  percent. 

Experimental  data  for  1  atm  (cf.  Figure  2)  at  temperatures  up  to 
~  450°K  (first  stage  of  the  reaction)  found  quite  good  agreement  with  the 
data  of  other  authors  [1,  3]  and  with  calculated  data  made  on  the  assumption 
of  local  chemical  equilibrium.  The  maximum  departure  from  calculated  values 
was  approximately  15  percent.  No  contribution  of  chemical  reaction  to 
thermal  conductivity  was  noted  in  the  second  stage  in  the  temperature  range 
450-780°K,  and  the  experimental  points  fitted  on  the  curve  of  the  "frozen" 
component. 

Experimental  values  of  X  at  a  pressure  of  50  atm  abs  in  the  temperature 
range  550-780°K  are  at  a  level  approximately  20  percent  higher  than  at 
1  atm  abs  (Figure  3).  The  influence  of  pressure  alone  cannot  account  for 
this  increase,  but  it  is  obviously  necessary  to  relate  it  to  the  contribution 
of  the  chemical  reaction. 

The  authors  express  their  gratitude  to  Academician  of  the  Academy  of 
Sciences  Belorussian  SSR  A.  K.  Krasin  for  constant  interest  in  the  study. 
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Figure  1.  Diagram  of  experimental  section.  1,  Ceramic  washer; 

2,  Thermocouples;  3,  Metal  washer;  Teflon  Inserts;  5,  Regulated 
heater;  6,  Three-section  heater;  7,  Spring-mounted  thermocouples 
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Figure  2.  Temperature  dependence  of  coefficient  of  thermal 
conductivity  of  nitrogen  tetroxide  at  P  -  1  atm  aba:  1,  Calculated 
with  the  contribution  of  the  chemical  reaction  taken  into  account; 
2,  Calculated  for  the  mixture  of  "frozen"  composition;  3,  Data  in 
[1];  A,  Data  in  [3];  5,  Authors'  data 


Mffi 


Figure  3.  Change  of  thermal  conductivity  coefficient  for  nitrogen 
tetroxide  as  a  function  of  temperature  and  pressure: 

1,  P  -  1  atm  abs;  2,  P  -  10  atm  abs;  3,  P  -  30  atm  ate; 

A,  P  -  50  atm  abs 
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THERMAL  CONDUCTIVITY  OF  POLYETHYLENE  AS  A  FUNCTION  OF  ABSORBED 

RADIATION  DOSE 

B.  A.  Brlskman  and  V.  D.  Bondarev 

Data  on  thermal  conductivity  X  of  polyethylene  as  a  function  of  temper¬ 
ature  and  the  value  of  absorbed  radiation  energy  are  of  interest  from  two 
points  of  view.  First,  the  behavior  of  polyethylene  as  an  insulator  is 
important  in  applied  problems.  Second,  study  of  the  X  of  polyethylene 
provides  additional  information  on  the  structure  of  this  polymer  and  on 
physicomechanlcal  and  chemical  transformations  occurring  when  it  is 
irradiated. 

The  X  of  both  high-pressure  and  low-pressure  polyethylene  as  a  function 
of  temperature  has  often  been  investigated,  for  example  in  the  temperature 
range  -190°  to  +110°-120°C,  that  is,  practically  up  to  the  melting  point  of 
uncross-linked  polyethylene  [1,  2].  It  has  been  shown  that  throughout  almost 
the  entire  temperature  range,  with  the  exception  of  very  low  temperatures, 
the  X  of  polyethylene  drops  otf  monotonically  with  rise  in  temperature,  while 
the  X  of  high-density  polyethylene  (low-pressure)  is  considerably  higher  than 
the  X  of  low-density  polyethylene  (high-pressure). 

As  we  know,  polyethylene  is  a  two-phase  polymer,  a  combination  of 
crystalline  and  amorphous  phases.  Bearing  in  mind  that:  1)  the  X  of  the 
crystalline  phase  is  considerably  higher  than  the  amorphous  phase;  2)  the  X 
of  the  crystalline  phase  drops  off  appreciably  with  temperature  rise;  3)  the 
X  of  the  amorphous  phase  rises,  but  more  slowly;  4)  the  extent  of  crystal¬ 
linity  drops  off  [3]  with  temperature  rise  all  the  way  to .complete  failure  of 
crystals  in  the  melting  point  region,  with  the  foregoing  in  mind  it  can  be 
shown  that  the  sign  of  the  temperature  effect  for  the  X  of  polyethylene 
depends  on  the  initial  phase  composition  of  polyethylene,  the  rate  of 
degradation  of  the  crystalline  phase,  and  on  the  temperature  coefficients  of 
the  c'vstalline  and  amorphous  phases. 
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Directing  attention  to  «.«mui.~eLc  values  of  (.lie  quanLii.ieu  listed  above, 
it  must  be  concluded  that  the  decrease  in  the  X  of  polyethylene  with 
temperature  rise  is  due  to  a  reduction  in  the  degree  of  crystallinity.  The 
mechanism  of  variation  of  X  of  polyethylene  exposed  to  radiation  is  more 
complex.  The  results  of  a  study  of  X  of  polyethylene  as  a  function  of 
absorbed  radiation  energy  have  been  presented  in  [4].  Under  the  effect  of 
radiation  [5],  radiation  cross-linking  of  the  molecular  polyethylene  chains 
takes  place,  which  leads  to  a  reduction  in  distances  between  carbon  atoms  and 
to  higher  polyethylene  density.  At  the  same  time,  a  decrease  in  extent  of 
crystallinity  occurs  [5] ,  leading  to  a  drop  in  polyethylene  density.  As  haB 
been  shown  in  [4],  the  X  of  polyethylene  drops  off  with  rise  in  absorbed 
radiation  energy  at  a  temperature  of  approximately  20°C,  but  at  relatively 
higher  temperatures  (approximately  150°C)  it  rises.  We  must  note  that  the 
absolute  values  of  polyethylene  (below  we  will  discuss  low-density  poly¬ 
ethylene)  vary  from  0.26  to  0.47  w/m*°C  (for  t  ■  20°C).  Such  a  wide  range  of 
polyethylene  X  values  is  accounted  for  by  the  dependence  of  its  physico- 
mechanlcal  properties  on  the  method  of  preparation,  storage,  etc.  For 
irradiated  polyethylene  account  must  be  taken  (if  the  irradiation  is 
conducted  in  air)  of  the  presence  of  oxygen  compounds  (carbonyl  groups, 
especially  on  the  specimen  surface,  and  peroxide  bridges  between  molecular 
chains).  It  is  very  important  to  allow  for  the  action  of  gaseous  hydrogen 
formed  upon  irradiation.  To  grasp  the  mechanism  of  thermal  conductivity  we 
must  bear  in  mind  also  the  presence  of  thermal  conductivity  anisotropy  [6]. 

The  aim  of  this  study  is  obtaining  of  additional  Information  on  the  X  of 

low-density  irradiated  polyethylene  (y  ■  0.920  g/cm  ).  Unatabillzed  block 

polyethylene  in  the  form  of  disks  of  30  mm  diameter  and  15  mm  thickness  was 

irradiated  in  a  nuclear  reactor  in  air  and  was  used  in  the  same  form  in  the 

experimental  apparatus  (Figure  1).  The  absence  of  impurities  in  polyethylene 

greater  than  10-^  percent  by  weight  was  corroborated  by  activation  analysis. 

The  measurements  were  made  by  a  steady-state  method  due  to  the  appreciable 

dependence  of  density  and  heat  capacity  of  polyethylene  on  temperature.  The 

-4 

experimental  apparatus  was  a  vacuum-treated  chamber  (pressure  of  10  - 

10" 5  mm  Hg)  in  which  the  specimens  (2)  were  placed  with  the  heater  (5) 
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located  between  them.  Symmetrical  removal  of  heat  was  achieved  by  silver- 
plated  massive  copper  refrigerators  (1).  Reduction  of  heat  losses  via 
radiative  heat  transfer  was  achieved  with  the  aid  of  ten  cylindrical  screens 
(4)  made  of  polished  stainless  steel  foil  (6  ■  0.1  mm).  Chromel-copel 
thermocouples  (three  for  each  spcelmen)  are  embedded  along  the  specimen 
axis.  The  spacing  between  the  hot  thermocouple  arms  is  measured  by  X  rays. 
The  measurements  were  made  in  the  temperature  range  25-155°C  and  the 
absorbed  dose  range  0-2300  Mrad.  Thermal  losses  did  not  exceed  0.5  percent 
of  heater  power.  The  maximum  experimental  error  was  *3.5  percent.  Results 
of  measurements  are  shown  in  Figures  2  ‘and  3. 

As  can  be  seen  from  Figure  2,  at  a  temperature  of  30#C  the  X  of  poly¬ 
ethylene  drops  off  sharply  with  rise  in  dose,  which  is  accounted  for  by  the 
decrease  in  extent  of  crystallinity.  Beginning  with  some  dose  value, 
depending  on  temperature,  X  is  reduced  considerably  more  slowly,  this  1b 
also  the  case  fcr  the  variation  of  X  with  dose  for  the  temperatures  40,  60 
and  80°C.  At  temperatures  of  100  and  200°C,  X  drops  off  somewhat  up  to  a 
dose  of  600  Mrad,  which  is  probably  explained  by  the  failure  of  the  laBt 
crystals.  Then  X  monotonically  rises,  revealing  a  trend  toward  stabiliza¬ 
tion.  At  a  temperature  of  150eC,  X  rises  throughout  the  dose  range.  We 
explain  the  course  of  all  these  curves  by  proceeding  from  the  following 
considerations.  Two  competing  processes  are  operative  with  increase  in 
dose  —  radiation  cross-linking,  leading  to  rise  in  X,  and  disruption  of 
polyethylene  macrostructure  owing  to  gas  release  upon  irradiation.  In  gas 
release  fine  pores  appear,  the  size  and  number  of  which  rise  with  increase  in 
received  dose  [7].  The  presence  of  vacuumed  pores  and  fissures  leads, 
especially  at  moderate  temperatures,  to  worsening  of  polyethylene  X.  The 
thermal  resistance  of  pores  drops  off  with  rising  temperature  owing  to  the 
ever  mounting  role  of  radiative  heat  transfer,  which  brings  about  a  rise  in  X 
with  increase  in  dose  at  relatively  high  temperatures. 

As  we  can  see  from  Figure  3,  X  as  a  function  of  temperature  is  fairly 
complex,  where  at  temperatures  up  to  50°C  clearly  pronounced  "valleys"  of  X 
occur,  the  reliability  of  which  is  confirmed  by  the  good  reproducibility 
of  the  experimental  points.  At  doses  of  0  and  100  Mrad,  X  drops  off  owing 
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to  the  reduced  degree  of  crystallinity.  At  a  dose  of  900  Mrad,  the  crystals 
disappear  already  at  a  temperature  of  40°C  and  thereafter  X  increases.  At  a 
dose  of  2300  Mrad  polyethylene  from  the  very  outset  is  an  amorphous  phase, 
whose  thermal  conductivity  rises  with  temperature.  Throughout  the  range  of 
absorbed  doses  there  is  a  correlation  with  the  specific  volume  of  poly¬ 
ethylene. 

In  the  future  it  will  be  of  great  Interest  to  study  the  effect  of 
irradiation  conditions  on  polyethylene  thermal  conductivity. 

We  are  duty  bound  to  express  our  gratitude  to  Professor  and  Doctor  of 
Chemical  Sciences  V.  L.  Karpov  for  participation  in  discussion  and  for 
valuable  counsel. 
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Figure  1,  Sketch  of  experimental  apparatus:  1,  Copper  refriger¬ 
ators;  2,  Polyethylene  specimen;  3,  Chromel-copel  thermocouple; 

4,  Screen;  5,  Heater 
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Figure  2.  Thermal  conductivity  of  polyethylene  as  a  function 
of  absorbed  radiation  dose 


Figure  3.  Thermal  conductivity  of  polyethylene  as  a  function 
of  temperature 
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TOUARD  A  STATISTICAL  THEORY  OF  THERMAL  CONDUCTIVITY  OF  CONDENSED  MEDIA 
V.  3.  Vikhrenko,  V.  B.  Nemtsov  and  L.  A.  Rott 

Investigation  of  the  thermal  conductivity  of  condensed  systems  is 
complex  owing  to  the  fact  that  there  are  certain  difficulties  in  describing 
the  perturbation  of  the  Hamiltonian  of  the  system,  owing  to  the  thermodynamic 
nature  of  the  perturbations. 

The  well  known  Kirkwood  theory  [1]  based  on  a  statistical  interpretation 
of  the  theory  of  Brownian  motion  males  it  possible  to  describe  the  thermal 
conductivity  phenomenon  with  satisfactory  rigor.  But  the  actual  calculation 
of  the  coefficient  of  thermal  conductivity  is  hampered  to  a  large  extent 
owing  to  the  absence  of  a  well  defined  expression  for  the  binary 
distribution  function. 

In  the  present  study  the  Kirkwood  method  is  applied  in  a  statistical 
scheme  of  conventional  distributions,  which  has  been  developed  earlier  (cf., 
for  example,  [2]). 

In  essence,  the  method  of  conventional  distributions  amounts  to  the 
following.  A  volume  V  occupied  by  system  of  N  particles  is  divided  into 

cells  of  volume  v  ■  V/N.  Then  the  probability  density  p,  p* . ,t) 

that  at  the  moment  of  time  t  t  particles  will  be  found  in  a  certain  specific 
cell  v, ,  k  particles  —  in  some  other  cell,  l  particles  in  a  third  cell,  and 
so  on  in  the  vicinity  of  specific  coordinates  q  and  q5  is  studied,  ....  these 
particles  exhibiting  the  momenta  p,  p  ,  ...  under  the  condition  that  in  the 
remaining  cells  the  remaining  N-i-k-J-...  molecules  are  distributed 
arbitrarily,  but  that  no  more  than  j  molecules  occupy  a  cell.  Here  the 
requirement  i,  k,  ?,  ...  <  J  is  satisfied. 
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An  aCtraccive  feature  of  the  method  applied  In  the  possibility  of  using 
successive  approximations.  For  example,  in  the  so-called  approximation 

employed  here,  only  those  system  states  are  taken  into  account  in  which 
each  cell  contains  one  molecule.  For  condensed  media  the  molecular  volume  v 
is  relatively  small  and  even  the  approximation  provides  satisfactory 
results  [4], 


Use  of  higher  approximations,  for  example,  F^2  and  F22  in  which  the 
system  states  are  taken  into  account  where  each  cell  can  simultaneously  have 
present  two  molecules  allows  us  even  more  precisely  to  describe  the  behavior 
of  the  system. 


Hydrodynamic  laws  of  conservation  were  obtained  by  using  the  method  of 
conventional  distributions  in  the  F  ^  approximation  as  well  as  an  expression 
of  the  heat  flux  density,  which  follows  from  the  equation  of  the  conservation 
of  energy, 


(1) 


Here  the  kinetic  component  of  the  heat  flux  density 


and  the  potential  component 


9  +  Kp,Pit)clpdjSt<tfi 


(3) 


where  r  *  q  -  q;  u  is  the  mean  molecular  velocity.  The  expression 

>  1  k  -V  + 

rr  represents  a  diad,  while  the  expression  r*r  is  the  scalar  product; 
m  is  the  mass  of  molecule,  t(r)  is  the  potential  of  interaction  between 
molecules.  The  interaction  is  presumed  to  be  central  and  pairwise. 
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Following  the  concept  of  Kirkwood  [5],  we  Introduce  distribution 
functions  smoothed  over  the  rime  r; 


(4) 


The  time  Interval  t  over  which  distribution  functions  are  smoothed  must 
be  large  In  comparison  with  the  relaxation  time  of  lntermolecular  forces,  but 
much  less  than  the  Poincare  cycle  of  the  system. 

Let  represent  the  Hamiltonian  of  a  system  of  N  particles.  Then  the 
distribution  function  of  the  dynamic  states  of  the  entire  system 

satisfies  the  Llouville  equation 


where  the  righthand  member  Is  the  Poisson  bracket. 


Averaging  equation  (5)  over  the  time  interval  t  and  integrating  it  with 
respect  to  coordinates  and  momenta,  as  shown  In  [2],  we  arrive  at  equations 
determining  the  smoothed  distribution  functions  %  (tf,  p it)  and 
in  the  form: 


(6) 
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*“1 ‘fltltlt 


2W$,  a  d$"Ui 


HI  SI  is  determined  in  a  similar  way.  Here  l  denotes  summation  over  all 

(p) 

possible  rearrangements  of  particles  in  cells  over  which  the  Integration  is 
performed.  fy)*-  is  tlie  *otce  acting  on  the  i-th  molecule 

on  the  part  of  the  j-tli  molecule. 

The  problem  consists  in  transforming  the  righthand  members  of  equations 
(6).  Assuming  that  the  vicinity  of  two  molecules  is  readily  deviated  from 
the  equilibrium  state  and  adopting  the  Kirkwood  method  in  [5]  for  the  first 
of  the  equations  (6),  we  obtain  an  equation  determining  the  function 

.  £-(§,  Pit). 


f  .|-f *$•/*  W*/ 


Here  the  parentheses  denote  averaging  over  the  equilibrium  state 

M 

of  all  particles,  with  the  exception  of  one.  K,  -  is  the  total  force 
acting  on  the  molecule,  *  is  nonequilibrium  correction  stemming  from 
the  deviation  of  the  vicinity  of  two  molecules  from  the  equilibrium  state. 

The  scalar  coefficient  of  friction 


3kT  J^C %(v '  S 


The  resulting  equation  (7)  agrees  with  the  Chandrasekar  equation  [6] 
well  known  in  the  theory  of  Brownian  motion.  However,  equation  (7)  was 
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\ 


obtained  by  a  statistical  method  «nd  the  coefficient  of  fraction  E, 

Is  expressed  here  via  parameters  of  intermolecular  interaction,  while  in  the 
Chandrasekar  equation  the  coefficient  of  friction  is  a  phenomenological 
constant. 

Assuming  that  the  vicinity  of  three  molecules  is  easily  deviated  from 
the  equilibrium  state,  and  similarly  considering  two  of  equations  (6),  we 
find  an  equation  determining  the  function 

|  A  q*  p‘;  t). 


The  bent  brackets  here  again  denote  statistical  averaging  oyer  the 
equilibrium  state  of  all  particles,  except  for  two. 

The  forces  and  are  due  to  deviation  of  the  state  of  the 

three-molecule  Vicinity  from  the  equilibrium  state. 

Let  ua  determine  the  kinetic  contribution  to  the  heat  flux  density.  For 
this  purpose,  we  will  multiply  equation  (7)  by  +  and  integrated  over 
the  momentum  space.  As  a  result  we  have: 

/A# 

where  the  mass  flux  density 
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4 


A  am /■&  ?»(§ .tit) dp. 


and  the  density  of  the  material 


M&O'm  fp.(q,  ft t) dp 


We  will  introduce  the  determination  of  temperature  according  to  the 
equality 


(§.&  wfi'irwYs 


where  I  is  a  unit  tensor. 

Comparing  (10)  and  (12)  and  neglecting  the  inertia  members,  and  also  the 
squared  members  in  terms  of  mean  velocity  u,  after  several  transformations  we 
find: 


Using  (13)  in  (10)  and  again  neglecting  the  inertia  members  and  the 
second  and  higher  order  members  in  terms  of  mean  velocity  u,  we  similarly 
obtain  an  expression  for  the  kinetic  component  of  the  heat  flux  density 

.  .  .  :  •  ■  (i4) 

- 

To  determine  the  contribution  of  interaction  between  molecules  to  the 
heat  flux  density,  we  integrate  the  expression  (3)  over  the  momentum  space 
and  take  into  account  that  <f'(r)  differs  from  zero  in  the  region  of  the 
order  of  intermolecular  distances.  Further,  we  write  the  expression  for  the 
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Using  equation  (9)  and  making  calculations  similar  to  those  in  [5],  we 
find  the  function: 


(19) 


Here  we  take  note  of  the  fact  that  in  homogeneous  media  we  can  assume 


Using  (19)  in  the  expression  for  the  potential  component  of  the  heat 
flux  density  (15),  we  have: 

{ jf<PM  -  >0>  V.  ** 

-  Jfi<p{rj - r<»,(i'J]rt( jfflct rj  v7 j  v) 

Using  the  Fourier  phenomenological  law  of  thermal  conductivity 

7*  -KrVT, 

we  find  an  expression  for  the  coefficient  of  thermal  conductivity  (according 
to  (14)  and  (20)) 

Kr  *  K*un  *Kmm. 

K™  mjjj$  j“  f[ Wrj-rV'(n/]nY(v£( (22) 

/V 

+  JcAVKrJ-rvWJr'fffiJdrJ 


(20) 


(21) 
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The  function  can  be  approximated  by  an  expression  in  [ 3 ] 


^  (rj « A  exp  (' 


(23) 


where  a  is  some  small  parameter  and  A  is  a  normalizing  multiplier,  determined 
by  the  normalizing  condition 


jr'/Ve  (r)dfi 
a 


(24) 


Adopting  the  approximation  (23)  and  integrating  the  first  term  in  the 
potential  part  of  the  thermal  conductivity  coefficient  by  parts,  we  arrive  at 
the  expression: 

jf$i  {J[rq>‘(rj.  <P(r)]na 

)  (25) 

-/ [r  H(r)  *  V ro>'(pJ-  60>Crj]na(j£)ol rf 
n , 


Here  the  derivative  with  respect  to  temperature  is  adopted  either  at  constant 
volume  or  at  constant  pressure,  depending  on  whether  we  need  to  know  the 
thermal  conductivity  at  constant  volume  or  at  constant  pressure. 

The  coefficient  of  friction  for  simple  nonpolar  liquids  has  been 
previously  determined  in  f  7  j 


2  ■><>•*  R.eec-i 


(26) 


where  e/k  and  a  are  parameters  of  the  Lennard-Joncs  potential  (o  ■  Angstroms, 
c/k  in  degrees,  and  M  -  molecular  weight). 

In  higher  approximations  (F^  and  F22^  calculetion  of.  the  thermal 
conductivity  coefficient  will  be  considered  at  a  later  date. 
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THERMOPHYSICAL  CHARACTERISTICS  OF  BINARY  SOLUTIONS  AS  A  FUNCTION 

OF  CONCENTRATION 


G.  V.  Grishchenko 

In  recent:  years  a  large  number  of  diverse  methods  for  determining 
thermophysical  characteristics  of  solids  have  appeared.  But  only  some  of 
these  have  been  used  for  determining  thermophysical  characteristics  of 
liquids.  In  this  respect,  the  method  of  the  second  temperature-time 
intervals  merits  attention  [1,  2],  The  theory  of  this  method  allows  us  to 
use  it  in  investigating  solids,  as  well  as  liquids.  An  Important  advantage 
of  this  method  is  the  fact  that  this  method  is  comprehensive,  that  is, 
all  thermophysical  characteristics  a,  X  and  c  can  be  determined  from  a  single 
experiment.  This  is  especially  important  in  studying  solutions  of  polymerB, 
since  both  such  solutions  undergo  physicochemical  changes  In  the  course  of 
heating,  and  therefore  it  is  necessary  to  determine  all  thermophysical  char¬ 
acteristics  during  a  single  experiment. 

Several  versions  of  the  method  of  second  temperature-time  Intervals 
have  been  theoretically  formulated.  For  all  these  versions  the  method  of 
measurement  is  the  same  [1].  Working  formulas  for  calculating  thermophysical 
characteristics  are  of  the  form: 


The  dimensionless  parameters  p  and  e  that  are  part  of  the  working  formulas 
are  calculated  from  worl.'ng  tables,  which  are  of  different  types  for 
different  versions  of  the  method  of  second  tempereture-time  intervals  [3]. 

A  difficulty  in  measuring  thermophysical  characteristics  of  liquids 
lies  in  elimination  of  convection.  This  requires  that  the  heat  flux  be 
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directed  from  the  top  downward.  In  our  experiments  we  used  a  coolant, 
placing  it  downward  from  the  layer  of  the  test  liquid  or  heater,  which  was 
placed  facing  upward.  A  stream  of  water  from  the  TS-15  thermostat  served  sb 
the  coolant,  and  as  a  heater  —  both  the  stream  of  water  and  an  electric  hot 
plate  with  an  electronic  temperature  control. 

The  thickness  of  the  test  liquid  layers  was  1.5-2  mm,  which  satisfied 
the  well  known  relationship  Gr*Pr  <  600. 

The  following  were  used  as  reference  materials  for  heat  receivers: 
polymethylmethacrylate,,  ebonite,  quartz  sand,  vulcanized  rubber  and  marble. 

For  the  one-dimensional  case  of  heat  propagation,  the  lateral  surface,  of 
the  cell  and  the  heat  receiver  were  heat-insulated. 

The  most  promising  method  for  determining  thermophysical  characteristics 
of  liquids  is  the  generalized  method  1.  Here  the  thermocouple  is  embedded  to 
a  certain  depth  in  the  heat-receiver  material,  which  makes  it  possible  to 
conduct  studies  in  very  thin  layers,  free  of  convection. 

Thermophysical  characteristics  of  solutions  of  sugar  and  water,  gasoline 
in  avtol,  and  solutions  with  low  stratifying  were  Investigated,  such 
as  the  following:  acetone-water,  ethanol-water,  pyridine-water,  aniline- 
propyl  alcohol.  The  index  of  refraction  for  these  compounds  is  measured  at 
the  same  time  as  a  function  of  solution  concentration. 

It  was  found  that  most  of  the  solutions  studied  reveal  a  deviation  from 
the  additivity  rule.  For  some  of  them,  well  defined  anomalies  on  the 
composition  versus  property  diagrams  were  revealed. 

Figure  1  presents  thermophysical  characteristics  of  the  acetone-water 
system.  Here  also  is  shown  variation  of  the  index  of  refraction  as  a 
function  of  concentration.  It  is  clear  that  the  nature  of  the  change  of  the 
quantity  1/n  is  very  close  to  the  change  of  the  quantity  a. 

Figure  2  presents  thermophyaical  characteristics  of  the  pyridine-water 
system.  We  see  a  well  defined  nonmonotonic  variation  of  the  thermal  diffus- 
J.vity  coefficient  as  a  function  of  solution  concentration. 


The  nonmonotonic  nature  of  thermophysical  characteristics  as  functions 
of  concentration  is  even  more  graphically  evident  in  the  aniline-propyl 
alcohol  system  (Figure  3). 

Symbols 

\  m  thermal  conductivity  of  test  compound;  a  ■  thermal  dlffuBivity  of 
test  compound;  c  -  specific  heat;  p  “  density;  h  ■  fair ? '-’ss  of  layer; 
r  »  experimentally  measured  delay  time  of  thermal  front;  n  ■  index  of 
refraction  of  solution;  b  ■  instrument  constant;  c  -  bulk  concentration. 
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Figure  3.  Coefficients  X  and  a  as  fui 
concentration  for  the  aniline-propyl  i 
II,  I  -  f2(o) 


-520- 


References 

1.  Vol’kenshteyn,  V.  S.  and  N.  N.  Medvedev,  TFZh ,  Vol.  2,  No.  10,  1959. 

2.  Vol'kenshteyn,  V.  S.,  Teplo-  i  Maeeopevenos,  [Heat  and  Mass  Transfer], 

Vol.  7,  Nauka  i  Tekhnika  Press,  Minsk,  1960,  [sic]. 

3.  Vol'kensnteyr ,  V.  S.  and  N.  N.  Medvedev,  Sb.  Isaledovaniya  po  Teplo- 

provodno8tij  [Collection:  Research  on  Thermal  Conductivity],  Nauka  i 
Tekhnika  Press,  Minsk,  1966. 


i 

I 


i 

i. 


THERMO PHYSICAL  PROPERTIES  OF  SEVERAL  NONFILLED  AND  FILLED  AMORPHOUS 
POLYMERS  IN  THE  GLASSY  STATE1 

V.  P.  Dushchenko,  V.  M.  Baranovskiy,  I.  A.  Uskov, 

V.  S.  Tytyuchenko,  Yu.  M.  Krasnobokiy,  N.  I.  Shut, 

V.  V.  Levandovskiy  and  Yu.  T.  Tarasenko 

Thermophysical  properties  of  polymers,  especially  filled  polymers,  have 
not  been  very  closely  investigated. 

We  used  the  nonsteady-state  method  of  dynamic  heating  for  the  investiga¬ 
tion  [1], 

Heat  transfer  in  polymers  is  determined  to  a  considerable  extent  by  the 
supermolecular  structure  and  nature  of  the  intra-  and  intermolecular  actions 
[2]. 

Heat  scattering  in  a  cross-linked  polymer  —  the  ED-5  epoxide 
resin  —  is  the  greatest,  a  lower  value  marks  the  branched  polystyrene  (PS), 
and  still  lower  is  linear  polyvinyl  chloride  (PVC).  Powders  of  copper 
(particles  of  irregular  shape)  and  iron  (particle  shape  approximately 
spherical)  were  used  as  fillers  for  these  polymers. 

The  coefficient  of  thermal  conductivity  X  of  the  amorphous  polymers 
investigated  in  the  glassy  state  region  increased  slightly  with  rise  in 
temperature  (Figure  1).  After  the  glassy  point  was  reached,  when  the  macro¬ 
molecules  acquire  flexibility,  X  decreases  sharply.  This  is  due  to  the 
gradual  thawing  out  of  its  segmentary  vibrations. 

These  results  are  in  agreement  with  the  Debye  phonon  theory,  if  we  bear 
in  mind  that  the  specific  heat  of  a  polymer  rises  sharply  with  temperature, 
while  its  density,  velocity  of  propagation  of  sound  and  the  mean  free  phonon 
path  decrease  slightly. 


*See  p.  526 
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When  metal  fillers  are  incorporated  into  a  polymer,  a  certain  rise  in 
thermal  conductivity  is  observed,  which  is  appreciable  even  at  relatively 
small  filler  concentrations.  Evidently,  this  is  largely  associated  with  the 
formation  in  the  polymer  of  domains  of  oriented  structure  exhibiting  height¬ 
ened  thermal  conductivity  under  the  action  of  the  filler,  and  also  with  the 
appreciable  thermal  conductivity  of  filler  particles. 

A  substantial  rise  in  the  thermal  conductivity  coefficient  was  observed 
for  PS  and  P VC  for  signif least  filler  contents  (40-50  percent).  Thi'j  can  be 
associated  with  a  rise  in  the  proportion  of  oriented  domains  around  the 
filler  particles  in  the  overall  polymer  bulk. 

The  anomaly  noted  in  variation  of  thermal  conductivity  coefficients  with 
temperature  becomes  evident  especially  forcefully  in  the  transitional  region 
from  the  glassy  state  to  the  highly  elastic  state.  Relaxation  effects  in 
amorphous  polymers  are  also  most  clearly  manifest  in  this  temperature  range 
[3].  It  is  obvious  that  this  anomaly  is  inherent  only  to  that  section  of  the 
polymer  which  undergoes  the  strong  variation  of  supermolecular  structure  as  a 
result  of  the  effect  of  the  solid  surfaces  of  filler  particles.  It  can  be 
assumed  that  in  compositions  with  high  filler  content  the  main  bulk  of  a 
polymer  proves  to  be  modified,  which  leads  to  a  rise  in  the  thermal  conduct¬ 
ivity  coefficient  with  temperature.  As  flexibility  of  macromolecules  and 
their  clustered  formations  grows  with  rise  in  temperature,  the  role  of 
structural  changes  in  the  polymer  decreases  and  the  filled  composition 
behaves  in  the  ordinary  way  —  the  coefficient  of  thermal  conductivity  passes 
through  a  maximum  and  begins  to  decrease. 

A  drop  in  the  thermal  diffusivity  coefficient  of  polymers  (a  •  X/cp) 
with  rise  in  temperature  (Figure  2)  is  accounted  for  by  the  temperature 
dependence  of  the  thermal  conductivity  coefficient,  the  specific  heat  c  and 
the  density  p. 

The  coefficient  of  thermal  diffusivity  of  the  filled  polymers  studied 
decreases  with  rise  in  temperature  (Figure  2) ,  which  is  directly  associated 
with  variation  of  the  thermal  dif fupivi ty  coefficients  of  the  components. 

The  coefficient  of  thermal  dlffue'  ity  rises  for  strongly  filled  compositions 
when  t  *  const  in  the  form  of  the  appearance  of  direct  contact  ~  "thermal 
bridges"  between  the  filler  particles. 
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Figure  1.  Coefficient  of  thermal  conductivity  of  polymers  as  a 
function  of  temperature:  a,  PVC  +  Fe;  b,  PVC  +  Cu;  c,  PS  +  Fe; 
d,  ED-5  +  Cu:  x  -  0  percent,  •  -  10  percent,  0-20  percent, 

A  -  30  percent,  A  -  40  percent.  ©  -  50  percent  (Fe  or  Cu) 


Figure  2.  Coefficient  of  thermal  diffusivity  as  a  function  of 
temperature:  a,  PVC  +  Fe;  b,  PVC  +  Cu}  c,  PS  +  Fe;  d,  ED-5  +  Cu: 
x  -  0  percent,  •  -  10  percent,  0-20  percent,  A  -  30  percent, 

A  -  40  percent,  0-50  percent  (Fe  or  Cu) 


Footnotes 
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MASS  TRANSFER  CHARACTERISTICS  OF  PLASTER  WITH  FILLER 
N.  P.  Zlobina 

The  intensity  of  moisture  transfer  in  construction  materials  character¬ 
ized  by  the  coefficient  of  the  potential  conductivity  of  mass  transfer  and 
the  thermal  gradient  coefficient,  which  depend  on  the  moisture  content  and 
temperature.  The  aim  of  this  study  has  been  to  determine  the  coefficient  of 
potential  conductivity  and  the  thermal  gradient  coefficient  of  plaster 
and  its  components  as  functions  of  moisture  content  at  a  temperature  of 
298®K. 

Study  of  the  coefficient  of  potential  conductivity  of  plaster,  gypsum 
rock  and  slag  sand  was  carried  out  by  a  momentum  method.  Essentially  the 
method  amounts  to  the  following:  the  end  of  a  rod  made  of  the  test  material 
with  moisture-insulated  side  surfaces  is  wetted  with  a  flat  pulsing  moisture 
source . 

After  a  certain  time  interval  of  moisture  pulse  action  the  time  of 
onset  of  the  moisture  content  maximum  in  the  material  layer  at  a  distance  x 
from  the  uoist  end  of  the  rod,  the  coefficient  of  the  potential  conductivity 
is  determined  from  the  following  relationship  [1,  2,  7]  from  analogy  with  the 
momentum  method  of  determining  the  thermal  diffueivity: 
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The  quantities  and  x^^  lire  determined  from  experiment*. 

Practically  speaking,  measurement  of  the  coefficient  of  potential 
conductivity  by  the  momentum  method  is  carried  out  as  follows.  In  a  cylinder 
made  of  organic  glass  6-10  cm  high,  3.4  cm  at  a  distance  x  from  the  moist 

end,  two  parallel  nichrome  wires  1  mm  in  diameter  serving  as  electrodes  are 
secured.  The  value  of  the  distance  x  and  the  cylinder  height  depend  on  the 
properties  of  the  material  and  are  selected  in  such  a  way  that,  first,  the 
condition  of  unboundedness  which  is  postulated  in  solving  the  differential 
equation  of  mass  conductivity  is  not  violated,  and,  second,  the  distance  has 
to  be  selected  in  such  a  way  that  the  time  of  onset  of  moisture  content 
maximum  in  the  selected  cross  section  and  the  time  of  moisture  pulse 

action  do  not  coincide,  since  in  this  case  the  function  $  is  not  determln- 

a 

ate.  Copper-constantan  thermocouples  are  placed  in  thia  same  cross  section 
x. 

The  cylinder  prepared  in  this  way  is  filled  with  the  test  materiel, 
which  is  moistened  to  the  desired  level  of  moisture  content.  Then  the  ends 
of  the  cylinder  are  sealed  hermetically  with  metal  sheets  fitted  with  rubber 
inserts. 

This  cylinder  is  placed  in  the  thermostat  at  T  -  298#K  until  uniform 
moisture  distribution  over  the  entire  specimen  is  achieved. 

In  carrying  out  the  experiment,  the  moist  end  of  the  cylinder  is  lowered 
into  water,  and  the  time  of  action  of  the  moiBture  pulse  x^  Is  established. 
Then  the  cylinder  is  hermetically  sealed  and  placed  in  the  thermostat. 
Electrodes  are  connected  to  a  self-recording  PSR-1-01  automatic  potentio¬ 
meter.  We  know  that  the  electrical  resistance  of  the  material  depends  on  Its 
moisture  content  [4].  Resistance  as  a  function  of  moisture  content  for 
almost  all  capillary-porous  bodies  can  be  expressed  by  a  power  function. 


where  A  and  k  are  positive  constants  dependent  on  the  test  material  and  the 


-529- 


conditions  of  measurement .  We  see  from  examination  of  the  function  (3)  that 
the  resistance  minimum  corresponds  to  a  moisture  content  maximum. 

Therefore,  the  time  of  onset  of  moisture  content  maximum  is  determined 
from  the  onset  of  moisture  content  maximum  Is  determined  from  the  onset  of 
the  electrical  resistance  minimum.  V,\s  muBt  note  that  only  the  nature  of  the 
variation  of  moisture  content  in  the  selected  cross  section  as  a  function  of 
time  is  established  by  this  method,  and  that  the  moisture  content  value  Is 
not  determined. 

The  nature  of  variation  of  electrical  resistance  in  some  selected  crosB 
section  as  a  function  of  time  for  a  given  time  interval  of  moisture  pulse 
action  ■  7  sec  for  plaster  Is  shown  in  Figure  1.  In  the  figure  the  time 
of  onset  of  the  electrical  resistance  minimum  Is  easily  seen,  which  corre¬ 
sponds  to  x  ,  the  time  of  onset  of  the  moisture  content  maximum  in  the 
max 

selected  cross  section. 

ThiB  method  was  used  to  determine  the  coefficient  of  potential  conduct- 

o 

lvity  of  plaster  (y  »  1300-1400  kg/m  )  of  gypsum  rock  (water-cement  ratio 
°  3 

0.5),  slag  sand  (yq  ■  1600  kg/m  )  as  a  function  of  moisture  content  at  a 
temperature  of  298°K  (cf .  Figures  2  and  3) . 

The  nature  of  variation  of  the  coefficient  of  potential  conductivity 
for  plaster,  gypsum  rock  and  slag  sand  correspond  to  studies  made  by 
other  authors  obtained  for  different  capillary-porous  bodies  [3,  4,  6]  and 
can  be  accounted  for  by  the  nature  of  distribution  over  radii  [3]. 

The  thermal  gradient  coefficient  of  these  materials  was  determined  by  a 
steady-state  method.  A  cylinder  3.4  cm  in  diameter  and  12  cm  high  made  of 
the  test  material  with  the  desired  level  of  moisture  content,  with  moisture- 
insulated  total  surface  and  heat-insulated  lateral  surface,  was  placed 
between  two  heater  elements.  In  our  case,  hot  plates  from  the  desks  of 
Doctor  Bokk,  the  Weiss  Company  for  determining  the  coefficient  of  thermal 
conductivity  were  used.  The  temperature  of  the  hot  plate  was  kept  constant 
at  80°C,  the  temperature  of  the  cold  plate  was  kept  at  15°C.  For  the 
temperature  drop  induced  in  the  specimen  of  test  material,  redistribution  of 
moisture  content  occurred  only  under  the  action  of  temperature.  After 
arrival  at  the  steady-state  the  temperature  field  was  recorded. 


Dividing  up  the  cylinder  by  layers,  the  moisture  content  of  each  layer 
wan  determined  by  a  weight  method.  Based  on  the  data  obtained,  the  following 
graphs  were  plotted:  T  -  f(x)  and  u  ■  f(x),  where  x  is  the  distance  of 
reference  points  from  the  hot  end  of  the  cylinder. 

The  temperature  gradients  and  the  moisture  content  gradients  were 
determined  from  the  graphs. 

The  thermo-gradient  coefficient  was  calculated  from  the  formula 
given  in  [8] 

(5  _  v  u  fcg/feg  ; 

V  t 

From  the  results  of  investigating  the  thermal  gradient  coefficient  of 
plaster  and  gypsum  rock  (Figure  4) ,  it  follows  that  in  the  region  of  sorption 
moisture  content,  the  thermal  gradient  coefficient  increases  and  reaches  a 
maximum  for  maximum  sorption  moisture  content,  but  in  the  super-sorption 
region  the  thermal  gradient  coefficient  drops  off  with  rise  in  molBture 
content.  This  kind  of  variation  in  thermal  gradient  coefficient  also 
corresponds  to  atudlea  made  by  other  authors  [3,  4]. 

Symbols 

aifl  -  coefficient  of  potential  conductivity;  ■  time  of  action  of 

moisture  pulse;  x  *  distance  from  moist  end;  t  ■  time  of  onset  of  moisture 
r  max 

content  maximum  in  selected  cross  section  x;  $a  ■  auxiliary  function; 

R  *»  electrical  resistance;  u  -  moisture  content;  yq  *  specific  weight  of 
material;  6  -  thermal  gradient  coefficient. 


-531' 


Figure  1.  Variation  of  electrical  resistance  as  a  function  of 
time  for  a  given  time  interval  of  moisture  pulse  action  T-  *  7  sec 
for  plaster 


Figure  2.  Coefficient  of  potential  conductivity  as  a  function  of 
moisture  content  for  T  ■  298#K  for:  1,  Plaster,  y  •  1300- 


Figure  3,  a^  as  a  function  of  moisture  content  for  gypsum  rock 
at  T  -  298°K 
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Figure  4.  Thermal  gradient  coefficient  as  a  function' of  moisture 
content  at  T  «  313#K  for:  1,  Plaster;  2,  Gypsum  rock 
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TRUE  DIFFUSION  COEFFICIENTS  OF  THE  He-C02  SYSTEM  AS  A  FUNCTION 

OF  TEMPERATURE 

N.  D.  Kosov  and  A.  F.  Bogatyrov 

• 

Heat  transfer  by  thermal  agitation  with  strict  observance  of  isobaricity 

•  of  conditions  is  charactererized  by  true  diffusion  coefficients  [1,  2].  The 
true  diffusion  coefficient  of  one  gas  in  another  is  not  equal,  in  the  general 
case,  to  the  true  diffusion  coefficient  of  the  latter  gas  in  the  former,  and 
according  to  Boltzmann  [1]  for  solid  spheres  with  account  taken  of  the 
persistency  of  velocities  [3,  4]  and  the  mean  free  path  as  a  function  of 
velocity  [3]  this  coefficient  is  determined  by  the  following  formula: 

D  ,  ■  1.051v/5kT 

1  3n jr^rrri^O  (i) 

The  true  diffusion  coefficient  of  the  second  gas  in  the  first  is  obtained 
from  formula  (1)  by  replacing  the  index  1  by  the  index  2. 

True  diffusion  coefficients  as  functions  of  temperature  have  as  yet  been 
poorly  investigated. 

♦  Beit..-  ’e  set  forth  the  results  of  measurements  of  true  and  mutual 
diffusion  coefficients  of  the  He-CC^  system  in  the  temperature  range  290- 

‘  430°K,  and  also  the  values  of  the  baro-effect  induced  when  there  is  diffusion 

within  closed  systems. 

Method  of  Measurements.  Results 

True  diffusion  coefficients  were  measured  by  a  steady-state  method  under 
isobaric  and  isothermal  conditions  [2].  The  difference  between  our  apparatus 
and  that  described  in  [2]  lies  mainly  in  the  thermostatlng  of  its  main 
assemblies  (diffusion  cell,  Interference  cells,  and  rheometers)  and  in  the 
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True  and  Mutual  Diffusion  Coefficients  of  the  He-CO  System  at  Different  Temperatures 
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3  percent.  The  diffusion  coefficients  of  carbon  dioxide  gas  in  helium, 
measured  and  calculated,  differ  from  each  other  by  35  percent  (precision  of 
Dj  5%  measurement).  However,  the  theoretical  curves  of  as  a  function  of 
temperature  followed  the  experimental,  but  was  shifted  downward  almost 
parallel  to  itself.  This  can  be  accounted  for,  firBt  of  all,  by  the  fact 
that  the  carbon  dioxide  molecule  is  far  from  a  solid  sphere.  Secondly,  as  we 
see  from  formulas  (2)-(4),  the  error  in  determining  the  effective  diameter 
had  a  considerably  greater  effect  on  the  precision  of  calculation  of  D^. 

Thus,  the  Boltzmann  theory  correctly  reflects  to  the  first  approximation 
the  course  of  the  temperature  dependence  of  true  diffusion  coefficients. 
However,  the  deviation  observed  (in  particular  for  the  diffusion  of  a  heavy 
gas  in  a  light)  requires  additional  investigation. 

The  Bato-ef feet  at  Different  Temperatures 

One  of  the  interesting  problems  directly  associated  with  investigation 
of  the  temperature  dependence  of  diffusion  coefficients  la  the  diffusion 
baro-effect  as  a  function  of  temperature.  According  to  the  Boltzmann 
diffusion  theory,  the  baro-effect  is  a  consequence  of  inequality  intrinsic  to 
diffusion  flows  of  the  first  and  second  gases.  In  closed  systems,  this 
leads  to  the  induction  of  pressure  on  the  part  of  ithe  gas  that  has  the  lower 
true  diffusion  coefficient,  which  in  turn  is  the  cause  of  the  hydrodynamic 
flow  of  the  gas  mixtures  toward  the  gas  that  has  the  higher  true  diffusion 
coefficient.  In  the  steady-state,  the  specific  hydrodynamic  flow  equals 
the  difference  between  the  diffusion  flows  and  for  the  single-dimensional 
case  can  be  written  as: 

»> 

For  given  gases  and  given  geometry,  the  hydrodynamic  flow  will  be  pro¬ 
portional  to  the  steady-state  pressure  difference.  Bearing  in  mind  equation 
(5)  and  taking  Ac  -  1,  we  will  obtain 


£k  p  D,  -  Dj_ 


(6) 


narrow  temperature  range  is  well  described  by  the  power  formula  in  [6],  From 
our  data,  the  exponent  for  the  mutual  diffusion  coefficient  was  found  to  be 
1.9,  which  within  limits  of  experimental  error  agrees  with  the  value  of  1.92 
found  in  [7].  This  betokens  the  satisfactory  performance  of  the  equipment  bb 
a  whole.  The  existing  small  difference  (within  the  limits  of  3  percent) 
between  our  data  and  the  data  in  [7]  is  accounted  for  by  the  fact  that  we  do 
not  Introduce  corrections  for  the  end  effect,  and  the  instrument  constant  was 
calculated  directly  from  geometrical  parameters  while  at  the  same  time  in  [7] 
the  geometric  parameter  was  determined  by  standardizing. 

The  true  diffusion  coefficients  of  the  function  of  temperature  are  more 
complex,  as  we  can  see  from  the  curves  in  Figure  1. 


The  diffusion  coefficient  measured  in  the  experiment  is  integral,  and 
that  given  by  formula  (1)  is  local.  The  true  integral  coefficient  of  diffu¬ 
sion  of  the  first  gas  in  the  second  obviously  equals t 


where 


C, 


1  c;  (6-c»Xr0«)  (i-Qn)c**a,i 

■n  20<n*  aws 

3ar0rrn, (1- to* ) p$ *  * yfn, p \ J1  )  > 


(2) 

C3) 


.  _  /  mt  r  mz_ 

11  2 /  V 


(4) 


The  true  integral  diffusion  coefficient  of  the  second  gas  in  the  first 
is  obtained  from  formulas  (2)- (4)  by  replacing  the  index  1  by  the  index  2. 

Coefficients  and  D2  calculated  from  formula  (2)  are  shown  in  Figure  1 
by  dashed  lines.  In  the  calculations  of  and  Dg,  the  values  of  the 

effective  diameters  were  calculated  from  the  self-diffusion  coefficients  at 
the  given  temperature  [6]  using  formula  (3).  As  we  can  see  from  the  curves 
in  Figure  1,  the  true  diffusion  coefficient  of  helium  in  carbon  dioxide, 
calculated  from  formula  (2)  differs  somewhat  (within  the  limits  of  5  percent) 
from  the  measured  values.  The  measurement  precision  of  the  coefficient  is 
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variation  of  the  geometric  parameters  of  the  diffusion  cell,  which  simul¬ 
taneously  also  fulfill  the  function  of  heat  exchanger.  The  thermostating 
liquid  for  the  diffusion  cell  was  dibutylphthalate,  and  for  the  rheometers 
and  interference  cells  the  temperature  of  which  was  always  kept  constant 
(298°K)  —  water.  Compressive  gaskets  were  made  of  Teflon.  Preheated  to 
temperatures  close  to  the  experimental  temperature,  the  gases  were  swept 
through  at  a  constant  rate  (the  hulk  rate  throughout  all  the  experiments  was 
kept  constant  at  172  *  10  m  /sec)  through  the  tubings  of  the  diffusion 
cell  in  which  they  took  on  the  desired  temperature.  The  tubings  were 
connected  with  a  set  of  capillaries  in  which  diffusion  was  carried  out.  The 

capillary  sets  were  made  of  approximately  identical  tubes  of  stainless  steel 

-5  2 

(injection  needles)  with  a  total  area  of  3.77  x  10  m  (180  tubes  in  one  and 

_2 

26  tubes  in  the  other).  The  tube  length  was  7.05  x  10  m.  The  tubings  were 
Blitlike  in  form  (height,  width  and  length,  respectively,  were  0.002  x  0.08  x 
x  0.5  m)  and  made  of  brass.  Special  grooves  were  ground  over  the  entire 
outer  surface,  covered  over,  and  served  to  channel  the  dibutylphthalate  from 
the  thermostat.  To  reduce  heat  losses  and  equalize  temperature,  the  diffu¬ 
sion  cell  with  the  capillary  set  was  carefully  heat-insulated.  Gases  after 
diffusion  were  cooled  to  the  temperature  of  the  analysis  in  a  special  heat- 
exchanger  and  then  sent  to  the  interference  cells  of  the  Interferometer. 
Establishment  of  the  isobaricity  of  the  conditions  and  measurement  of  the 
concentrations  of  the  diffusing  compounds  were  carried  out  as  in  [2]. 

The  coefficients  of  mutual  diffusion  were  measured  on  this  same 
equipment  using  the  steady-state  method  described  in  [5]  for  the  same 
temperature  range.  For  this  purpose  the  pressure  was  varied  in  one  of  the 
tubings  until  the  hydrodynamic  flow  induced  as  a  consequence  equaled  the 
overall  transport  of  a  particular  gas. 

The  results  of  measuring  true  diffusion  coefficients  of  helium  in 
carbon  dioxide  gas  and  of  carbon  dioxide  gas  in  helium  D^,  and  also  the 
mutual  diffusion  coefficients  are  presented  in  Table  1. 

Figure  1  presents  in  logarithmic  scale  the  diffusion  coefficients  D., , 
and  of  the  system  Investigated  as  a  function  of  temperature.  The 
temperature  dependence  of  the  mutual  diffusion  coefficients  of  gases  for  a 


-539- 


I 


Another  point  of  view  on  the  nature  of  the  diffusion  baro-effect  [8,  9] 
for  the  sane  conditions  loads  to  a  function  of  the  form 

&p  ~jd14  <7> 

The  true  diffusion  coefficients  and  the  mutual  diffusion  coefficient 
varied  dissimilarly  with  temperature.  Therefore,  fozmulas  (6)  and  (7)  give 
the  baro-effect  es  different  functions  of  temperature.  Therefore  it  was  of 
definite  interest  to  Investigate  the  baro-effect  at  different  temperatures. 

The  same  apparatus  used  in  measuring  mutual  diffusion  coefficients  was 

employed  in  measuring  the  value  of  the  baro-effect  at  different  temperatures 

-4 

(mean  capillary  radius  was  equal  to  (5.8  *  0.5  *  10  m).  In  the  stdady- 

state  method  the  baro-effect  value  Ap  corresponds  to  that  pressure  value  at 
which  the  overall  flows  of  a  given  gas  are  equalized.  The  pressure 
difference  was  measured,  as  in  [2],  from  the  change  In  the  height  of  the 
level  of  gas  exiting  from  the  tubing.  Figure  2  presents  the  results  of 
measurement  of  the  baro-effect  at  different  temperatures  for  the  He-COg  and 
He-Ar  systems,  and  the  diffusion  coefficients  of  the  latter  as  a  function  of 
temperature  have  been  investigated  earlier  [10].  Aa  we  can  see  from  Figure  2, 
the  baro-effect  varied  linearly  with  temperature.  Figure  3  presents  the 
quantity  as  a  function  of  temperature  plotted  from  the  data  in 

Table  1  and  in  [10],  this  function  also  proving  to  be  linear.  The  mutual 
diffusion  coefficient  is  proportional  to  Tr,  where  r  >  1.5  (for  the  He-COj 
system  r  -  1.92,  for  the  He-Ar  system,  r  »  1.75  [6]).  Therefore,  according 
to  formula  (7),  Ap  must  be  proportional  to  Tc,  which  was  not  observed  in  the 
experiment.  ' 

Symbols 

and  Dj  ■  true  local  diffusion  coefficients;  and  Dj  “  true 
integral  diffusion  coefficients;  and  •  coefficients  of  mutual 
diffusion  and  of  self-diffusion;  u  -  persistency  of  velocities;  m  and 
a  ■  mass  and  effective  diameter  of  the  molecule;  n  ■  number  of  molecules  per 
unit  volume *  relative  concentrations  c£-C£  -  concentration 
difference  at  ends  of  capillary  of  length  L;  k  •  Boltzmann  constant;  p  and 
T  -  pressure  and  absolute  temperature. 
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Figure  1.  Diffusion  coefficients  of  the  He-COj  system  as  a  function 

of  temperature:  1,  For  the  true  diffusion  coefficient  of  helium 
in  carbon  dioxide;  2,  For  the  true  diffusion  coefficient  of 
carbon  dioxide  in  helium;  3,  For  the  mutual  diffusion  coefficient. 
The  solid  line  is  experimental  values;  the  dashed  line  is  theor¬ 
etical  values. 


Figure  2.  Diffusion  baro-effect  aa  a  function  of  temperature  for 
the  He-Ar  system  (1)  and  for  the  Ho.-COj  system  (2) 


Figure  3.  Difference  of  true  diffusion  coefficients  as  a  function 
of  temperature  for  the  following  systems:  1,  For  He-Ar;  2,  For 
He-C02 
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NONISOTHERMAI  DIFFUSION  IN  GASES 
N.  D.  Kosov,  A.  F.  Bogs ty rev  and  L.  I.  Kurlapov 

The  diffusion  of  two  gases  at  constant  pressure  under  Isothermal  condi¬ 
tions  is  caused  both  by  a  difference  in  the  density  of  the  number  of  mole¬ 
cules  as  well  as  by  a  difference  in  their  thermal  velocities.  Both  of 
these  factors  act  simultaneously,  therefore  they  must  also  be  taken  into 
account  simultaneously  when  examining  nonisothermai  diffusion.  Usually 
nonlsothermal  diffusion  of  a  binary  mixture  of  gases  at  constant  pressure  is 
called  thermodiffusion  [1].  If  under  nonisothermai  conditions  molecular 
transport  caused  by  a  difference  in  their  population  densities  and  a 
difference  in  thermal  velocity  is  taken  into  account,  then  even  in  clean  gas 
a  dlffuslonal  molecular  flux  from  the  cold  region  to  the  hot  becomes  manifest 
[2].  This  molecular  transport  is  conveniently  called  thermo-self-diffuslon 
(from  analogy  with  self-diffusion  under  Isothermal  conditions).  The 
thermo-self-diffuslon  flux  leads  to  a  rise  in  the  molecular  population 
density  in  the  hot  region  (compared  with  the  equilibrium  region  at  a  given 
temperature)  and  to  a  decrease  in  the  cold  region,  that  is,  to  the  appearance 
of  a  pressure  difference  between  these  two  gas  domains.  The  baro-effect 
Induced  when  there  is  nonisothermai  diffusion  1b  conveniently  called  the 
thermodiffusive  baro-effect  [3]. 

Thermo-self-diffusion  graphically  accounts  for  the  mechanism  of 
thermodiffusive  separation  of  gases  as  the  diffusion  of  each  component  from 
the  cold  region  to  the  hot  at  its  own  velocity. 

The  present  study  examines  a  steady-state  method  of  measuring  the 
coefficient  o,f  thermo-self-diffuslon  and  the  results  of  experiments  analyzing 
the  physical  sense  of  thermodiffusive  separation. 
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Steady-State  Method  of  Measuring  the  Thnrmo-Self-Dif fusion  Coefficient 

The  thermo-self-diffusion  flow  (and  interms  of  it  the  coefficient  of 
thermo-self-diffusion)  can  be  measured  on  an  apparatus  designed  to  measure 
the  thermodlffuslve  baro-effect  in  gases  [2]  and  shown  schematically  in 
Figure  1,  modifying  somewhat  the  method  of  observation. 

Let  the  left  vessel  be  kept  at  a  constant  temperature  T^  (hot  region  of 
gas),  and  the  right  at  a  temperature  T ^  (cold  region  of  the  gas).  The 
vessels  are  connected  by  a  wide  capillary  (2r  ►  X),  such  that  the  wall 
effects  can  be  neglected.  Owing  to  thermo-self-diffusion  in  section  A  of  the 
capillary  (Figure  1),  an  increased  pressure  is  Induced,  and  In  section  B  —  a 
reduced  pressure.  These  pressure  changes  give  rise  to  hydrodynamic  gas  flows 
in  the  capillary  from  the  cross  section  A  to  both  sides  and  to  the  section  B, 
and  also  lead  to  the  movement  of  a  drop  of  lnvlscous  silicon  oil  in  the 
horizontal  tube  from  the  hot  vessel  to  the  cold  at  a  specific  velocity  V. 

In  the  steady-state  the  number  of  gas  molecules  intralned  by  the  hydro- 
dynamic  flow  to  both  sides  from  section  A  equals  the  number  of  gas  molecules 
arriving  at  the  section  of  the  cold  region  owing  to  thermo-self-diffusion, 
that  is, 


The  number  of  molecules  transported  per  unit  time  by  the  hydrodynamic 
flow  through  the  capillary  from  region  A  to  region  B  (regarding  the  movement 
in  the  capillary  as  of  the  Polseuille  type),  we  will  have 


Apr1  S’ 

t8* 


(2) 


An  analogous  number  of  molecules,  but  transported  through  the  manometric 
tube,  can  be  written  as  follows: 


6P*. 


(3) 
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length  In  the  steady  state  vary  linearly,  as  shown  in  the  upper  part  of 
Figure  1,  then 


Ap  -  Ap*  ♦  A p,  4  Apa  . 


(4) 


The  pressure  difference  Apk  under  the  effect  of  which  the  drop  is  moved 
can  be  found  from  its  rat^  of  movement  through  the  force  of  friction 
f  -  oV 


s. 


„.oCV 

“X 


(5) 


Determining  Ap.^  and  Ap2  from  equation  (3)  and  substituting  these  values 
as  well  as  Apfc  in  equation  (4)  ve  find 

jtf[t*T^(h<‘T'*h‘T4v-  «) 

Substituting  (6)  in  (2),  and  then  (2)  and  (3)  in  (1),  we  will  find  the 
final  expression  for  the  thermo-self-diffusion  flow 


The  thermo-self-diffusion  flow  can  be  calculated.  By  using  the  ordinary 
method  of  the  kinetic  theory  of  gas  [3],  Laranjeira  [4]  obtained  the 
following  expression  for  the  overall  flow  T of  molecules  of  the  species  i 
per  unit  area  per  unit  time: 

H  niWi[AiVdnct4'AiV^np-(At-|S)vtnT].  (8) 

Neglecting  baro-d if fusion,  we  obtain  from  equation  (8)  the  following 
expression  for  the  thermo-self-diffusion  flow: 


(9) 


Jt 


T)  sdlnr 

"tS  dx 


Equating  the  rlghthand  members  of  equations  (7)  and  (9),  we  find  the 
value  of  the  thermo-self-diffuBion  coefficient. 

To  measure  thermo-self-dif fusion  coefficients,  two  units  of  equipment 

were  built  with  different  geometric  parameters  (for  one  unit:  r  •  3.12  x 

/  -2  —2  —6  2 
x  10  m{  L  ■  9.5  x  io  mj  l.  »  *  2.2  x  10  ms  S_  ■  4.15  x  10  m  ;  for 

_4  2  _o  0  -2 

the  second:  r  ■  3.42  x  10  mj  L  *  8.8  *  10  n;  Z  *  J.  ■  2.5  x  10  m} 

—6  2  X  z 

S  ■  2.83  x  io  m  ).  The  measurements  were  conducted  at  T,  ■  370°K, 

0  -2^-22 
T2  -  285°K,  and  at  different  pressures  (from  90  x  10  to  480  x  10  N/m  ). 

The  measured  values  of  the  thermo-self-dif fusion  coefficient  prove  to  be 

(0.9  *  0.1)  x  io  4  and  (1.6  *  0.2)  x  io-4  m2/sec,  respectively,  for  hydrogen 

and  helium. 

According  to  Jeans  [3]  X'  ■  0.745  X.  Then  the  self-diffusion 
coefficients  of  hydrogen  and  helium  calculated  according  to  formula  (9)  in  . 
terms  of  the  thermo-self-diffuslon  coefficients  prove  to  be,  respectively, 

1.4  x  io~4  and  2.5  x  10  4  m2/sec  for  a  mean  temperature  of  328°K,  which  is  in 
agreement  with  literature  data  [5], 

Physical  Meaning  of  Thermodiffusive  Separatirn  of  Gases 

As  we  know,  no  well  defined  explanation  of  thermodiffusion  has  yet  been 
given  [6].  Thermo-self-dif fusion  allows  us  to  account  for  the  phenomenon 
of  thermodiffusive  separation  as  diffusion  under  nonisothemal  conditions  of 
each  component  from  the  cold  region  of  the  gas  to  the  hot  at  its  own 
velocity.  This  leads,  as  in  the  case  of  thermo-oelf-dif fusion,  to  the 
Induction  of  the  thermodiffusive  baro-effect  in  a  gas  mixture,  whose  value 
depends  on  the  hydrodynamic  resistance  between  the  hot  and  cold  regions  of 
the  gas.  In  a  closed  system  under  steady-state  conditions  transport  by 
thermodiffusion  of  each  component  from  the  cold  region  is  equilibrated  by 
reverse  transport  by  the  hydrodynamic  flow  from  the  hot  region  of  the  gas. 

The  thermodiffusive  baro-er'fect,  as  a  phenomenon  that  is  diffusional  in 
nature,  must  be  inversely  proportional  to  pressure.  This  does  occur  for  a 
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clean  gas  [2].  Measurements  of  the  thermodif fusive  baroreffect  were  made  on 
the  two  equipment  units  described  above  for  different  pressures  for  a  mixture 
of  68  percent  helium  with  carbon  dioxide  gas.  The  results  of  the  measure¬ 
ments  are  shown  in  Figure  2.  As  we  can  see  from  Figure  2,  the  product  Ap',p 
remains  constant,  which  emphasizes  the  dlffusional  nature  of  the  thermo- 
diffusive  baro-effect  also  for  a  gas  mixture.  The  mean  deviation  from  the 
constant  is  13  percent  for  the  upper  iine  (first  unit)  and  5  percent  for  the 
lower  (second  unit). 

The  hydrodynamic  flow  Induced  as  a  consequence  of  the  thermodif fuBive 
baro-effect  worsens  gas  separation.  The  shift  In  concentrations  will  be 
greater  if  by  some  method  the  hydrodynamic  flow  could  be  eliminated.  In  the 
case  of  isothermal  diffusion  the  hydrodynamic  flow  can  be  eliminated  by 
maintaining  the  pressure  at  capillary  ends  identical  and  continually  sweeping 
gas  through  them  [7].  We  used  an  analogous  diffusion  cell  in  studying  the 
concentration  shift  as  a  function  of  intensity  of  hydrodynamic  gas  flow 
through  the  capillary. 

Let  the  test  mixture  sweep  through  the  tubings  of  the  thermodif fusive 
cell  (cf.  Figure  3)  at  a  specific  bulk  flow  rate  q  (in  the  upper  tubing  at  a 

d 

temperature  of  T^  and  in  the  lower  —  at  a  temperature  T2;  we  set  T^  >  Tj). 

The  tubings  are  connected  by  capillary  (to  Increase  measurement  precision  a 

set  of  51  capillaries  is  used)  in  which  thermodiffusion  and  macroscopic 

movement  of  the  gas  mixture  occurs.  The  set  of  capillaries  was  made  of  brass 

tubes  with  a  mean  internal  diameter  of  3.0  *  10~^  m  and  a  length  of  0.15  m, 

therefore  a  slight  pressure  difference  is  Induced  at  its  ends,  according  to 

2 

calculations  of  the  order  of  several  hundredths  of  a  dyne/cm  .  It  was  not 
possible  to  employ  replacement  of  the  capillary  set  to  establirh  the  condi¬ 
tion  Ap  ■  0,  as  was  done  in  [7].  Therefore  the  value  of  the  hydrodynamic 
flow  £f  gas  mixture  through  the  capillary  was  varied  by  sucking  in  the  gas  at 
the  outlet  from  the  thermodif fusive  cell  initially  from  the  lower  tubing,  and 
then  --  from  the  upper.  Further,  the  bulk  flow  rate  (using  the  displacement 
method)  at  the  outlet  from  the  given  tubing  and  the  concentration  of 
components  in  the  mixture  at  the  outlet  from  both  tubings  were  measured. 
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During  all  the  experiments  and  Tn  remained  constant  (411°K  and  286aK, 

respectively),  and  the  bulk  flow  rates  at  the  inlet  to  the  thermodif fuaivc 

_6  3 

cell  in  both  tubings  were  identical  and  equal  to  0.908  *  10  m  /sec. 

Figure  4  shows  the  variation  in  the  helium  concentration  in  the  upper 
tubing  as  a  function  of  change  in  bulk  flow  rate  <qQ  -  q)  (curve  1)  and  the 
variation  in  carbon  dioxide  concentration  in  the  lower  tubing  (curve  2) .  We 
note  that  the  bulk  flow  rate  at  the  inlet  to  the  thermodif fueive  cell  and  at 
the  outlet  from  the  tubings  were  always  measured  at  the  same  temperature. 
Curves  1  and  2  intersect  at  a  point  for  which  qQ  *  q  (within  limits  of 
experimental  error).  Under  these  conditions  thermodlffusion  occured  in  such 
a  way  that  the  bulk  flow  rates  at  the  outlet  from  the  thermodif fuaive  cell 
remained  the  same  as  at  the  inlet.  This  was  possible  in  two  cases. 

First  when  the  mixture  components  under  the  effect  of  a  temperature 
gradient  diffuse  to  opposite  sides  (the  light  component  to  the  hot  region, 
and  the  heavy  component  to  the  cold)  in  equal  amounts  (by  numbers  of 
molecules  transported)  in  the  absence  of  a  hydrodynamic  flow.  In  this  case 
thermodiffusive  flows  will  be  at  a  maximum  and  the  concentration  shifts  will 
be  the  highest.  If,  however,  thermodiffusive  flows  are  not  identical  in 
magnitude  but  are  directed  to  opposite  sides,  equilibrating  them  is  possible 
only  by  means  of  the  hydrodynamic  flow,  as  in  the  case  of  isothermal  mutual 
diffusion.  Flows  equalized  in  this  way  will  be  characteristics  of  the 
overall  transport  —  by  thermodif fusion  and  by  the  hydrodynamic  flow.  In 
this  case  the  concentration  shift  will  not  attain  a  maximum,  since  the 
hydrodynamic  flow  will  mix  the  mixture  being  divided.  If  we  eliminate  the 
hydrodynamic  flow  In  some  manner  (that  is,  conduct  an  experiment  under 
strict  compliance  with  isobaricity  of  conditions)  the  intrinsic  thermo¬ 
diffusive  flows  of  components  will  prove  to  be  different  in  value. 

Second,  the  equality  of  bulk  flow  rates  at  the  inlet  and  outlet  is 
possible  also  in  the  case  when  both  components  diffuse  to  the  sane  side,  but 
at  different  rates,  and  the  hydrodynamic  flow  moves  in  the  opposite  direction, 
reducing  the  effectiveness  of  mixture  separation.  If  an ■ experiment  is 
carried  out  under  conditions  of  strict  isobaricity  (that  is,  in  the  absence 
of  hydrodynamic  flow),  the  bulk  flow  rate  at  t'.  e  outlet  from  the  upper 
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rate  at  the  outlet  from  the  upper  ("cold")  tubing  is  reduced.  In  this  case 
the  separation  will  be  at  a  maximum  and  will  be  governed  by  the  difference 
between  the  thermodiffusive  component  flows. 


The  quantity  of  component  equal  to  the  product  (c *  -  c)q  was  found  from 
data  in  Figure  4  aa  a  function  of  velocity  at  the  tubing  outlet.  Figure  5 
presents  this  function  for  the  light  component  in  the  upper  tubing  (curve  1) 
and  for  the  heavy  component  in  the  lower  tubing  (curve  2).  The  curves  have 
maxima  that  are  symmetrically  shifted  from  the  intersection  point.  The 
intersection  point,  as  we  see  from  the  figure,  corresponds  to  the  condition 
q  ■  qQ.  The  maxima  correspond  to  the  greatest  shift  of  concentration  between 
the  hot  and  cold  regions,  which  is  possible  only  in  the  absence  of  hydro- 
dynamic  flow  (that  is,  when  Ap  -  0) , for  any  assumption  of  the  direction  of 
thermodiffusive  flows.  From  the  k-axis  Intercept  of  the  maxima  it  is  clear 
that  the  velocity  at  the  output  from  the  lower  tubing  decreases  precisely  by 
as  much  as  the  velocity  rises  at  the  upper  tubing.  This  unambiguously  points 
to  the  same  direction  (from  the  cold  region  to  the  hot)  aa  the  thermo¬ 
diffusive  flows. 


Symbols 

J,j,  ■  thermo-self-dif fusion  flow;  Tq,  T^,  T2  ant*  T  ■  absolute  temper¬ 
atures  of  the  manometrlc  tube,  the  hot  and  cold  vessels,  and  the  mean 
temperature,  respectively;  nQ,  n^,  n2  and  n  -  number  of  molecules  per  unit 
volume,  respectively,  for  Tq,  T^,  Tj  and  X;  p  ■  pressure;  r,  S  and  L  *  radius, 
area  and  length  of  capillary;  SQ  ■  cross-sectional  area  of  manometric  tubing; 
n^,  n2  and  n  ■  coefficient  of  viscosity  at  the  temperatures  T^,  T2  and  T, 
respectively;  v  *  linear  velocity  of  drop;  f  -  force  of  friction; 
ot  ■  coefficient  of  proportionality  between  f  and  V;  k  -  Boltzmann  constant; 

W  ■  velocity  of  macroscopic  flow;  u  -  mean  velocity  of  thermal  agitation; 

X  and  X'  -  mean  free  paths  for  transport  of  numerical  density  and  momentum; 
c  ■  relative  concentration;  and  ■  coefficients  of  self-diffusion  and 
thermo-self-dif fusion;  qQ  and  q  ■  bulk  flow  rates  at  inlet  and  outlet  from 
tubing. 
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Figure  1.  Diagram  of  the  apparatus 


Figure  2.  Thermodiffusive  baro- effect  of  the  He-C02  mixture  as  a 
function  of  pressure 


|  I 

Figure  3.  Diagram  of  thennodif fusive  cell;  1,  2,  "Hot"  and 
ii cold"  tubings;  3,  Capillary 
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Figure  4. .  Variation  of  concentration  in  tubing  as  a  function  of 
change  in  velocity  at  outlet:  1,  Change  in  helium  concentration 
in  the  "hot"  tubing;  2,  Change  in  concentration  of  carbon  dioxide 
in  the  "cold"  tubing 


Figure  5.  Amount  of  absorbed  component  as  a  function  of  velocity 
at  outlet:  1,  For  helium  in  the  "hot"  tubing;  2,  For  carbon 
dioxide  in  the  "cold"  tubing 
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THERMAL  CONDUCTIVITY  OF  LIQUID  ORGANIC  COMPOUNDS 
G.  Kh.  Mukhamedzyanov  and  A.  G.  Usmanov 

A  review  of  methods  of  determining  coefficients  of  thermal  conductivity 
given  in  [1,  2]  points  to  the  lack  of  reliable  recommendations  for  precision 
calculation  of  X  with  temperature  change  for  a  wide  class  of  compounds.  At 
present  there' is  no  developed  system  of  concepts  which  would  allow  us  to  find 
coefficients  of  thermal  conductivity  even  within  the  limits  of  a  single 
homologous  series. 

Therefore,  we  conducted  measurements  of  X  of  various  classes  of  organic 
compounds  within  wide  limits  of  variation  in  number  of  carbon  atoms  in  the 
compound  molecule  n  and  over  wide  temperature  ranges,  on  the  basis  of  which 
we  propose  to  solve  the  following  problems: 

1)  to  establish  variation  of  X  as  a  function  of  temperature,  and  also  to 
discover  the  effect  of  various  functional  groups  on  thermal  conductivity; 

2)  to  derive  interpolation  equations  for  calculating  thermal  conduct¬ 
ivity  coefficients  of  the  series  of  compounds  studied. 

By  way  of  a  solution  of  these  problems  we  have  given  in  previous  n;.vV  :» 

results  of  investigations  on  the  thermal  conductivity  of  homologous  serioia  o,. 

saturated  hydrocarbons  normal-C^H^^j  (n  ■  6-24)  [3,  4];  normal  alcohol1. 

normal-C  H„  .,0H  (n  ■  1-18)  [5]}  saturated  monobasic  acids  C  H„  ,,COOH 
n  2n+l  n  2n+l 

(n  ■  1-18)  [61 }  complex  esters  of  saturated  monobasic  acids  ^nH2n+iC00CnH2n+l 
(n  ■  3-22)  [71;  simple,  mixed  esters  cnH2n+l0CnH2n+l  "  8-1^  [8]}  and 
individual  representatives  of  other  classes  of  compounds  [9]. 

This  study  is  a  continuation  of  this  series  of  measurements  and  includes 

experimental  studies  of  yet  another  series  of  classes  of  organic  liquids: 

unsaturated  hydrocarbons  of  the  olefin  series  CnH2n  (n  ■  6-14);  aldehydes 

C  H_  , .  CHO  (n  -  2-12);  ketones  CH,  ..COC  H_  (n  -  3-13);  aromatic 
n  Zn+1  n  zn+i  n  zn+i 


-554- 


C  H,  ,  (n  -  6-12)  and  halogen-substituted  hydrocarbons  C  H,  ,X  (n  -  2-9) 
n  i'n-o  n  Zuta 

where  X  -  I,  Br  and  Cl. 

Measurements  of  the  thermal  conductivity  of  these  chemically  pure 
liquids  were  also  conducted  at  atmospheric  pressure  by  the  most  carefully 
elaborated  method  of  the  heated  filament  [3].  Special  attention  was  paid  to 
eliminating  possible  effect  of  convection  on  the  value  of  X.  For  this 
purpose  the  experiments  were  conducted  In  a  measuring  tube  of  very  small 
diameter  (less  than  1  mm) .  The  results  of  X  measurements  for  different 
temperature  drops  in  the  layer  of  the  test  liquids  from  4  to  10®C  gave  good 
convergence . 

Lacking  the  opportunity  to  dwell  in  detail  on  the  method  of  measurement, 
we  note  only  that  In  the  experiments  possible  errors  were  taken  into  account. 
Controlled  measurements  of  X  for  several  compounds,  the  thermal  conductivity 
of  which  have  been  studied  by  many  authors,  revealed  that  the  deviations  were 
within  the  limits  of  1  percent.  "This  is  in  agreement  with  the  maximum 
error  of  our  experiment  under  the  most  unfavorable  conditions,  equal  to 
*1.5  percent. 

Tables  1-7  present  smoothed  experimental  data  on  the  organic  liquids  we 
investigated. 

The  thermal  conductivity  coefficients  of  particular  classes  of  compounds 
investigated  have  been  studied  by  a  number  of  authors.  Table  B  presents 
comparisons  of  the  results  of  our  measurements  of  thermal  conductivity 
coefficients  of  these  compounds  at  a  temperature  of  30*C,  X^q,  with  the  data 
of  other  authors,  X'^q.  The  X^g  values  taken  for  comparison  from  [10]  were 
obtained  as  the  result  of  averaging  experimental  data  of  not  less  than  two 
authors.  Deviation  of  our  data,  as  seen  in  Table  5,  does  not  exceed 
experimental  error.  Tables  1-7  show  that  the  thermal  conductivity  coeffi¬ 
cients  decrease  with  temperature  rise. 

The  nature  of  the  change  of  the  temperature  coefficient  of  thermal 
conductivity  a  for  normal  and  associated  liquids  has  been  established  lu 
[3-9]. 


t 


A  similar  picture  haa  been  observed  also  for  measurement  data.  The 
temperature  coefficient  a  for  a  given  homologous  series  depends  only  on  the 
number  of  carbon  atoms  in  the  compound  molecule.  A  decrease  of  a  with  rise 
in  n  has  been  observed  for  all  normal  compounds. 

Identical  values  oi  a  have  been  observed  for  the  ketones  K-CO-R'  apart 
from  any  dependence  on  distribution  in  the  hydrocarbon  radicals  R  and  R'  for 
equal  n. 


Table  1.  Thermal  Conductivity  of  Unsaturated  Olefinlc  Hydrocarbons 


Compound 

Chemical 

formula 

Temperature 
range  of 
measurement 
t,  °C 

Coefficients  of 
Thermal  conductivity, 
X,  w/m*deg 

Hexene  -  I 

C6H12 

0  +  60 

0.1260  +  0.1070 

Heptene  -  I 

C7H14 

0+80 

0.1290  +  0.1062 

Octene  -  I 

C8H16 

0  +120 

0.1316  +  0.0995 

Decene  -  1 

C10H20 

0  +160 

0.1363  +  0.0976 

Dodecene  -  I 

C12H24 

0  +200 

0.1402  +  0.0978 

Tetradencene  -  I 

C14H28 

0  +200 

0.1439  +  0.1063 

Table  2.  Thermal  Conductivity  of  Ketones 


Compound 

Chemical 

Temperature  range 

Coefficients  of 

formula 

of  measurement 

Thermal  conductivity, 

t.  °C 

X.  w/m*deK 

Acetone 

ch3coch3 

0 

+  40 

0.1678 

+ 

0.1541 

Methylethylketone 

ch3coch2Ch3 

0 

+  80 

0.1536 

+ 

0.1304 

Methylbutylketone 

ch3co(ch2)3ch3 

0 

+120  ' 

0.1452 

+ 

0.1164 

Methylamylketoue 

ch3co(ch2)4ch3 

0 

+140 

0.1432 

+ 

0.1100 

Methylhexylketouc 

CH3C0(CH2)5CH3 

0 

+  160  ! 

0.1410 

+ 

0.1044 

Diethylketone 

C2H5COC2H5 

0 

+  100 

0.1488 

+ 

0.1228 

Ethylbuty Ike tone 

c2h5co(ch2)2c2h5 

0 

+  140 

■ - - - f- - 

0.1426 

+ 

0.1096 

Dlamylketone 

C5H11C0C5H11 

20 

+  200 

0.1390 

+ 

0.1010 

Dihexylketone 

C6H13C0C6H13 

40 

+  200 

0.1396 

+ 

0.1078 

Methylpropylketone  CH^CO (CHj) 2CH3 

0 

+  100 

0.1486 

+ 

0.1226 

Table  3.  Thermal  Conductivity  of  Aldehydes 


Compound 

-  Chemical 
formula 

Temperature  range 
of  measurements 
t,  °C 

Thermal  conductivity 
coefficients , 

A,  w/m*deg 

Acetaldehyde 

ch3cho 

0+20 

0.1890  +  0.1804 

Butyraldehyde 

ch3(ch2)2cho 

0+60 

0.1551  +  0.1368 

Valeraldehyde 

ch3cch2)3cho 

0  +  100 

0.1460  +  0.1201 

Enantaldehyde 

ch3(ch2)5cho 

0  +  140 

0.1448  +  0.1120 

Decylaldehyde 

CH3(CH2)gCHO 

0  +  200 

0.1492  +  0.1092 

Lauraldehyde 

ch3(ch2)iScho 

60  +  200 

0.1454  +  0.1197 

Table  4.  Thermal  Conductivity  of  Bromine-Substituted  Saturated  Hydrocarbons 


Compound 

Chemical 

formula 

Temperature  range 
of  measurements 
t,  ®C 

Thermal  conductivity 
coefficients, 

X,  w/m*deg 

Ethyl  bromide 

C2H5Br 

0  +  30 

0.1046  +  0.0986 

Propyl' bromide 

C3H7Br 

0+70 

0.1062  +  0.0927 

Butyl  bromide 

C4H9Br 

0  +  100 

0.1080  +  0.0900 

Amyl  bromide 

C5HuBr 

0  +  120 

0.1096  +  0.0896 

Hexyl  bromide 

C6H13Br 

0  +  140 

0.1114  +  0.0891 

Nonyl  bromide 

C9H10Br 

0  +  200 

0.1161  +  0.0907 
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Table  5.  Thermal  Conductivity  of  Chlorina-Substi tuted  Saturated 

Hydrocarbons 


Compound 

Chemical 

formula 

Temperature  range 
of  measurements 
t,  °C 

Thermal  conductivity 
coefficients, 

X,  w/m*deg 

Ethyl  chloride 

c2h5ci 

0+10 

0.1264  +,  0.1242 

Propyl  chloride 

c3h7ci 

0+40 

0.1224  +  0.1128 

Amyl  chloride 

C5H11C1 

0  +  100 

0.1253  +  0.1048 

Amyl  iodied 

C6H13C1 

0  +  120 

0.1265  +  0.1041 

Hexyl  chloride 

C7H15C1 

0  +  140 

0.1277  +  0.1022 

Nonyl  chloride 

c9h19ci 

0  +  200“ 

0.1314  +  0.0980 

Table  6.  Thermal  Conductivity  of  Iodine-Substituted  Saturated  Hydrocarbons 


Compound 

Chemical 

formula 

Temperature  range 
of  measurements 
t,  °C 

Thermal  conductivity 
coefficients, 

X,  w/m'deg 

Ethyl  iodide 

w 

0  +  60 

0.0900  +  0.0792 

Propyl  chloride 

w 

0  +  100 

0.0916  +  0.0744 

Butyl  iodide 

Vs1 

0  +  120 

0.0933  +  0.0736 

Amyl  iodide 

C5H111 

0  +  140 

0.0956  +  0.0740 

Hexyl  iodide 

C6H13T 

0  +  180 

0.0978  +  0.0716 

Heptyl  iodide 

C7H15r 

0  +  200 

0.1001  +  0.0723 

Octyl  iodide 

0  +  200 

0.1029  +  0.0762 

Nonyl  iodide 

W 

0  +  200 

0.1056  +  0.0789 

» 


A  similar  picture  has  been  noted  earlier  also  for  homologous  series  of 
complex,  simple  and  mixed  esters  [7,  8J.  Maintenance  of  the  constancy  of 
temperature  coefficients  a  was  observed  also  for  identical  values  of  n  for 
different  halogen-substituted  hydrocarbons  CnH2n+^X  (X  ■  I,  Br  and  Cl). 
Naturally,  different  functional  groups  also'  affect  thermal  conductivity. 

The  presence  of  a  double  bond  (olefins)  and  also  replacement  of  a 
hydrogen  atom  by  the  halogens  I,  Br  and  Cl  leads  to  a  decrease  of  thermal 
conductivity.  In  the  latter  case  the  effect  ie  enhanced  with  rise  in 
molecular  weight  of  the  halogens. 
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A  decrease  ot  thermal  conductivity  was  noted  when  primary  and  secondary 
alcohols  are  oxidized  to  aldehydes  and  ketones.  This  can  be  easily  seen  by 
comparing  data  given  In  Tables  1-7  with  the  results  of  measurements  published 
earlier  [3-5]. 

The  experimental  material  obtained  on  thermal  conductivity  coefficients 
of  organic  liquids  allows  us  to  examine  the  possibility  of  using  different 
equations  in  calculating  X. 

We  compared  experimental  data  with  X  values  calculated  by  the  equations 
of  Weber  [11],  Smith  [12],  Palmer  [13],  Sakiadis  and  Coates  [14]  for  liquids 
representing  various  calsses  of  compounds. 

The  mean  error  of  calculations  using  the  above-cited  equations  lies 
within  the  limits  10-20  percent  and  often  rises  to  50  percent.  Similar 
deviations  have  been  noted  by  Reid  and  Sherwood  [1]  and  by  Brets’./ ayder  [2]. 


Table  9 


. —  ■“  -r— — — * -  *~  - .  i  . — - 1  —  '■  - — 

Smith  [12]  20  0.1320  0.1309  -  0.9  40  0.1508  0.1816  20.4 

80  0.1183  0.1352  14.3  80  0.1448  0.1837  27.0 


Palmer [13]  20  0.1320  0.1377  4.3  40  0.1508  0.1489  -1.3 

80  0.1163  0.1414  19.7  80  0.1448  0.1502  3.6 


Sakiadis  and  20  0.1320  0.1402  6.2  40  0.1408  0.1816  20.4 

Coates  £14]  80  0.1183  0.1539  30  80  0.1448  0.1782  23.3 


» 


In  Table  9,  by  way  of  example  a  comparison  Is  presented  of  our  exper¬ 
imental  data  on  ethylbenzene  and  normal-butyl  alcohol  [5]  with  data 
calculated  from  the  equations  In  [11-14], 

We  note  that  all  the  equations  considered  in  Table  9  require  knowledge 
of  precise  values  of  heat  capacity,  density  and  viscosity.  Determination  of 
the  latter,  in  turn,  Involves  experimental  difficulty  and  we  do  not  have 
reliable  data  for  most  of  the  compounds  of  the  classes  investigated. 

Therefore,  on  the  one  hand,  the  X  values  from  the  equations  in  [11-14] 
are  reliable  to  the  same  extent  as  are  the  values  of  heat  capacity,  density 
and  viscosity.  On  the  other  hand,  larger  deviations  are  evidently  due  to  the 
fact  that  these  methods  were  not  derived  by  a  reliable  theoretical  approach 
and,  as  was  noted  in  [15] ,  in  the  case  of  liquids  it  is  impossible  to 
establish  the  scope  of  their  application. 

A  method  based  on  similarity  considerations  was  used  in  [16,  17]  to 
determine  the  temperature  dependence  of  thermal  conductivity  coefficients. 

Its  application  does  not  require  introduction  of  any  corrections  both  in  the 
case  of  normal  and  associated  liquids,  as  well  as  their  mixtures  [17].  It  is 
pointed  out  in  [16,  17]  that  generalization  of  experimental  data  by  this 
method  in  the  coordinates  X/Xg,  and  S/S^  have  made  it  possible  to  obtain 
relative  thermal  conductivity  as  unique  functions  of  the  criterion  S/S^. 

For  normal  liquids  and  their  mixtures  this  ratio  is  of  the  form 

A  -  A,,  (l.'fS  *  0,7®  ■$]  ,  (1) 

for  associated  liquids  and  their  mixtures 

A  -  A,,  [  1.875  -  0,376^].  (2) 

In  Table  10,  by  way  of  example,  comparisons  are  given  of  experimental 
and  calculated  data  from  equations  (1)  and  (2)  for  individual  representatives 
of  the  classes  of  compounds  we  examined. 
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Table  10 


Liquids 


Chemical  Temper-' 

forumla  atrure 

t,°  C 


xliaiuml  conductivity 
coefficients 


X  .  X  , 

exp  calc 


40 

0.1305 

tmmwma 

-0.4 

Normal-decane 

oC10H22 

80 

160 

0.1209 

0.1020 

0.1209 

0.1028 

0.8 

80 

0.1506 

0.1506 

- 

Normal- te traces ane 

faC24H50 

160 

200 

0.1373 

0.1306 

0.1372 

0.1300 

-0.1 

-0.5 

46 

0.1363 

0.1362 

BBm 

Tetradecene-1 

C14H28 

120 

200 

0.1214 

0.1063 

0.1206 

0.1055 

1  1 

o  o 

0 

0.1520 

0.1518 

-0.1 

Ethyl  ester  of 

CH„C00C„He 

40 

0.1402 

0.1408 

0.4 

acetic  acid 

60 

0.1346 

0.1356 

0.7 

80 

6.1725 

0.1728 

0.1 

Normal-octadecyl  C, oH,_0H 

•  i  •  n  J.0  j/ 

alcohol 

120 

200 

0.1649 

0.1497 

0.1662 

0.1517 

0.8 

1.4 

0 

0.1678 

0.1678 

- 

Acetone 

ch3coch3 

20 

40 

0.1609 

0.1541 

0.1620 

0.1556 

-0.7 

-0.9 

0 

0.1504 

- 

Benzene 

C6H6 

40 

80 

0.1404 

0.1302 

0.1401 

0.1298 

-0.2 

-0.3 

Normal-butyric  acid 

C3H7C00H 

0 

80 

0.1454 

0.1356 

0.1450 

0.1355 

0.2 

160 

0.1260 

0.1270 

-0.8 

0 

0.1426 

0.1426 

Diethyl  eattjr 

C2H5OC2H5 

20 

0.1358 

0.1364 

-0.4 

Acetaldehyde 

CH3CH0 

0 

i  20 

0.1890 

0.1804 

0. 1899 
0.1804 

-0.5 

It  is  clear  from  Table  10  that  deviations  of  experimental  data  from  data 
calculated  using  equations  (1)  and  (2)  do  not  exceed  the  possible  exper¬ 
imental  error.  These  equations  can  be  recommended  for  determination  of 
thermal  conductivity  of  various  classes  of  organic  liquids  throughout  the 
range  of  the  liquid  state  at  normal  pressure.  This  requires  that  we  have 
data  on  X  even  though  for  just  one  temperature. 

We  have  al»o  applied  the  method  of  corresponding  states  [8,  17}.  The 
relationship  we  proposed 


A.  A*t(ut) 


(3) 


r 


for  calculation  of  thermal  conductivity  of  all  homologous  series  considered 
gives  good  agreement  with  experiment. 

However,  we  did  not  use  this  method  in  calculating  the  X  of  low 
representatives  of  the  series  in  [ 17 J .  Therefore,  use  of  equation  (3) 
requires  introduction  of  corrections  that  allow  for  the  effect  of  the  number 
of  carbon  atoms  n  for  low  representatives  of  a  series.  Additionally, 
equation  (3)  also  requires  knowledge  of  empirical  coefficients  that  take 
account  of  various  functional  groups  and  bonds,  and  for  associated 
liquids  —  the  effect  of  the  degree  of  association. 

This  fact  restricts  the  scope  of  application  of  (3)  in  determining  X  and 
places  obstacles  in  the  path  of  its  use  in  finding  the  thermal  conductivity 
coefficients  of  mixtures. 
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Footnote » 

559  Data  on  benzene,  toluene,  1,4-diiBethylbenzene  are  given  in 
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EXPERIMENTAL  STUDY  OF  THE  SORFT  EFFECT 


A.  M.  Narbekov  and  A.  G.  Usmanov 


Experimental  data  on  thermal  diffusion  in  the  liquid  phase  are  fully 
described  by  the  Soret  coefficient 


(1) 


which  is  a  measure  of  the  separation  effect.  Knowing  S  and  the  diffusion 
coefficient  D^t  we  can  calculate  other  thermodif fusion  characteristics  (a,. 

and  D^)  which  are  necessary  in  solving  several  problems  of  the  practical 
use  of  the  thermod if fusion  method  in  Industry. 

Experimental  study  of  the  Soret  effect  has  given  specific  material  to 
verify  the  theories  of  thermal  diffusion  in  liquids.  The  expression  for  the 
thermodiffusion  factor  a  Includes  the  so-called  heat  of  transfer  Q*^; 
however,  the  overall  function  for  Q*  is  lacking,  therefore  interpretation 
of  Q*  is  usually  made  from  the  results  of  measurement  of  the  Soret  effect. 

For  these  reasons,  obtaining  reliable  experimental  data  on  thermal 
diffusion  in  liquid  mixtures  has  remained  as  before  an  urgent  problem. 

Experimental  investigation  of  the  thermal  diffusion  process  in  liquids 
is  carried  out  in  cells  of  two  types.  The  first  type  of  cell  [1]  consists 
of  two  reservoirs  divided  by  a  porous  partition;  one  reservoir  is  heated,  the 
other  is  cooled,  and  a  temperature  gradient  is  Induced  in  the  partition 
and  through  it  occurs  molecular  diffusion.  In  order  to  sustain  the  same 
temperature  and  concentration  in  each  reservoir,  the  liquid  in  the  reservoirs 
is  stirred  by  magnetic  mixers.  The  cell  described  has  the  disadvantage  that 
it  is  difficult  to  maintain  uniform  temperature  and  concentration  in  the 
reservoirs  and  to  precisely  determine  the  temperature  difference. 
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The  second  type  of  cell  [2-4]  consisted  of  two  horizontal  plates,  between 
which  the  liquid  mixture  being  separated  resides.  Tha  horizontal  liquid 
layer  thus  obtained  is  subjected  to  the  action  of  a  vertical  temperature 
gradient.  The  ultimate  concentration  gradient  is  determined  optically  [2,  3, 
5-7],  conductometrically  [4,  8,  9],  or  by  sampling  at  different  levels  [10]. 

Usually,  in  cells  of  this  type  the  upper  plate  is  heated  and  the  lower 
cooled;  In  this  case  a  thermal  conductivity  regime  always  exists  in  the 
liquid  layer  regardless  of  the  value  of  the  temperature  difference  applied. 
Cells  of  the  second  type  are  Inapplicable  for  this  direction  (from  bottom  to 
top)  of  temperature  gradient,  if  the  heavier  component  of  the  binary  mixture 
accumulates  at  the  heated  wall,  since  reverse  mixing  of  the  liquid  takes 
place.  We  know  that  when  the  vertical  temperature  gradient  runs  in  the 
opposite  direction,  a  thermal  conductivity  regime  also  exists  in  the  liquid 
layer  under  specific  conditions.  For  experiments  on  thermal  diffusion  a 
reliable  thermal  conductivity  regime  in  the  horizontal  layer  of  liquid  is 
determined  by  the  following  value  of  the  dimensionless  product 

Gr*Pr  <  1500. 

The  permissible  temperature  difference  is  determined  from  the  equation 


Calculation  shows  that  temperature  difference  that  can  be  used  in  this 
case  is  slight  (for  example,  for  aqueous  solutions  if  the  cell  height  is 
0.5  cm,  the  temperature  difference  Is  about  1*C);  the  separation  effect  is 
very  slight,  that  is,  use  of  methods  of  study  sensitive  enough  for  the  small 
changes  in  concentration  of  the  mixture  investigated  is  necessary. 

Thus,  in  cells  with  a  vertical  temperature  gradient  any  temperature 
gradient  direction  can  be  used  depending  on  the  nature  of  the  separation. 
Further,  there  must  be  the  possibility  to  establish  concentration  changes 
in  the  course  of  thermal  diffusion  with  requisite  precision. 


This  requirement  can  be  satisfied  when  optical  methods  of  investigation 
are  used.  In  research  on  the  Soret  effect,  usually  two  optical  methods  are 
employed:  the  interference  method  and  the,  method  of  image  shifting.  We 
selected  the  interference  method  as  providing  the  possibility  of  a  relatively 
simple  quantitative  calculation  of  the  sought-for  characteristics  of  the 
process  under  study.  The  IAB-451  shadow  Instrument  was  used  as  the  interfer¬ 
ometer,  and  the  LG-75  optical  quantum  generator  served  as  the  light  source. 
Conversion  of  the  IAB-451  into  a  two-beam  diffraction  interferometer  is 
described  in  [11]:  transparent  diffraction  grids  R^  and  were  placed  in 
the  focal  planes  of  the  objectives  and  L2  (Figure  1,  a)  instead  of  a  slit 
and  a  knife,  and  in  the  plane  of  the  object  was  placed  the  diaphragm  D  with 
two  windows  used  to  form  the  two-beam  Interference  field.  A  beam  of  light 
from  the  optical  quantum  generator  was  focused  onto  the  plane  of  the  diffrac¬ 
tion  grid  R^  with  the  aid  of  the  short-focal  lens  L^. 

The  cell  for  investigating  thermal  diffusion  was  designed  (Figure  1,  b) 
to  consist  of  two  plates  (1)  arranged  horizontally  and  kept  at  different 
temperatures,  the  inserts  (2),  the  partition  (3),  and  two  flat-parallel 
plates  (4).  These  parts  formed  two  cells.  The  test  solution  was  poured  into 
one  of  the  cells,  and  the  solvent  into  the  other.  The  vessel  was  placed  in 
front  of  the  diagram  in  such  a  way  that  the  solution-containing  cell 
stood  in  front  of  one  of  its  windows,  and  the  solvent-containing  cell  in 
front  of  the  other.  Both  of  the  solutions  in  which  thermal  diffusion  took 
place  as  well  as  the  solvent  were  subjected  to  the  same  temperature  gradient, 
therefore  the  interference  picture  that  resulted  was  a  summation  of  two 
effects:  variation  of  the  index  of  refraction  owing  to  the  concentration 
gradient  established  upon  thermal  diffusion,  and  that  arising  from  the 
different  temperature  dependence  of  the  indices  of  refraction  of  the  solution 
and  of  the  solvent.  The  influence  of  the  second  effect  is  not  hard  to  take 
into  account  in  interpretation  of  the  interference  patterns. 

The  interference  patterns  were  obtained  both  when  the  interference  field 
was  adjusted  for  bands  of  infinite  width  (uniformly  illuminated  fields)  as 
well  as  for  adjustment  for  bands  of  finite  width.  When  the  bands  were 
adjusted  for  infinite  width  variation  of  the  index  of  refraction  associated 
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with  change  in  concentration  of  the  teat  «n1«Hnn  caused  interference 
bands  to  appear  in  the  pattern  (Figure  2,  b)  from  the  arrangement  of  which 
the  concentration  distribution  in  the  cell  could  be  calculated.  When  the 
interference  patterns  were  adjusted  fox  finite-width  bands  (Figure  2,  a) 
the  concentration  changes  could  be  calculated  from  the  band  shifts. 

The  interference  pattern  was  recorded  on  moving  motion  picture  film  K. 
The  interval  between  photographs  was  30  sec  and  was  provided  by  the  rotating 
disk  B. 

The  vessel  we  used  was  of  the  following  dimensions:  height  0.494  cm 
and  light  path  traveled  3.5  ciu.  The  lower  and  upper  plates  of  the  vessel 
were  kept  at  different  temperatures  with  the  aid  of  water  supplied  from 
thermostats.  Constancy  of  vessel  wall  temperature  was  controlled  by  copper- 
constantan  thermocouples,  the  signals  from  which  were  recorded  by  the 
R-307  potentiometer. 

Test  liquids  were  aqueous  solutions.  The  solutions  were  prepared  by 
weighing  on  analytical  balances;  the  index  of  refraction  was  determined  on 
the  IRF-23  refractometer. 

The  process  of  separation  was  carried  out  until  the  steady-state  was 
achieved.  The  variation  of  concentration  was  measured  from  the  beginning  of 
the  experiment  according  to  the  equation 

t difi/Oc.  (3) 

Variation  of  the  number  of  bands  (or  their  displacement)  was  determined  by 
measurement  on  the  MIR-12  comparator. 

We  must  note  the  following:  our  method  profitably  differs  from  the 
early  used  method  in  that  the  interference  pattern  obtained  more  fully 
reproduces  the  separatory  process.  The  interference  pattern  (Figure  2,  a) 
shows,  for  example,  that  the  distribution  of  concentration  by  height  in 
separation  is  uneven:  different  concentration  change  occurs  at  approximately 
the  central  line  of  the  cell. 


-570- 


Tiic  _eiuijei.ui.iue  difference  applied  on  the  layer  of  test  liquid  Is 
measured  optically.  For  this  purpose,  the  diaphragm  is  placed  in  such  a 
way  that  one  of  its  windows  proves  to  be  in  front  of. the  cell,  and  the  other 
remains  free,  and  the  beam  of  ].lght  passing  through  it  forms  a  comparison 
branch.  The  interference  pattern  obtained  here  consists  of  a  series  of 
horizontal  bands  characterizing  the  temperature  field.  The  temperature 
difference  is  determined  from  an  equation  analogous  to  (3): 


aT  * 


1  X 

TfijdT 


(4) 


The  Soret  coefficient  is  obtained  from  equation  (1).  The  values  of  S 
for  single-molar  solutions  of  KC1  and  KBr  (Table  1)  agree  well  with  the  data 
in  [2]s  for  KC1,  S  **  0.94*10  ^  1/deg,  and  for  KBr,  K  ■  1.33*10  ^  1/deg 
(AT  =  14°C) . 


Table  1.  Experimental  and  Calculated  Data  on  Determination  of  the 
Soret  Coefficient  for  Aqueous  Solutions  of  KCl  and  KBr 


System 

Concentr- 
ration  of 
second  com¬ 
ponent,  m 

1  . 

Temperature  °K 

Am*  102 

s.ioa 

1/deg 

Cold 

¥aU 

Temperature 

difference 

A  T 

Mean 

Temperature 
.  Ta™  . 

Water-KCl 

1 

290.5 

15 

298 

1,008 

1.53 

-1.02 

Water-KBr 

1 

290.5 

15 

296 

1.47 

2.12 

-1.41 

These  measurements  were  control  values  to  verify  the  method  used.  Later 
the  Soret  coefficients  were  determined  for  four  aqueous  solutions  (Table  2), 
for  two  of  which  (water-glucc>6e  and  water-propyl  alcohol)  literature  data 
are  lacking. 
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Thus,  by  using  the  IAB-451  shadow  instrument  employed  as  a  two-beam 
diffraction  interferometer,  the  Soret  coefficients  were  obtained  for  six 
aqueous  solutions  of  different  compounds.  The  interference  pattern  of  the 
separation  process  was  quite  easily  interpreted:  necessary  information  on 
kinetics  of  the  separation  process  can  be  obtained  from  the  nature  of 
variation  of  the  interference  pattern  with  time. 

Table  2.  Experimental  and  Calculated  Data  on  the  Determination  of  the 


Soret  Coefficient  for  Aqueous  Solutions  of  Sugars  and  Alcohols 


System 

Concentration 

Tempers 

kture  °K  I 

S'*  10  3 
1/deg 

of  second 
component,  c 

X 

Cold 

wall 

T1 

Differ¬ 

ence 

A  T 

Mean 

»gr 

av 

&l°3 

AC.102 

Water- 

saccharose 

Water- 

5 

289 

8 

293 

1.48 

3.5 

-0.92 

glucose 

Water- 

5 

294 

8 

298 

1.465 

2.4 

-0.63 

ethyl  alcohol 
Water- 

10 

288 

10 

293 

0.71 

33.75 

-3.72 

propyl  alcohol  10 

288 

10 

293 

0.925 

14.04 

-1.56 

Symbols 

S  ■  Soret  coefficient;  m  ■  molar  concentration;  T  ■  temperature; 

c  -  weight  concentration;  a  -  thermodif fusion  factor;  ■  thermodiffusion 

ratio;  DT  ■  thermodiffusion  coefficient;  "  coefficient  of  concentration 

diffusion;  Q*  ■  heat  of  transfer;  Gr  ■  Grashof  criterion;  Pr  *  Prandtl 
^  2 

criterion;  g  -  9.8  a/sec  ;  h  «  height  of  cell;  l  ■  extent  of  cell  along  light 
path  traveled;  p  ■  density;  0  ■  coefficient  of  volumetric  expansion; 
u  ■  coefficient  of  viscosity;  a  ■  coefficient  of  thermal  diffusivity; 
k  »  number  of  interference  bands;  X  «  wavelength  (6328*10  m) ; 


n  -  index  of  refraction;  the  index  "st"  Indicates  that  the  expression  in  the 
parentheses  was  taken  in  the  steady  state. 


Figure  1.  Diagram  of  apparatus:  a,  Optical  layout;  b,  Vessel 


1  2  S 


Figure  2.  Interference  patterns  of  the  thermal  diffusion  process 
of  an  aqueous  solution  of  normal  propyl  alcohol:  a,  Adjustment 
to  bands  of  finite  width  (1,  prior  to  applying  temperature 
difference;  2,  60  min  after  applying  temperature  difference; 

3,  after  completing  the  separatory  process) ;  b,  Adjustment  for 
bands  of  infinite  width  (1,  2,  at  different  time  intervals  from 
the  beginning  of  the  separatory  process;  3,  after  completing  the 
separatory  process) 
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EFFECT  OF  THERMOPHYSICAL  PROPERTIES  ON  MASS  WEAR  PARAMETERS  OF 

GLASSLIKE  MATERIALS 


Yu.  V.  Polezhayev  and  Yu.  V.  Khramov 

Glasslike  materials  (glass-plastics)  ere  characterized  by  the  fact  that 
exposed  to  a  high-temperature  gas  flow  they  disintegrate  with  the  formation 
of  a  melt  film,  the  viscosity  of  which  decreases  exponentially  with  rising 
temperature  T: 


This  leads  to  the  fusion  parameters  as  a  function  of  the  temperature 
distribution  within  the  body  and,  consequently,  to  the  fusion  parameters  as 
a  function  of  the  thermophysical  properties  of  the  condensed  phase. 

It  must  be  noted  that  direct  measurement  of  viscosity,  thermal  conduct¬ 
ivity.  heat  capacity  and  even  the  density  of  the  melt  at  temperatures  above 
2500°K  is  practically  impossible  for  a  long  series  of  reasons,  and  in 
calculations  several  extrapolated  values  as  a  rule  are  employed.  Though 
compared  to  viscosity  variation  of  the  remaining  properties  with  temperatures 
not  as  great  (such  that  they  need  not  be  taken  into  account  in  the  first 
approximation),  the  results  of  calculation  depend  materially  on  assumptions 
mad  relative  to  their  averaged  values. 

To  determine  this  effect  a  large  series  of  parametric  calculations  were 
made  in  which  the  function  of  the  mass  fusion  rate  G  and  the  surface 
temperature  were  studied  with  thermophysical  properties  varied  relative  to 
certain  standard  (marked  by  the  index  ()o  values).  Further,  density  and 
heat  capacity  varied  by  more  than  twofold,  thermal  conductivity  by  more  than 
fivefold,  but  the  exponent  in  the  viscosity  law  (1)  —  by  1.5-2  times. 


* 


Calculation  of  fusion  parameters  was  made  following  the  method  in  [1] 
applied  to  the  vicinity  of  the  critical  point  of  a  blunt  body  placed  in  an 
axially  symmetrical  air  current  with  the  following  retartlng  parameters: 
temperature  Tg  ■  6000°K,  pressure  pe  from  0.001  to  10  bar.  This  choice  of 
gas  flow  parameters  was  based  on  the  fact  that  the  value  of  the  retarding 
pressure  has  a  substantial  effect  on  the  ratio  of  the  quantity  of  material 
entrained  in  liquid  and  gaseous  form,  that  Is,  intensifies  or  weakens  the 
relationship  between  fusion  processes  and  heat  transfer  in  the  solid  phase. 
The  variation  of  the  retardingtemperature  by  2000°R  to  either  side  did  not 
lead  to  very  appreciable  rearrangements  in  the  nature  of  the  disintegration 
and  therefore  variation  of  this  parameter  is  not  of  Interest. 

Results  of  calculating  the  nonsteady-state  fusion  of  glasslike  material 
showed  that  of  all  the  physical  parameters  entering  into  the  differential 
equation  and  the  boundary  conditions,  the  course  of  the  nonsteady-state 
curve  G(t),  where  t  -  time,  is  affected  only  by  the  thermal  conductivity 
coefficient  and,  somewhat  more  weakly,  by  the  density.  All  the  remaining 
thermophysical  properties  mainly  determine  the  value  of  quasi-steady-state, 
settled  velocities  and  temperatures  of  disintegration.  Taking  account  of 
heat  capacity  C  >as  a  function  of  temperature  did  not  lead  to  differences 
either  in  the  instantaneous,  or  in  the  quasi-steady-state  values  of  the 
fusion  rate  compared  to  calculation  for  a  constant,  mean- integral  value 
over  the  entire  temperature  drop  in  the  body.  The  latter  fact  somewhat 
simplifies  engineering  methods  of  calculation. 

The  mass  fusion  rate  and  the  surface  temperature  T^  depend  very  weakly 
on  the  density.  However,  it  is  natural  that  this  parameter  has  a  more 
decisive  effect  on  calculation  of  linear  dimensions  of  the  entrained  layer. 

Thus,  two  parameters  remain:  thermal  conductivity  and  viscosity, 
knowledge  of  which  with  adequate  precision  represents  one  of  the  indispens¬ 
able  factors  of  calculation  reliability. 

Figure  1  shows  how  a  difference  in  the  value  of  the  thermal  conductivity 

\  from  the  standard  value  \  can  affect  the  relative  fusion  rate  G/G  , 

o  o 

where  G  corresponds  to  calculations  w.  th  X  ■*  \  .  The  different  curves 
o  o 

correspond  to  variation  of  pressure  of  the  retarding  flow  Pe  from  0.001  to 
10  bar. 
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Figures  2  and  3  present  similar  functions  for  the  characteristics  of  the 
effect  that  the  parameter  B  In  the  viscosity  low  (1)  has  on  the  relative 
values  of  the  fusion  rate  G/Gq  (Figure  2)  and  the  temperature  of  the 
disintegrating  surface  Tw/Tq  (Figure  3). 

It  is  interesting  to  note  that  not  one  of  the  thermophyslcal  coeffi¬ 
cients,  with  the  exception  of  viscosity,  has  an  influence  greater  than 
*3  percent  on  the  surface  temperature  T^,  and  therefore  this  fact  can  be 
used  in  estimating  thermophysical  properties  of  glass-plastics  at  very  high 
temperatures  (on  the  order  of  3000°K).  Actually,  if  in  comparing  the 
fusion  parameters  of  a  certain  glass-plastic  with  a  reference  material  (for 
example,  quartz  glass)  it  turns  out  that  the  surface  temperatures  are  close, 
and  the  wear  rates  differ,  the  difference  must  be  sought  for  in  the  nature  of 
the  variation  of  thermal  conductivity  or  in  the  value  of  the  heat  capacity. 
And,  in  contrast,  if  the  disintegration  surface  temperatures  differ  sharply 
for  the  glasslike  materials,  then  the  reason  for  this  lies  above  all  in  the 
difference  of  the  viscosity  laws  of  the  melts  and  cannot  be  explained  by 
other  factors. 

If  in  examination  nonBteady-state  characteristics  of  fusion  are 
additionally  brought  into  consideration,  by  comparing  parametric  calculations 
with  experimental  data  we  can  unambiguously  determine  at  least  two  out  of 
three  thermophysic ql  parameters:  viscosity  and  thermal  conductivity. 


different  pressures 


Figure  2.  Relative  fusion  rate  G/Gq  as  a  function  of  8/0  at  dif¬ 
ferent  pressures 


Figure  3.  Disintegration  surface  temperature  Tw/To  as  a  function 
of  8/ Bo  at  different  pressures 
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THERMO PHYSICAL  PROPERTIES  OF  MULTICOMPONENT  DISPERSE  SYSTEMS 
CONTAINING  DIAMOND  GRAINS 

A.  N,  Reznikov,  V.  K.  Lyakhov  and  Yu.  I.  Ivanov 


In  solving  problems  relating  to  diamond  treating  processes,  we  must  know 
the  thermophysical  properties  of  diamond  wheels  and  bars.  The  diamond¬ 
bearing  layer  of  a  wheel  or  a  diamond  bar  consists  of  disperse  mixtures  made 
up  of  several  components:  binder,  filler,  diamond  grains  and  air  (pores). 

The  relative  volume  of  diamond  grains  and  fllle  amounts  to  50  percent  In  the 
diamond-bearing  layer  of  wheels. 

This  study  Investigates  binders  on  a  bakellte  base  containing  as 
fillers  --  boron  carbide  and  iron  powder  —  and  wheels  with  bakellte  binder 
containing  diamond  grains.  Since  various  wheel  grades  differ  in  concentra¬ 
tion  of  components,  their  thermophysical  properties  vary  appreciably. 

Determination  of  thermophysical  characteristics  of  disperse  mixtures  is 
of  general  Interest,  but  it  involves  a  number  of  difficulties. 

Investigations  have  shown  that  the  bulk  weight  and  specific  heat 
capacity  of  mixtures  can  be  calculated  With  adequate  precision  based  on  the 
law  of  conservation  of  mass  and  the  additivity  rule: 


i 


C„ 


Tcm 


These  formulas  are  valid  if  in  the  mixing  process  no  new  chemical  compound  Is 

I 

formed. 

Several  theoretical  formulas  based  on  different  principles  are  available 
in  calculating  the  coefficient  of  thermal  conductivity  of  the  mixture  at  the 
present  time.  Calculations  based  on  these  formulas  give  results  that  differ 
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appreciably.  Therefore,  In  pat ui lei  wic’n  the  calculation  we  conducted 
experimental  studies  of  the  thermal  conductivity  coefficients  of  mixtures. 

Here  a  nonsteady-state  method  of  the  a-calorlmeter  and  the  steady-state  flat 
layer  method  were  used.  Both  binary  and  ternary  mixtures  were  investigated. 

For  binary  mixtures  bakelite-iron  powder  and  bakelite-boron  carbide,  five 
specimens  were  prepared  with  filler  concentration  in  the  range  0-50  percent. 

Graphs  (Figure  1,  a,  b)  were  plotted  from  experimental  results,  and  the 
graphs  simultaneously  presented  calculation  results  based  on  the 
V.  I.  Odelevskiy  formula  [1]  for  a  statistical  mixture  and  a  metrical  system  and 
results  from  calculations  using  the  K.  Likhteneker  formula  [3]. 

Results  show  that  within  the  limits  of  the  concentrations  examined  the 
formulas  of  other  authors  [2-3],  etc.)  for  a  metrical  system  give  approx¬ 
imately  identical  results,  which  within  the  limits  of  additive  concentration 
of  30-50  percent  give  values  lower  than  the  experimental  data.  The  best 
agreement  with  experiments  iB  given  by  the  K.  Likhteneker  formula 


(1) 


Experiments  with  the  ternary  mixture  bakelite-boron  carbide-iron  were 
also  conducted  with  specimens  of  different  concentrations.  Experimental 
results  are  shown  in  a  graph  (Figure  1,  c).  The  iron  powder  concentration  in 
the  specimens  was  varied  in  such  a  way  that  the  overall  content  of  iron 
powder  and  boron  carbide  grains  remained  constant  at  50  percent.  This 
variation  in  concentrations  was  characteristic  for  compositions  of  various 
binders  of  diamond  wheels.  In  the  graph  (cf.  Figure  1,  c)  experimental  data 
were  compared  with  the  results  of  calculation  according  to  the  K.  Likhteneker 
formula  (1). 

Calculations  based  on  formulas  for  the  metrical  system  and  a  statistical 
mixture  give  results  that  differ  appreciably  from  experimental  data. 

Experiments  with  the  ternary  system  bakelite-boron  carbide-diamond 
powder  were  conducted  with  specimens  (diamond  wheels)  of  two  concentrations. 
The  results  of  these  experiments  and  their  comparison  with  calculation  based 
on  the  K.  Likhteneker  formula  are  shown  in  the  table. 


Table 


No. 

Wheel  Characteristics 

Composition  in  £  of 
Bakelite  carB¥de 

volume 

Diamond 

-1 

mw*  cm 

deg-1 

■  Calculated 

1 

AS40-BI-50 

50 

37.5 

12.5 

15.6 

17.3 

2 

AS12-BT-50 

50 

37.5 

12.5 

17.1 

17.3 

3 

AS4-BI-50 

50 

37.5 

12.5 

19.8 

17.3 

4 

AS12-BI-100 

50 

25.0 

25.0 

19.5 

22.1 

Specimens  1,  2  and  3  differed  by  diamond  and  boron  carbide  grain  sizes. 
This  was  accounted  for,  in  particular,  by  deviation  in  experimental  data  of 
the  first  three  specimens  in  which  the  percentage  composition  of  the  compon¬ 
ents  was  identical.  Specimens  2  and  4  differed  in  diamond  and  boron  carbide 
content .  Since  the  percentage  content  of  diamond  in  the  fourth  specimen  la 
greater  than  in  the  second,  and  the  thermal  conductivity  coefficient  of 
diamond  is  greater  than  the  thermal  conductivity  coefficient  of  boron 
carbide,  with  rise  in  diamond  concentration  and  with  retention  of  the  overall 
diamond  and  boron  carbide  concentration,  the  equivalent  thermal  conductivity 
coefficient  rises. 

Calculation  of  the  equivalent  thermal  conductivity  coefficient  for  these 
tables  was  made  with  a  five  percent  air  pore  content  taken  into  account  (the 
table  lists  the  percentage  content  of  components  without  including  porosity). 

Conclusions 

1.  It  follows  from  our  studies  that  calculations  according  to  the 
K.  Likhteneker  formula  O )  give  satisfactory  agreement  with  experimental 
data. 

2.  Calculations  based  on  the  formulas  of  other  authors,  both  for  a 
statistical  mixture,  as  well  as  for  a  metrical  system,  give  results  that 
differ  considerably  from  experimental  data. 

Only  in  the  low-concentratiou  region  (Figure  1)  is  the  deviation 
between  experimental  data  and  calculations  based  on  the  formulas  of  various 
authors  small. 
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Symbols 


Ymlx*  cmlx‘  ^3  "  ^ulk  w®ight,  specific  hest  and  equivalent  thermal 
conductivity  coefficient  of  tha  mixture;  Vj  »  and  \  »  as  above,  for 

indivldual  components;  p^  and  p2  ■  bulk  concentrations  of  components. 
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Figure  1.  Results  of  comparison  of  experimental  and  calculated 
data:  1,  Experimental  results;  2,  Calculation  based  on  the 
K.  Likhteneker  formula  (1);  3,  Calculation  based  on  the 
V.  I.  Odelevskiy  formula  [1]  for  a  statistical  mixture; 

4,  Calculation  based  on  the  V.  I.  Odelevskiy  formula  [1]  for  a 
matrlcal  system;  a)  Bakelite-iron  powder  mixture;  b)  Bakelite- 
boron  carbide  mixture;  c)  Bakelite-iron  powder-boron  carbide 
mixture 
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STATISTICAL  DETERMINATION  OF  THE  DIFFUSION  COEFFICIENT 
L.  A.  Rott,  V.  B.  Nemtsov  and  V.  S.  Vikhrenko 

The  possibility  of  expressing  kinetic  coefficients  of  a  condensed  system 
by  autocorrelation  functions  is  of  unquestioned  and  stimulating  interest  in 
the  statistical  theory  of  nonequilibrium  processes. 

This  approach  was  used  in  calculating  diffusion  coefficients,  which  is 
the  main  problem  area  of  the  present  study.  But  first  of  all  we  will 
consider  how  we  cun  express  viscous  properties  of  the  medium  in  the  language 
of  autocorrelation  functions  (using  the  well  known  Kubo  method  [1]).  It  is 
precisely  an  expression  for  the  coefficient  of  shear  viacoalty  that  la 
essential  for  a  definitive  calculation  of  the  diffusion  coefficient. 

First  we  will  examine  the  general  case  of  an  anisotropic  homogeneous 
medium.  Let  a  weak  mechanical  perturbation  be  applied  on  the  Isolated  system 
in  which  its  Hamiltonian  H  Is  varied  by  AH  ■  -AF(t).  Then  the  mean  value 
AB(t)  of  the  variation  in  the  dynamic  variable  B(t)  is  determined  according 
to  [1]  with  the  aid  of  the  expression 

<1) 

where  0  ■  kT,  B  dB/dt,  and  the  brackets  denote  the  averaging  over  the 
equilibrium  canonical  assemblage. 

To  explore  viscous  properties  of  the  system  we  will  examine  Its  behavior 
when  a  weak  deformation  is  applied 
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STATISTICAL  DETERMINATION  OF  THE  DIFFUSION  COEFFICIENT 
L.  A.  Rott,  V.  B.  Nemtsov  and  V.  S.  Vikhrenko 

The  possibility  of  expressing  kinetic  coefficients  of  a  condensed  system 
by  autocorrelation  functions  is  of  unquestioned  and  stimulating  interest  in 
the  statistical  theory  of  nonequillbrium  processes. 

This  approach  was  used  in  calculating  diffusion  coefficiente,  which  is 
the  main  problem  area  of  the  present  study.  But  first  of  all  we  will 
consider  how  we  can  express  viscous  properties  of  the  medium  In  the  language 
of  autocorrelation  functions  (using  the  well  known  Kubo  method  [1]).  It  Is 
precisely  an  expression  for  the  coefficient  of  .ihear  viscosity  that  is 
essential  for  a  definitive  calculation  of  the  diffusion  coefficient. 

First  we  will  examine  the  general  case  of  an  anisotropic  homogeneous 
medium.  Let  a  weak  mechanical  perturbation  be  applied  on  the  isolated  system 
in  which  its  Hamiltonian  H  is  varied  by  AH  ■  -AF(t).  Then  the  mean  value 
AB(t)  of  the  variation  in  the  dynamic  variable  B(t)  is  determined  according 
to  [1]  with  the  aid  of  the  expression  ... 

a  £>(t)  -  B'j<  .  (1) 

where  6  ■  kT,  B  ^  dB/dt,  and  the  brackets  denote  the  averaging  over  the 
equilibrium  canonical  assemblage. 

To  explore  viecous  properties  of  the  system  we  will  examine  its  behavior 
when  a  weak  deformation  is  applied 

X*  (2) 
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the  deformation,  and  is  after  the  deformation;  is  the  component  vector 
of  particle  displacement.  The  particle  momenta  here  are  transformed 


according  to  the  law  of  canonical  transformations  in  the  Lagrange  formalism: 


(3) 


We  note  that  the  general  approach  developed  here  differs  from  the 
particular  case  based  on  the  use  of  the  Feinwan  ideas  [1,  2], 

With  a  precision  that  includes  neglecting  the  squared  values  and  the 
higher  powers  of  deformation,  we  can  write 


(4) 


For  a  system  of  N  particles  with  pairwise  interaction  described  by  the 
lntermoleculap  potential  $(r),  the  Hamiltonian  becomes 


(5) 


Then  la  determined  by  the  following  equality  with  the 

above-indicated  degree  of  precision 


i » 


(6) 


where 


-f**' 


ht 


Taking  into  account  the  symmetry  of  the  tdnsor  n^,  we  will  represent 
the  expression  for  dH  as 


-588- 


4 


(7) 


where 


deformation  tensor. 
Thus, 


AU.fi> 

—  *  ♦  i  rmn  U/nn  > 


'ft*mhrnn  t  fi •  Umn 

Selecting 

8*  fnt«(s)ds  , 

we  obtain  based  on  the  relationship  (1) 

■aB*®  J *  flmn  (o)niK(i)>umn(t-sjcis 


(8) 


(9) 


Let  the  deformation  be  cyclic 


'  Umn*Umn(o)t 


iult 


Then 


(10) 


Ab  *  e’iui<  (lmAMni«U)>dt]  (11) 


We  will  compare  this  equality  wijth  the  well  known  phenomenological 
relationship 


dS*  ^tMmn  Umn  (&*  mt^ininn  UmnJ  > 


(12) 
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we  obtain  an  expression  for  the  tensor  of  viscosity  coefficients 


q(ttmn  m  &''Jemiwt<  /?«,  CO)  f)i„  (t)  >dt  (13) 

It  is  essential  to  note  that  the  flow  II  Is  determined  with  a  precision 

inn 

up  to  any  value,  the  divergence  of  which  equals  zero.  To  eliminate  this 
indeterminacy,  we  can  according  to  [3]  rewrite  the  expression  (13)  in  the 
following  form  (the  canonical  assemblage  is  used  for  the  averaging): 

(ftnviCO)  -  m)  (hiK(i)~h*M)>dt  .  <U> 

ar 

where 


Here,  is  the  equilibrium  mean  stress  tensor  and  E  Is  the  mean  value 
of  total  energy  of  system  H^. 

The  expressions  presented  agree  with  the  results  obtained  in  [4,  5)  by 
another  method. 

Fo  .  the  isotropic  liquid  a  familiar  result  for  shear  viscosity  [1]  is 
obtained  from  the  general  expression 

?(«)  *  co)htiCt) >cct  as) 

0 

The  case  of  slow  immersion  corresponds  to  the  condition  co  ■  0. 

The  expression  for  bulk  viscosity  also  follows  from  (14). 

For  the  isotropic  liquid,  if  we  use  the  Stokes  model,  we  can  determine 
the  coefficient  of  friction  of  a  molecule  C  via,  the  coefficient  of  shear 
viscosity  n.  At  the  same  time  the  coefficient  €,  according  to  Kirkwood  [6], 
can  also  be  represented  by  the  integral  of  the  autocorrelation  function 


(16) 


j  *  3kT I  <  K (°JK  ($)»ds 

oJ 

Here  K  Is  the  force  acting  at  moment,  s  on  a  selected  molecule  of  the  liquid 
on  the  part  of  the  remaining  molecules  In  the  system.  Since  the  system  Is 
isotropic,  then  the  integral  encloses  a  scalar  product. 

The  self- diffusion  coefficient  (ve  are  considering  a  single-component 
system)  can  be  expressed  by  the  integral  of  the  autocorrelation  function  of 
molecular  velocity,  and  in  the  simplest  model  is  equal  to  It  [1],  that  Is, 


T>‘~ij~?toj?(slicls  <17> 

Here  we  also  have  a  scalar  product  of  vectors. 

In  all  the  cases  cited,  attainment  of  ultimate  numerical  results 
involves  great  and  as  yet  uneliminated  difficulties. 

Though  integral  equations  determining  autocorrelation  functions  have 
been  obtained  (cf.,  for  example  [7,  8]),  the  solution  of  these  equations 
sf1 1 1  remains  unclear.  Therefore  to  some  extent  we  can  surmount  the 
difficulties  by  attempting  to  represent  the  integral  of  the  autocorrelation 
function  as  a  product  of  the  mean  square  of  the  quantity  under  consideration 
and  relating  to'  a  single  arbitrary  moment  of  time  for  its  mean  relaxation 
time. 

Here  we  must  forthwith  point  to  the  considerable  difference  in  relax¬ 
ation  characteristics  for  quantities  dependent  on  spatial  coordinates 
(forces)  and  quantities  dependent  on  velocity.  Thus  the  coefficient  of 
friction  of  an  individual  molecule  can  be  represented  (for  :the  system  states 
far  removed  from  the  critical  point)  as 


(1«) 


We  let  represent  the  mean  relaxation  time  for  the  force  K  dependent  only 
on  spatial  coordinates. 
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The  Integral  for  the  diffusion  coefficient  can  be  written  In  approximate 


terms  as  follows 


i?vz f 


Here  t  la  the  mean  relaxation  time  of  molecular  velocity  (momentum). 

P 

Using  the  concept  of  relaxation  times  t  and  t  ,  we  can  determine  in  an 

*1  P 

approximate  way  the  autocorrelation  function  entering  into  the  formula  for 
the  shear  viscosity  coefficient  (15)  when  w  ■  0.  The  members  of  the 
expression  for  JI^  relates  here  to  a  single  moment  of  time  and  the  averaged 
quantity  can  be  reduced  to  binomial  form  with  the  aid  of  the  procedure 
described  by  Zwanzig  19]. 

The  final  resul'.  is 


Thus,  the  problem  of  determining  kinetic  characteristics  breaks  down 

into  two  subproblems.  One  is  associated  with  the  need  to  carry  out  the 

procedure  of  statistical  averaging  in  the  equilibrium  state  of  the  system, 

and  the  other  with  determination  of  relaxation  time,  and  the  relaxation  times 

t  and  t  differ  from  each  other  by  two  orders  of  magnitude. 

9  P 

The  statistical  determination  of  relaxation  time  t  (for  simple  liquids 
-12  " 

it  is  of  the  order  of  10  sec)  has  been  given  in  [10] .  The  values  of  t 
agree  with  thermal  estimates  [11].  Here  however  we  will  dwell  on  determin¬ 
ation  of  Tp.  For  this  purpose  we  will  use  the  Langevin  equation,  the 
statistical  approach  of  which  was  initially  given  in  the  already  cited 
study  by  Kirkwood.  If  we  do  not  take  into  account  the  member  expressing  the 
rapidly  oscillating  (random)  force,  the  Langevin  equation  becomes 
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where  p  Is  particle  momentum,  m  Is  particle  mass;  £  is  equal  to  the  above¬ 
given  expression  (16).  From  this  it  is  clear  that  the  characteristic 
momentum  relaxation  time  will  be  ^2) 

r/>*TL 

Following  the  general  Kirkwood  approach,  we  can  repeat  the  derivation  of 
the  Langevin  equation,  but  now  with  use  of  the  two-index  distribution 
functions  (statistical  method  of  conventional  distributions).  The  friction 
coefficient  for  the  molecule  will  then  be  equal,  according  to  (18),  if  we 
carry  out  the  statistical  averaging  procedure  [11],  to  the  following: 


!  °  A*fv~J  fo'(r)]  *V(r)  r*dr 


Here  4>(r)  is  the  pairwise  intermolecular  potential  and  v  is  the  molecular 

3 

volume  (4/3  irr  ■  V). 

o 


fUJ  »  V 


We  will  recall  that  in  the  statistical  scheme  we  used  denotes  the 

probability  density  that  in  the  two  selected  m  lecular  volumes  v1  and  v, 

1  2  1  * 

correspondingly  surrounding  the  coordinates  q  cv^  and  q  cvj  any  two 

molecules  are  present,  and  in  the  remaining  volumes  there  are  also  no  more 

than  one  molecule  (this  approximation  is  called  the  approximation).  The 

probability  density  is  similarly  determined  for  a  single  molecule  F^(q). 

From  (24)  it  is  clear  that  $(r)  designates  the  conventional  distribution 
2  1  1 

function  F^(q  /q  )•  Below,  the  fixed  point  q  is  chosen  in  the  center  of 
a  sphere  having  a  volume  V. 

This  approach  was  used  in  [11]  for  calculating  the  friction  coefficient 
5,  and  then  also  the  shear  viscosity  coefficient  for  simple  nonpolar  liquids 
on  the  phase  transformation  line  ultimately  equal  to 


'•fs\ {¥^rv*,0~' 


g  -1  -3 

g*cm  • nee 
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where  m  is  Che  molecular  weight,  V  la  the  volume  of  a  single  mole  In  era, 
e/k  and  o  are  parameters  of  the  Lennard- Jones  potential  (c  given  In 
Angstroms  and  e/k  in  degrees). 

Using  the  determination  of  according  to  (22)  with  (23)  :aken  into 
accbunt,  we  will  find  the  final  expression  for  the  self-diffusion  coefficient 
of  a  single-component  system.  Since  the  mean  square  of  the  velocity  is 
3kT ,  then 


or,  using  the  expression  for  the  shear  viscosity  coefficient  (25), 


(O'*  CH*C*K' 


Below  are  presented  results  of  calculating  the  Belf-dif fusion  coeffi¬ 
cient  according  to  (27)  for  a  number  of  compounds  (argon,  hydrogen, 
nitrogen  and  carbon  dioxide)  along  the  phase  transformation  line  (we 

4 

individually  also  present  the  results  of  calculating  the  momentum  relaxation 
time). 


Compound 


D  1O0 


X  10l4sec 
P 


An 

2 

An 

2 

An 

2 

An 

2 

An 

2| 

Hi 

2i 

Hi 

2, 

/Vi 

2, 

M 

2, 

ut  no  little  interest  is  the  tact  that  the  relaxation  time  r  and  ▼  ara 

P  ^ 

Inversely  related,  since,  by  combining  (18)  and  (22),  we  can  write 


It  appears  to  the  authors  that  this  fact  la  very  significant  in 
accounting  for  characteristics  in  the  behavior  of  a  compound  In  its  critical 
state,  when  diffusion,  just  like  other  properties,  is  manifest  In  a  specific 
way  (cf.,  for  example  [12-13]). 

Taking  into  account  the  approximations!  nature  of  the  expression  (27) 
adopted,  comparison  with  experimental  data  [1]  must  be  regarded  as  wholly 
satisfactory. 

By  knowing  the  coefficient  of  friction  of  a  molecule  we  can  determine 
the  coefficient  of  molecule  mobility  b. 

Phenomenologically,  the  mobility  coefficient  is  defined  as  the 
coefficient  of  proportionality  between  the  diffusion  flow  and  the  thermo¬ 
dynamic  driving  force  (the  gradient  of  the  chemical  potential).  Corre¬ 
spondingly,  the  mobility  coefficient  can  be  defined  via  the  diffusion 
•coefficient 


where  p  is 
If  we 


the  pressure. 

take  the  self-diffusion  coefficient  as  the  coefficient  D,  then 


6* 


(20) 


The  transition  to  the  binary  isotropic  system  is  based  on  use, 
appropriately  enough,  of  a  system  of  two  Langevin  equations  (for  the  mole¬ 
cules  of  the  a  and  b  species).  The  friction  coefficients  were  determined 
with  the  aid  of  the  sum  of  two  autocorrelation  functions 
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The  expression  (14)  for  Che  tensor  of  Che  viscosity  coefficients  is 
easily  generalised  for  the  case  of  a  binary  mixture.  In  this  case  the 

A 

operator  becomes: 


■"*?+ 


I hi  >•/ 


(32) 


Here  is  the  number  of  molecules  of  species  a,  n2  is  the  number  of 
molecules  of  species  b,  and  4^  la  the  potential  of  pairwise  molecular 
interaction  of  the  corresponding  species. 

The  averaging  leading  to  an  expression  analogous  to  (20)  is  conducted 
here  with  the  aid  of  two  momentum  relaxation  times  of  the  corresponding 
5a  and  5b. 
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EFFECT  OF  THE  COMPOSITION  OF  CLOTH  AND  THE  FORM  OF  BINDING 
OF  ABSORBED  MOISTURE  ON  ITS  THERMOPHYSICAL  CHARACTERISTICS 

N.  I.  Salivon  and  M.  F.  Kazanskiy 

The  main  Chermophy3lcal  characteristics  of  cloth  determining  their  heat- 
protective  properties  is  thermal  conductivity  X,  thermal  dlffusivlty  a  and 
heat  capacity  c. 

Since  the  values  of  these  characteristics  depend  decisively  on  moisture 
content  bound  to  the  solid  phase  of  the  cloth  in  different  ways,  we 
conducted  a  study  aimed  at  discovering  how  moisture  in  different  forms  and 
types  of  binding  with  cloth  of  different  physicochemical  character  affects 
the  thermophysical  properties  of  the  cloth. 

The  objects  of  the  study  were  specimens  of  100  percent  lavsan  cloth  and 
woolen  combed  suit  cloth  of  No.  32/2  weave  in  the  base  and  in  the  neck,  with 
Interweaving  of  2/2  serge  for  an  almost  identical  density  and  weight  with 
10,  30  and  50  percent  lavsan  content,  and  also  a  specimen  of  cloth  of  the 
same  structure  made  of  100  percent  wool,  woven  in  the  Central  Scientific 
Research  Institute  of  the  Garment  Industry. 

The  water-retaining.  properties  of  the  fabrics  were  determined  by  means 
of  Isothermal  drying  thermograms  [1]  and  with  the  aid  of  water  vapor 
adsorption  isotherms.  The  results  of  analysis  of  the  forms  in  which  water  is 
bound  with  these  fabrics  are  shown  in  the  Table. 

Thermophysical  coefficients  X,  a  and  c  were  determined  on  an  apparatus 
we  developed  [2]  using  the  method  of  two  temperature-time  points  [3]  for 
specimens  2~3  mm  thick. 

The  results  of  a  series  of  experiments  with  each  specimen  are  shown  in 
Figures  1  and  2  in  the  form  of  curves  describing  the  thermal  coefficients  of 
functions  of  moisture  content.  In  the  figures  the  vertical  dashed  lines 
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denote  the  boundaries  of  moisture  of  various  forms  and  types  of  binding  with 
fabric:  A  is  the  maximum  amount  of  adsorbed  moisture;  T  is  the  moisture  of 
maximum  hygroscopic  state. 


Table.  Amount  of  Water  in  Various  Forms  of  Binding  with  Fabric 


Speci¬ 

^content  of  lavsan 

Moisture  in 

Maximum  amount 

Amount  of 

men 

in  mixture  with 

mono-mole  cular 

of  absorbed 

hygroscopic 

No. 

wool 

adsorption  % 

moisture 

mois  ture 

1 

0 

5.8 

9.6 

28.5 

2 

10 

5.4 

8.7 

24.6 

3 

30 

4.0 

6.7 

17.4 

4 

50 

2.7 

5.1 

15.0 

5 

100 

- 

0.3 

3.8 

At  low  moisture  contents  (Figure  1,  curves  1),  when  only  adsorbed 
moisture  is  found  in  the  specimen,  the  coefficient  of  thermal  conductivity  X 
remains  almost  unchanged.  An  appreciable  rise  in  the  thermal  conductivity 
coefficient  is  observed  in  the  moisture  content  interval  from  the  adsorbed  to 
the  maximum  hygroscopic  moisture  content,  when  micropores  fill  up  with 
moisture  and  thermal  contacts  between  fibers  are  established  via  formation 
of  continuous  water  films.  Upon  further  rise  in  moisture  content,  the 
increase  in  thermal  conductivity  coefficient  is  less  intense  and  is  char¬ 
acterized  by  a  tendency  toward  saturation. 

The  coefficient  of  thermal  diffusivity  a  also  remains  almost  unchanged 
at  low  moisture  contents  (Figure  2,  curve  III),  The  appreciable  rise  in  the 
thermal  diffusivity  coefficient  begins  from  the  moment  that  capillary 
condensation  sets  in,  when  micropores  of  radius  r  <  10  ^  cm  begin  to  be 
filled  in  With  wnter,  An  intensive  rise  in  the  thermal  diffusivity  coeffi¬ 
cient  with  rise  in  moisture  content  is  observed  all  the  way  to  attainment  of 
maximum  hygroscopic  state  of  the  specimen. 

The  thermal  diffusivity  maximum  for  all  specimens  corresponds  to  the 
moisture  content  close  to  the  maximum  hygroscopic  state,  when  all  specimen 
pores  are  filled  with  moisture,  and  when  water  films  are  formed  between 
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Individual  fibers  producing  theiw.1  contacts.  This  transition  from  the 
three-component  system  (fabric-water-air)  to  a  two-component  system  (fabric- 
water)  is  characterized  by  a  drop  in  the  thermal  diffusivity  coefficient. 

Due  to  the  diversity  of  forms  and  types  of  binding  of  fabric-absorbed 
water,  in  each  of  the  curves  the  functions  of  the  thermal  conductivity  and 
thermal  diffusivity  coefficients  can  be  divided  into  three  sections.  Here 
the  value  of  the  linear  initial  section  corresponding  to  the  quantity  of 
adsorbed  moisture  is  reduced  with  rise  in  lavsan  addition  both  on  the 
thermal  conductivity  curve  and  on  the  thermal  diffusivity  curve,  but  for  the 
specimen  of  pure  lavsan  there  is  no  such  section.  The  maximum  on  the  thermal 
diffusivity  curve  is  shifted  to  the  side  of  lower  moisture  content  with  rise 
in  percentage  lavsan  content,  this  shift  corresponding  to  the  lower  moisture 
capa-ity  of  lavsan  compared  to  wool. 

The  heat  capacity  for  all  specimens  rises  linearly  with  riBe  in  moisture 
content  (curve  II). 

Just  as  for  typical  capillary-porous  and  colloidal  capillary-porous 
bodies,  the  thermophysical  coefficients  as  functions  of  moisture  content 
for  fabrics  with  different  lavsan  content  are  determined  above  all  by  the 
diversity  of  forms  and  types  of  bonding  of  fabric-absorbed  moisture.  In 
this  case  this  diversity  depends  on  the  amount  of  lavsan  added  to  the  wool. 
These  conclusions  on  the  effect  of  lavsan  addition  on  thermophysical  coeffi¬ 
cients  of  fabrics  can  be  used  in  designing  suits  with  desired  heat-protective 
properties. 

Symbols 

A  -  thermal  conductivity  coefficient;  a  ■  thermal  diffusivity  coeffi¬ 
cient;  c  ,■  specific  heat  coefficient;  r  -  capillary  radius;  A  ■  boundary  of 
maximum  amount  of  adsorbed  moisture;  T  ■  boundary  of  moisture  in  maximum 
hygroscopic  state. 
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EFFECTS  OF  FILLING  ON  THERMAL  CONDUCTIVITY,  ELECTROCONDUCTIVITY 
AND  DIELECTRIC  CONSTANT  OF  THIN-LAYER  POLYMER  COATINGS 

F.  M.  Smekhov,  L.  N.  Novichenok, 

L.  V.  Nitsberg  and  A.  I.  Nepomnyashchly 

Variation  of  thermophysical  and  electrical  properties  of  filled  polymer 
coatings  upon  variation  of  filler  content  was  investigated. 

Epoxide  resin  E-SSCM^  ~  1000)  solidified  by  the  FO-200  polyamide  resin 
(M^  ~  3000)  with  formation  of  a  reticulate  structure  was  used  as  the  film¬ 
forming  agent. 

The  rutile  modification  of  TiOj  of  the  R-l  grade  served  as  the  filler, 
and  its  concentration  in  the  binder  varied  from  10  to  50  percent  by  volume. 

Experimental  Method 

The  thermal  conductivity  of  coatings  X  was  measured  by  a  nonsteady-state 
method  for  a  single  temperature-time  point  in  the  specimen  [1]  on  an 
instrument  at  the  Institute  of  Heat  and  Mass  Transfer  of  the  Academy  of 
Sciences  Belorussian  SSR.  The  electroconductivity  *  was  determined  for 
direct  current. 

The  dielectric  permeability  e*  and  the  Iobs  angle  of  slope  tan  S  were 
measured  in  the  frequency  range  5-200  KHz.  All  these  characteristics  were 
determined  at  normal  temperature.  To  investigate  the  effects  of  filling, 
relaxation  processes  in  the  polymer  were  measured  as  well  as  the  temperature 
dependence  of  e'  and  tan  5  in  the  temperature  range  20-250*0.  The  specimen 
temperature  was  measured  by  a  chromel-copel  thermocouple  with  an  error  not 
greater  than  *lcC.  The  error  in  determining  e’  and  tan  5  did  not  exceed 
±3  percent,  and  for  X  and  x  —  *12  percent. 
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Figure  1  presents  the  curves  of  heat  and  electrical  conductivity  of 
coatings,  and  also  their  dielectric  characteristics  as  functions  of  filler 
concentration. 

An  increase  in  the  first  two  indices  with  rise  in  TiOg  content  is 
accounted  for  by  the  fact  that  A  and  x  of  the  filler  are  approximately  two 
orders  of  magnitude  higher  than  those  for  the  polymeric  binding  agent.  In 
addition,  introduction  of  filler  evidently  leads  to  a  change  in  the  mechanism 
of  thermal  conductivity  and  charge  transfer.  Thus,  according  to  modern 
concepts,  heat  transfer  in  a  polymer  occurs  mainly  via  phonon  transport  [2], 
to  which,  with  addition  of  filler,  is  added  heat  transfer  by  weakly  bound 
electrons  of  the  oxide  crystals.  The  process  of  charge  transfer  in  a 
polymer,  in  the  view  of  many  authors  [3],  follows  an  ionic  mechanism;  in  the 
case,  however,  of  titanium  dioxide  electronic  conductivity  occurs.  Thus, 
introduction  of  filler  leads  to  complication  of  heat  and  charge  transfer 
processes. 

The  dielectric  permeability  of  the  coating  rises  with  rise  in  Ti02 
content,  but  the  loss  slope  angle  decreases,  since  for  these  indices  as  well 
the  difference  in  absolute  values  between  the  filler  and  the  binder  is 
one  to  two  order  of  magnitude. 

Along  with  experimental  investigation  of  thermophysical  and  electrical 
characteristics  of  coatings  as  functions  of  the  filler  content,  an  analytical 
calculation  was  made  of  these  same  parameters  using  the  formulas  proposed  by 
V.  I.  Odelevskiy  for  two* phase  systems: 

a)  for  metrical  distribution  with  nonelongated  inclusions: 


b)  for  statistical  distribution  of  the  disperse  phase: 


wnere  ana 
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phase.  It  has  oeen  ascertained  that  the  beat  approximation  to  experimental 
data  is  provided  by  the  relationship  (1),  cf.  the  table. 


Experimental  and  Calculated  Values  of  e'  of  Coatings 


Value  of  | 

Bulk  Content  of  Disperse  Phase. 

q  r 

9 

• 

10 

;  »  * 

30  | 

40  {  50 

Experimental 

Data 

4,2 

5,3 

6,2 

7,9  9,1 

Calculated  Data 
(Matrical  Dis¬ 
tribution 

4,12 

5, IB 

5,66 

* 

6.9  10.7 

Calculation 

(Statistical . 

4,3 

6,5 

n.o 

19,0  31,0 

Distribution) 
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It  was  shown  earlier  that  the  change  of  effective  thermal  conductivity 
of  linear  polymer  (acrylic  copolymer  AC)-based  coatings  with  rutile  filling 
corresponds  to  the  metrical  distribution  [4].  ThiB  has  also  been  confirmed 
for  the  polymer  of  reticulate  structure  (Figure  1). 

In  the  calculation  filler  characteristics  taken  from  literature  data 
were  employed.  The  difference  in  the  rate  of  change  of  effective  electro¬ 
conductivity  —  both  experimental  and  calculated  --  is  evidently  due  to 
the  low  degree  of  precision  in  determining  the  x  of  rutile  [5]. 

The  temperature  dependence  of  c'  and  tan  6  of  coatings  at  a  frequency  of 
200  KHz  is  shown  in  Figure  2.  The  nature  of  the  variation  of  dielectric 
permeability  with  temperature  is  identical  for  coatings  with  different  Ti02 
content.  The  course  of  curves  describing  the  temperature  dependence  of  tan  6 
is  precisely  analogous  both  for  unfilled  and  for  filled  polymers:  dielec¬ 
trical  losses  rise,  reaching  a  relaxation  dipole-segmental  maximum. 


A  fut liter  increase  in  losses  is  caused  by  consumption  of  the  electric  field 
energy  for  continuous  electroconductlvlty.  These  functions  were  derived  when 
e'  and  tan  6  ware  measured  at  other  frequencies.  A  decrease  In  the  absolute 
value  of  tan  6  both  in  the  maximum  region  and  in  the  region  of  conductivity 
losses  is  observed  with  rise  In  filler  concentration.  Additionally,  it  was 
established  that  Introducing  small  amounts  of  TIOg  (10  percent  by  volume) 
into  the  binder  produces  a  shift  in  the  relaxation  loss  maximum  toward  the 
side  of  lower  temperatures.  Reduction  of  the  vacancy  rate  of  the  three- 
dimensional  network  by  changing  the  resin-solidifier  ratio  leads  to  a  similar 
result.  Evidently,  in  the  solidifying  process  some  of  the  relaxation-capable 
groups  of  resin  and  solldlflar  are  screened  by  filler  particles.  For  a 
20-percent  rutile  content  the  shift  of  the  maximum  is  less  pronounced,  A 
further  rise  in  the  filler  content  leads  to  a  marked  drop  in  the  quantity 
tan  6,  and  the  position  of  the  maximum  is  difficult  to  fix.  It  can  be 
assumed  that  for  an  appreciable  filler  content  the  reduction  in  the  number  of 
transverse  intermolecular  bonds  ("cross-linking")  is  compensated  by  the 
adsorptive  powers  of  Interaction  of  the  molecule  and  the  solid  particle.  As 
a  result,  the  segmental  mobility  does  not  change  and  the  shift  of  the  dipole- 
segmental  loss  maximum  is  not  observed. 

Symbols 

-  mean  molecular  weight  of  the  polymer;  \  ■  thermal  conductivity; 

A  ■  generalized  conductivity;  x  f  electroconductlvlty;  G  ■  bulk  fraction; 
s'  -  dielectric  permeability;  £  ■  angle  of  losses. 
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Figure  1.  Effect  of  titanium  dioxido  content  on  X,  Xi  e'  and 
tan  6  of  coatings:  1  and  l1,  Thermal  conductivity  X;  2  and 
2^,  Electroconductivity  x»  3  and  3*,  Dielectric  permeability  e * ; 

4  and  4^t  Loss  slope  tan  6  at  frequency  f  ■  3  KHz.  1,  2,  3, 

4  are  experimental  data;  l1,  21,31,  41  are  calculated  data  using 
the  relationship  (1) . 


Figure  2,  Temperature  dependence  of  e*  and  tan  {  of  coatings  at 
a  frequency  of  200  KH2:  O,  Without  filler;  □,  10  percent  Ti02; 

A,  20  percent  I102;  *•  30  percent  Ti02;  x,  40  percent  T102 
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THERKAL  CONDUCTIVITY  OF  AIR  AT  TEMPERATURES  UP  TO  800 *K  AND 
PRESSURES  UP  TO  1000  BAR 

A.  A.  Tarzimanov  and  V.  S.  Lozovoy 

In  recent  years  thenaophyslcal  studies  have  gained  Increasing  attention. 
However,  investigations  dealing  with  measurements  of  thermal  conductivity  X 
of  gases  at  P  >  500  bar  and  T  >  400®K  are  very  few  in  number  [1,  2].  We  have 
devised  a  new  experimental  apparatus  using  the  heated  wire  method  to  measure 
the  X  of  several  gases  and  liquids,  above  all,  air.  The  method  of  measure¬ 
ment  is  analogous  to  that  adopted  In  [3]. 

A  diagram  of  the  experimental  apparatus  is  shown  In  Figure  1.  The 
measuring  tube  made  of  "Pyrex"  glass  was  placed  within  a  thick-walled  auto¬ 
clave  made  of  the  stainless  steel  lKhl8N9T.  The  autoclave  was  equipped  with 
three  external  electric  heaters  supplied  with  stabilized  alternating  current. 
Control  of  temperature  field  uniformity  over  the  autoclave  height  was 
maintained  by  five  chromel-alumel  thermocouples.  The  lead-out  of  the 
measuring  wires  from  the  high-pressure  region  made  use  of  Teflon  sleeves. 

A  thermocompressor,  whose  operating  principle  is  analogous  to  that 
described  In  [A],  was  used  to  produce  high  pressure  In  the  system.  By  means 
of  the  thermocompressor  it  was  possible  to  raise  the  pressure  from  120- 
150  bar  to  approximately  500  bar.  A  further  pressure  rise  to  1000  bar  and 
regulation  in  the  course  of  the  experiments  was  achieved  with  the  aid  of 
a  high-pressure  oil  press  and  a  U-shaped  compressing  vessel  filled  with 
mercury. 

2 

Model  manometers  at  600,  1000  and  1600  kg/cm  were  used  to  measure 
pressure,  and  to  their  readings  corrections  were  added  based  on  calibration 
of  the  manometers  with  the  MP-2500  piston  manometer  of  the  0.02  clasB. 
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The  design  of  the  measuring  tube  is  similar  to  that  used  in  [3].  Its 
main  dimensions  are  as  follows:  inner  diameter  of  capillary  —  0.928  mm, 
outer  capillary  diameter  —  1.876  mm,  diameter  of  wire  of  inner  resistance 
thermometer  —  0.100  mm,  length  of  measuring  section  for  the  first  series  of 
experiments  —  88.9  mm,  and  after  their  assembly  —  89.5  mm.  The  data  of 
control  measurements  of  thermal  conductivity  of  air  at  atmospheric  pressure 
at  temperatures  up  to  800°K  gave  satisfactory  agreement  with  known  reliable 
data  of  other  researchers. 

For  the  investigations  clean,  dried  air  taken  from  a  high-pressure 
still  for  air  separation  was  used.  Two  drier-filters  filled  with  alumogel 
were  additionally  included  in  the  layout  (Figure  1).  Measurements  of  the 
X  of  air  under  pressure  are  presented  on  the  isotherms.  In  all  the 
experiments  the  value  of  the  product  of  the  criteria  Gr*Pr  did  not  exceed 
1000,  which  doubtless  points  to  the  absence  of  convective  currents  within 
the  capillary. 

In  interpreting  the  measurement  results  corrections  were  introduced  that 
allowed  for  the  removal  of  heat  from  the  heater  endB  in  the  temperature  drop 
at  the  measuring  tube  walls,  for  variation  of  geometric  dimensions  with 
temperature,  and  for  radiative  heat , transfer.  In  view  of  the  effect  that 
elongation  of  the  spring  has  on  the  inner  thermocouple  readings,  the  latter 
was  calibrated  during  the  course  of  the  experiments  against  the  outer 
resistance  thermometer.  To  verify  the  effect  that  high  pressures  and 
temperatures  have  on  platinum  resistance  thermometer  readings,  a  series  of 
special  measurements  was  made.  Well  defined  pressure  functions  of  platinum 
resistance  were  obtained  in  the  range  T  ■  300-800°K  and  p  ■  1-1000  bar.  In 
calculating  the  X  of  air,  a  correction  was  introduced  for  variation  of 
platinum  resistance  with  temperature.  The  possible  experimental  error  in  our 
estimation  did  not  exceed  1.5  percent.  Figure  2  presents  isobars  with 
experimental  pointB  given  for  these  pressures.  The  table  lists  results  of 
experiments  where  the  following  symbols  were  adopted:  p  ■  pressure, 

T  ■  mean  air  temperature,  AT  ■  temperature  difference  in  air  layer,  and 
X  ■  coefficient  of  thermal  conductivity  with  all  the  corrections  taken  into 
account. 
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It  must  be  noted  that  the  thermal  conductivity  of  air  as  a  function  of 
pressure  has  thus  far  not  been  adequately  studied ,  which  is  evidenced  by  the 
meagemess  of  experimental  data  [5(  6];  further,  systematic  discrepancy  [7] 
la  observed  for  the  data  [5,  6]. 

The  data  we  obtained  at  pressures  up  to  1000  bar  and  temperatures  up  to 
800*K  relate  mainly  to  the  region  of  parameters  of  state  uninvestigated  thus 
far.  Their  comparison  with  the  results  of  other  researchers  is  most 
conveniently  presented  in  the  form  of  excess  thermal  conductivity  A  -  XT 
(A^  is  the  thermal  conductivity  at  p  *  1  bar)  as  a  function  of  density  p. 

As  we  can  see  from  Figure  3,  our  experimental  points  and  the  data  of 
I.  F.  Golubev  fit  along  the  averaged  curve  with  relatively  Bmall  scatter, 
the  exception  here  being  the  196*K  isotherm  [6],  the  points  of  which  are 
overstated  by  up  to  8-10  percent.  The  results  in  [5]  are  understated  by 
5-7  percent. 

Close  study  of  the  data  shows  that  in  the  coordinates  (X  -  A^)  -  f(p), 

strictly  speaking,  no  single  curve  is  obtained.  The  isotherms,  though  very 

weakly,  are  stratified  (Figure  4).  At  a  density  close  to  the  critical  point 
o 

(p  -  313  kg/m  ),  stratifying  of  the  300  and  786 °K  isotherms  amounts  to 
3  percent.  A  similar  state  of  affairs  has  been  obtained  for  the  thermal 
conductivity  of  nitrogen  [1]. 
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Figure  1.  Diagram  of  experimental  apparatus:  1,  Thermocouples; 
2,  Heat  insulation;  3,  Electric  heaters;  4,  Autoclave  housing; 

5,  Measuring  tube;  6,  Cooler;  7,  Assembly  for  leading  out  of 
measuring  wires;  8,  Magnet;  9,  Compressing  vessel;  10,  High- 
pressure  oil  press;  11,  Piston  manometer;  12,  Gas  cylinder; 

13,  Filter-driers;  14,  Thermocompressor  cylinders;  15,  Valves; 
16,  Electric  heater;  17,  Lens  gaskets 


Figure  2.  Thermal  conductivity  of  air  as  a  function  of  temper¬ 
ature  (iBobars) 


Figure  3.  Excess  thermal  conductivity  of  air  as  a  function  of 
density:  1,  [5];  2,  [6];  3,  Authors'  data 


Figure  4.  Excess  thermal  conductivity  of  air  as  a  function  of 
density  for  the  300°  and  786#K  isotherms 
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THERMAL  CONDUCTIVITY  OF  THIN  BOUNDARY  FILMS  OF  LIQUIDS  STABILIZED  BY 

THE  SURFACE  OF  A  SOLID 

G.  T.  Timoshenko  and  M.  S.  Mats lk 

Preliminary  experlmenta  by  one  of  the  authors  established  the  fact  of 
the  anomalously  high  thermal  conductivity  of  films  of  water  on  the  surface  of 
mica  crystals  [3].  Ws  made  a  further  careful  examination  of  this  effect. 

A  method  of  determining  thermal  conductivity  of  thin  liquid  fllmB  has 
been  formulated  [1-4]  baaed  on  the  general  regular  regime  theory. 

Saturating  a  pile  of  a  large  number  of  identically  thick  sheets  of  any 
material  with  the  teat 'liquid  and  compressing  the  pile  in  a  vise,  we  can  get 
the  corresponding  thickness  of  the  liquid  films..  The  thermal  conductivity  of 
the  film  liquid  is  calculated  from  the  formula: 


Here ; 

X,  Y  And  Z  are  the  dimensions  of  the  pile,  measured  in  a  direction 
perpendicular  to  the  sheet  plane; ] 

X"  and  X1  are  the  components  of  the  thermal  conductivity  of  plates  in 
the  longitudinal  and  normal  directions,  respectively  (given  the  presence  of 
thermal  conductivity  anisotropy); 

c  and  p  are  the  specific  heat  and  density  of  the  pile  with  absorbed 
liquid  taken  into  account; 

mj  and  m2  are  cooling  rates  of  the  pile  in  the  regular  regime  for  a  heat 
transfer  coefficient  a  ■*  «  for  thf  corresponding  Z ^  and  Z2  given  by  compres¬ 
sion. 


Knowing  Che  cooling  rate  and  the  height  2^  of  Che  pile  tC  maximum 
compression  (in  Che  limiting  case,  a  continuous  specimen,  for  example  a 
monocrystal  of  mica) ,  we  can  determine  Che  mean  value  of  thermal  conductivity 
of  liquid  films  of  any  thickness,  by  setting  Che  corresponding  pile  height  Z 2 
and  experimentally  determining  the  value  of  Che  cooling  rate  m2*  Since 
according  to  the  conception  of  Che  method  the  thermal  conductivity  of  a  pile 
Zj  in  height  saturated  with  liquid  must  be  compared  with  the  thermal 
conductivity  of  a  continuous  specimen  of  the  same  dlmenalon,  a  correction 
allowing  for  reduction  of  the  cooling  rate  of  the  continuous  specimen  when 
its  height  is  Increased  to  Z2  is  introduced  into  the  measured  value  of  the 
cooling  rate  m^.  Studies  were  made  of  the  thermal  conductivity  of  film 
liquids  on  mica.  The  following  method  was  employed. 

The  working  arm  of  a  differential  copper-conatantan  thermocouple 

made  of  0.1-mn-diameter  wires  without  reinforcement  was  placed  in  the  center 

of  a  previously  cleaved  mica  crystal  (muscovite  or  phlogopite)  cut  in  the 

2 

form  of  a  rectangular  parallelepiped  with  e  base  4  *  5  cm  and  about  4.5  cm 
in  height. 

Then  the  crystal  was  compressed  in  a  vise  and  according  to  the  regular 
regime  method  [5]  Its  cooling  rate  was  determined  for  5s  in  the  thermostat 
whose  temperature  was  kept  unchanged  to  a  precision  of  0.02  percent, 
with  the  specimen  heat  transfer  coefficient  a  •*  <*>,  A  mirror  galvanometer 
with  a  sensitivity  of  0.5*10”*’  v/mm  was  employed.  The  time  intervals  were 
fixed  with  a  precision  of  0.2  sec.  The  cooling  rate  was  determined  with  an 
error  of  0.002-0.005  sec. 

Later  the  mica  crystal  was  converted  into  a  pile  of  sheets.  For  this 
purpose  a  muscovite  crystal  was  accurately  cut  into  30-40  layers  1.0-1. 5  mm 
thick,  without  varying  the  mutual  positions  of  the  sheets  in  tho  pile,  the 
sheets  then,  in  turn,  were  carefully  stacked  to  a  height  of  20-30  sheets. 
However,  small  (3-5  mm  )  contact  sites  remained  along  the  margins  between 
the  sheets,  making  it  impossible  to  shift  them  relative  to  each  other.  With 
careful  shaping  and  using  the  appropriate  compressive  force,  this  pile  became 
the  original  crystal  without  visible  change  in  the  appearance  of  the  side 
surfaces. 
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To  stratify  a  phlogopite  crystal  the  lov-temperature  swelling  effect  was 
uaed  (6,  7].  A  swollen  phlogopite  crystal  Is  an  aggregate  of  a  large  number 
of  plates  fastened  to  each  other  along  the  margin  on  any  one  side.  Upon 
being  cooled,  this  crystal  is  converted  into  its  initial  state  without  any 
change  in  the  end  surfaces. 

Saturation  of  the  pile  with  the  test  liquid  was  carried  out  in  vacuum. 
The  cooling  rate  of  the  pile  for  a  different  bulk  concentration  of  the 
Impregnating  liquids  was  determined  un^er  the  same  set  of  conditions  as  for 
the  Integral  crystal.  The  and  Z2  values  were  determined  with  a  micrometer 
with  a  precision  of  up  to  5  p .  The  values  of  the  thermal  characteristics  of 
mica  were  taken  from  [2].  Interpretation  of  the  results  of  measuring 
thermal  conductivity  of  liquid  films  followed  the  Kornfel’d  method  [8]  with 
a  reliability  of  90  percent  out  of  not  less  than  30  experiments  for  each 
film  thickness.  The  relative  error  In  thermal  conductivity  measurements  was 
10-15  percent. 

Distilled  water,  ethyl  alcohol,  carbon  tetrachloride  and  transformer 
oil  were  Investigated.  To  reduce  the  probability  of  the  teat  liquid  being 
washed  out  of  the  pile  with  the  thermostat  water,  the  experiments  with 
alcohol,  carbon  tetrachloride  and  transformer  oil  were  carried  out  for 
relatively  small  film  thicknesses. 

From  the  tabulated  data  it  follows  that  the  thermal  conductivity  of 
these  nonpolar  liquids,  such  as  carbon  tetrachloride  and  transformer  oil, 
within  limits  of  experimental  error  do  not  change  with  variation  of  film 
thickness  and,  consequently,  no  anomalous  thermal  conductivity  was  found  for 
them.  The  thermal  conductivity  values  found  from  formula  (1)  for  these 
liquids  is  approximately  one  order  of  magnitude  different  from  the  values 
presented  in  handbook  literature. 

A  different  state  of  affairs  is  observed  for  films  of  polar  liquids.  As 
we  see  from  the  table,  the  effect  of  anomaloua  thermal  conductivity  is  quite 
distinctly  evidenced  for  ethyl  alcohol  and  especially  for  water.  In 
Figure  1,  curve  1  describes  the  variation  of  thermal  conductivity  of  water 
films  with  film  thickness  on  fresh  surfaces  of  mica  sheets.  The  thermal 


0 
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conductivity  of  water  rises  raoidly  with  decrease  in  film  thickness  and 
reaches  a  value  of  60-65  w/m*dag  for  thicknesses  less  than  0.1  u,  exceeding 
the  thermal  conductivity  of  bulk  water  (0.6  w/m*deg)  by  two  orders  of 
magnitude  and  is  close  to  the  thermal  conductivity  of  metals. 

The  thermal  conductivity  of  thin  layers  of  ethyl  alcohol  on  mica 
surfaces  is  also  ten  times  higher  than  its  thermal  conductivity  In  the  bulk 
(0.16  w/m*deg)  [9].  The  possibility  of  accounting  for  the  anomalous  thermal 
conductivity  effect  seen  in  alcohol  films  by  its  partial  replacement  with 
thermostat  water  was  precluded  by  the  results  of  experiments  with  carbon 
tetrachloride. 

The  magnitude  of  the  effect  cannot  be  accounted  for  by  possible 
variation  of  side  surface  roughness  in  the  transition  from  crystal  to  the 
pile.  The  very  method  of  obtaining  the  pile  from  the  crystal  provides  for 
the  invariancy  of  the  mutual  arrangement  of  the  individual  parts  of  the 
specimen.  Additionally,  in  previously  conducted  studies  [3]  the  effect  of 
anomalous  thermal  conductivity  of  water  film  was  clearly  observed  when 
the  pile  was  made  of  individual  identically  fixed  sheets  of  mica  and  when  a 
pile  with  greatest  compression  was  taken  as  the  original  state.  Under 
these  conditions  it  is  obvious  that  the  roughness  of  the  end  surfaces 
remained  unchanged.  Data  given  in  Figure  2  also  evidences  the  small  role  of 
roughness.  These  data  make  clear  that  the  cooling  rate  of  a  water-saturated 
pile  passes  through  a  maximum  with  increasing  film  thickness  and  for  thick 
enough  films  (maximum  roughness)  becomes  close  to  the  cooling  rate  of  the 
uncleaved  crystal. 

The  effect  of  anomalous  thermal  conductivity  of  thin  films  of  polar 
liquids  depends  substantially  on  the  state  of  the  crystal  surface.  The 
thermal  conductivity  of  water  films  in  a  pile  was  investigated  after  the  mica 
sheets  were  rinsed  in  boiling  water.  This  treatment  leachea  out  the  surface 
and  thus  weakens  the  surface  field,  which  leads  to  a  considerable  decrease 
of  thermal  conductivity  of  water  films  (Figure  1,  curve  2). 

When  a  pile  made  of  polyethylene  sheets  is  saturated  with  water, 
practically  no  anomaly  in  thermal  conductivity  is  observed  (curve  3). 


It  is  evident  thet  the  effect  of  anomalous  thermal  conductivity  of  thin 
films  rises  with  decrease  in  temperature,  and  it  can  be  assumed  that  it  must 
be  most  clearly  evident  that  temperatures  close  to  0°C.  The  thermal 
conductivity  of  water  films  drops  off  rapidly,  however,  with  rise  in 
temperature,  and  at  60-65BC.  becomes  equal  to  the  thermal  conductivity  of  bulk 
water. 

Results  of  the  experiments  described  confirm  the  concept  of  the  quasi¬ 
crystalline  structure  of  thin  films  of  polar  liquids  that  are  in  interface 
with  the  surface  of  a  solid  [lti-12]. 

As  we  know,  in  the  crystalline  state  compounds  are  more  heat-conductive 
compared  to  the  amorphous  state.  Crystalline  quartz,  for  example,  has  a 
thermal  conductivity  that  is  an  order  of  magnitude  higher  than  the  thermal 
conductivity  of  fused  quartz  [13].  It  is  obvious  that  structure-ordered 
films  of  polar  liquids,  and  above  all  water,  also  exhibit  an  exceptional 
capacity  to  transmit  thermal  agitations  along  films.  However,  with  rise  in 
temperature  their  structure  begins  to  break  down,  the  water  films  apparently 
become  amorphlzed.  At  a  temperature  of  GS*C  total  amorphization  sets 
in,  to  which  corresponds  the  rapid  drop  in  thermal  conductivity  down  to  the 
values  characteristic  of  bulk  water. 

A  similar  temperature  dependence  of  thermal  conductivity  has  been 
observed  for  crystalline  compounds  in  thelri  transition  from  'he  solid  state 
to  the  liquid  [13]. 

Based  on  the  studies  made,  the  following  conclut:  .v?  <  m  be  drawn: 

1.  Experiments  with  cleaved  mica  crystals  graphic]  ly  demonstrate  the 
existence  of  the  anomalous  thermal  conductivity  effect  Cor  films  of  polar 
liquids  on  mica  surfaces. 

2.  The  magnitude  of  the  effect  depends  on  the  nature  of  the  liquids  and 
rises  with  decrease  in  film  thickness;  water  films  about  0.1  u  in  thickness 
on  mica  surfaces  have  a  thermal  conductivity  that  la  one  order  of  magnitude 
less  than  that  of  metal. 
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Thermal  Conductivity  of  Film  Liquids  on  Mica  Crystals 


Cooling 
Rate  of 

.Cooling  ,  .Pile  In- 


.cooling  ,  rixe  .  /Mean  Thsr- 
Bulk  Thickness  *?te  of  SI!?’**1®'5  «•!  Ccnduc 

C.  .  Concentre- of  Liquid  JlSuid  *  tlvlty  of 

Llquld  tion  of  Layers,  £ry«al,  Liquid 

Liquid  in  2h,  u  “l  10  »  m2  10  '  Layer,  X, 

P 


2,20 

0, 

Distilled 

jjf 

0, 

Water 

S' 

0. 

oils 

0 

oils 

o! 

Ethyl  : 
Alcohol 


Distilled  „  g.|9 
Water  •  q!36 

Oil? 


0,70 

Transformer  0,40 
Oil  °'20 


IE  8jS  I 


90  0,42  .  I 

s$  n  i 

17  0,08  I 


Distilled  : 

?: 

gg  :  i 

Water 

0 

50  -  I 

0 

30  -  I 

& 

60 

0 

15  i! 

o 

1-1 

•d 

A- 

Transformer 

Oil 

11  :  i: 

30  -  I, 

# 


88  i 

IS?  i 

? 

610 

I 

395 

400 

I 

•  I 

405 

I 

405 

I 

410 

I 

-623- 


No  anomalv  In  t-hamal  nnmi.<.k4.^>-  j  «  .  . 

'  - - — x.  uuuarvea  lor  nonpolar  liquids. 

3.  The  size  of  the  effect  depends  on  the  state  of  the  surface  and  the 
nature. of  the  solid.  The  effect  is  most  graphically  evident  on  "fresh" 

surfaces  of  mica  crystals.  Leaching  out  the  surface  leads  to  a  sharp  drop  in 
thermal  conductivity  of  film  water. 

4.  The  anomalously  high  thermal  conductivity  of  films  of  polar  liquids 
is  accounted  for  by  the  improved  quasi-crystalline  structure  acquired  as  a 
result  of  the  action  of  these  liquids  on  molecules  of  a  strong  surface  layer 
of  mica.  With  temperature  rise,  amorphization  of  film  structure  sets  in  and 
film  thermal  conductivity  drops  off  rapidly  down  to  values  characteristic  of 
the  bulk  state. 

Symbols 

X  -  coefficient  of  thermal  conductivity;  c  -  specific  heat;  p  «  density 
(bulk  weight);  m  ■  cooling  rate  of  specimen  in  regular  regime;  a  -  heat 
transfer  coefficient  at  specimen  surface;  X,  Y  and  Z  -  linear  dimensions  of 
specimen;  2h  -  thickness  of  liquid  film  between  two  flat  surfaces. 
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Figure  1.  Variation  of  thermal  conductivity  of  water  filmn  with 
film  thickness:  1,  On  "fresh"  surfaces  of  mice  crystals; 

2,  On  surfaces  of  mica  crystals  treated  In  boiling  water; 

3,  On  polyethylene  surfaces 


Figure  2,  Cooling  rate  of  specimens  as  a  function  of  their 
height:  1,  For  continuous  mica  crystal  at  constant  density; 
2,  For  cleaved  crystal  saturated  with  water;  3,  For  cleaved 
crystal  saturated  with  carbon  tetrachloride 


toe 


Figure  3.  Thermal  conductivity  of  film  water 
as  a  function  of  temperature 
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THERMAL  PROPERTIES  OF  SOLID  AND  MOLTEN  METALS  AT  KT3H  TEMPERATURES 

L.  P.  Filippov,  A.  V.  Arutyunov,  I.  N.  Makarenko, 

I.  P.  Mardykln,  L.  N.  Trukhanova, 

B.  N.  Khusainova  and  R.  F.  Yurchak 

The  article  is  a  brief  description  of  studies  aimed  at  experimental 
investigation  of  the  basic  regularities  underlying  the  behavior  of  thermal 
conductivity,  thermal  dlffusivity,  heat  capacity  and  electroconductivity 
of  solid  and  molten  metals  and  in  part  of  semiconductors  at  temperatures 
from  1000°  to  3000°K. 

A  necessary  stage  in  carrying  out  the  projected  research  program 
involved  an  effort  aimed  at  improving  existing  methods  of  measurement  and 
devising  new,  high-temperature  methods.  The  main  focus  here  was  given  to 
developing  comprehensive  methods  making  It  possible  to  obtain  simultaneously 
the  entire  totality  of  main  thermal  characteristics.  At  present  we  have 
four  methods  of  this  type  available  to  us.  These  methods  have  been  based  on 
the  use  of  the.  so-called  third-orde,.  regular  .thermal' regime  — •  a  steady-state 
periodic  process.  The  main  merit  of  this  process  is  the  diversity  of 
information  obtained  in  the  experiment.  (The  following  information 
sources  were  used:  field  of  mean  temperatures,  field  of  amplitudes  and 
phases  of  several  harmonic  components  of  temperature  fluctuations  [1,  2].) 
Below  follows  a  brief  description  of  these  methods. 

The  work  aimed  at  devising  a  method  for  simultaneous  determination  of 
thermal  dlffusivity  and  thermal  conductivity  of  methods  based  on  use  of 
temperature  fields  propagated  along  the  axis  of  a  cylindrical  specimen  has 
been  published  earlier  [3].  Usually  this  kind  of  process  is  employed  In 
measuring  only  thermal  dlffusivity  (the  Angstrom  method).  The  possibility  of 
determining  also  thermal  conductivity  here  Involves  use  of  variable  heating 
by  electronic  bombardment  making  it  possible  to  determine  quite  precisely  the 


value  of  the  heating  Intensity.  The  method  allows  us  to  measure  the 
parameters  indicated  for  the  "long"  (several  centimeters)  and  "short"  (about 
one  centimeter)  specimens  with  a  precision  of  6-8  percent  at  temperatures  up 
to  approximately  1500®.  Of  late,  a  comparative  variant  of  the  method  has 
also  been  devised  making  it  possible  to  bypass  electronic  heating.  Study  is 
under  way  on  using  the  method  for  molten  metals  and  semiconductors. 

The  second  of  the  comprehensive  methods  developed  is  based  on  use  of 
temperature  fluctuations  of  cylindrical  symmetry  —  radial  temperature  waves. 
We  earlier  [3]  reported  several  versions  of  the  use  of  radial  temperature 
waves  in  measuring  thermal  diffusivity,  Use  of  electronic  heating  in  this 
method  allows  us  to  determine,  along  with  thermal  diffusivity,  also  heat 
capacity  and,  as  a  consequence,  thermal  conductivity.  We  used  a  modification 
of  the  method  in  which  heating  was  conducted  within  a  hollow  specimen  for 
these  comprehensive  measurements.  For  this  purpose,  the  incandescent  cathode 
imitating  electrons  is  placed  along  the  specimen  axis.  A  measured  voltage  of 
several  hundred  volts  is  applied  between  the  cathode  and  the  specimen.  The 
temperature  fluctuations  at  the  outer  specimen  surface  are  recorded. 

Thermal  diffusivity  is  determined  from  the  phaBe  difference  between  heating 
intensity  fluctuations  and  temperature  fluctuations;  heat  capacity  1b  found 
from  the  value  of  the  amplitude  of  temperature  and  heating  intensity  fluctu¬ 
ations  [4],  The  maximum  error  in  determining  thermal  diffusivity  is  usually 
4-6  percent,  while  the  maximum  error  in  heat  capacity  measurement  varies  from 
3  to  6  percent.  A  contactless  photoelectric  method  is  used  in  the  latter, 
high-temperature  version  of  the1  apparatus  for  recording  temperature 
fulctuations;  maximum  temperatures  here  are  approximately  1800°K. 

The  method  described  was  used  by  the  authors  mainly  in  Investigating 
molten  metals.  Liquid  metal  was  placed  in  a  vessel  made  of  two  thin-walled 
(approximately  0.1  mm)  tantalum  tubes  coaxial  with  aach  other.  Thermal 
properties  of  tin,  lead,  bismuth,  cadmium,  antimony,  copper  and  germanium 
(5,  6]  were  studied  by  this  method. 

One  of  the  main  results  of  this  investigation  was  a  conclusion  of  the 
weak  dependence  of  thermal  conductivity  of  the  materials  ptudied  on  temper¬ 
ature.  Thus,  for  tin  thermal  conductivity  of  a  temperature  range  covering 
1100  degrees  centigrade  varied  by  not  more  than  10  percent.  The  behavior  of 
the  Lorentz  number  of  liquid  metals  is  remarkable.  At  temperatures  close 


Co  the  melting  point  the  Lorentz  number  L,  as  a  rule,  is  dose  to  the  theor- 

"8  2 

etlcal  value  of  2.45*10  w ohms /deg .  As  the  temperature  rises,  however,  a 

marked  progressive  drop  In  L  is  observed,  which  provides  the  necessary 
formulation  of  the  problem  of  what  causes  underlie  this  regularity. 

The  third  of  the  comprehensive  methods  we  developed  was  also  based  on 
use  of  radial  temperature  waves.  In  contrast  to  the  former  method,  these 
waves  are  Induced  by  heating  a  cylindrical  specimen  in  an  induction  furnace, 
the  power  of  which  is  varied  periodically  (modulated) .  Owing  to  the 
existence  of  a  skin-effect,  liberation  of  heat  in  the  specimen  occurs  In  a 
relatively  thin  layer  near  the  surface;  accounting  for  the  bulk  nature  of 
heat  release  1b  nonetheless  still  necessary.  In  the  version  of  the  method 
now  used,  the  temperature  fluctuations  of  the  outer,  heated  surface  are 
recorded,  though  it  is  possible  that  there  is  another  modification  of  the 
method  —  with  temperature  fluctuations  recorded  at  the  specimen  axis.  A 
contactless  photoelectric  method  of  recording  temperature  fluctuations  is 
employed.  As  in  the  previous  method,  in  determining  the  thermel  dlffuslvity 
of  the  specimen  it  Is  enough  to  know  the  phase  difference  between  power 
fluctuations  and  temperature  fluctuation.  In  determining  the  power  a  coll 
made  up  of  several  turns  of  wire  surrounding  the  specimen  is  used.  The  value 
of  the  electromotive  force  recorded  from  this  coil  can  be  associated  with  the 
field  Intensity  within  the  inductor  and,  ultimately,  with  the  heater  power. 
The  measurement  method  described  allows  us  to  determine  the  thermal  diffusiv- 
lty,  thermal  conductivity  and  heat  capacity  in  the  temperature  range 
1000-2500°K  with  errors  of  about  6-8  percent  and  8-10  percent,  respectively. 
The  method  is  used  for  systematic  measurements  of  thermal  properties  of  solid 
refractory  metals.  t 

The  fourth  type  of  the  comprehensive  experiment  was  designed  for 
specimens  in  wire  or  foil  form.  It  combines  measurements  of  heat  capacity  by 
the  alternating  current  heating  method  with  thermal  conductivity  measured  by 
the  steady-state  method.  Both  methods  separately  have  been  described  by  us 
[3]  (cf.  also  [7-9]).  Lately,  these  methods  have  been  combined  so  that  both 
heat  capacity  and  thermal  conductivity  are  determined  on  the  same  specimen  by 
means  of  a  single  common  apparatus,  the  basis  of  which  is  the  earlier 


described  differential  optical  pyrometer  [10,  11].  The  electroconductivity 
of  the  specimen  and  its  radiative  characteristics  are  also  described  with  the 
same  device  (the  absolute  black  body  models  are  used  for  this  purpose).  The 
error  of  heat  capacity  determination  is  4-6  percent,  and  6-8  percent  for 
thermal  conductivity;  the  temperature  range  is  1000-3000*K. 

The  comprehensive  methods  for  measuring  thermal  characteristics 
described  do  not  exhaust  the  arsenal  of  means  of  Investigation  used  in  our 
work.  Thus,  in  particular,  we  used  the  ordinary  Angstrom  method  for  molten 
metals.  Recently  we  have  developed  a  steady-state  method  of  measuring 
thermal  conductivity  for  molten  metals  —  the  small  bridge  method,  earlier 
used  by  Katler  and  Cheney  [6,  12,  13]. 

Let  us  proceed  to  presenting  several  new  specific  results  obtained  only 
most  recently.  Figures  1  and  2  present  the  results  of  the  measurement  of 
thermal  conductivity  and  heat  capacity  of  niobium  obtained  by  different 
methods.  The  filled-ln  symbols  in  these  figures  refer  to  measurements  of  a 
massive  specimen  (rod)  of  niobium  of  the  following  composition:  99.2  percent 
Nb,  0.3  percent  Ta,  0.08  percent  Ti,  0.04  percent  Fe  and  0.04  percent  Si 
(specimen  1) .  The  circles  denote  values  obtained  by  the  variable  induction 
heating  method,  and  the  triangles  —  values  obtained  by  the  electronic 
heating  method,  The  unfilled-in  circles  designate  the  wire  specimen1.  Its 
composition  was  as  follows:  98,4  percent  Nb,  0.6  percent  Ta,  0.86  percent  W 
and  0.015  percent  Mo  (specimen  2). 

All  the  data  obtained  for  thermal  conductivity  agree  quite  well  with 
each  other,  and  also  with  results  of  measurements  in  [14],  Agreement  on  the 
results  of  our  heat  capacity  measurements  are  also  quite  satisfactory.  The 
data  obtained  differed  from  the  results  of  measurements  in  [15]  by  4-7  per¬ 
cent,  from  the  data  in  [16]  by  -1.5— H, 5  percent,  from  the  data  in  [17]  by 
-4-0  percent,  and  from  the  data  in  [18]  by  -3-4*1  percent.  The  results 
in  [18]  are  6-7  percent  below  ours. 

Figures  3  and  4  present  the  results  of  thermal  conductivity  and  heat 
capacity  measurements  for  antimony  (Sb  96.8  percent,  Pb  2.63  percent,  Ag 
0.50  percent,  Cu  0.05  percent  and  Fe  0.02  percent),  and  for  germanium 
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(specific  resistance  at  room  temperature  approximately  10  ohms* cm)  In  the 
solid  and  liquid  states  obtained  by  the  method  of  radial  temperature  waves 
(electronic  heating).  In  addition  to  these  data,  here  also  are  presented 
results  obtained  for  solid  germanium  (specific  resistance  at  room  temperature 
0.005  ohm*cm,  alloyed  with  Sb  impurity)  by  the  first  of  the  above-described 
comprehensive  methods2,  and  the  data  for  liquid  gallium  (99.9997  percent). 

The  values  found  for  the  thermal  conductivity  of  gallium  were  secured  by  the 
low  bridge  method.  The  data  are  in  satisfactory  agreement  with  the  results 
of  measurement  in  [20]  obtained  up  to  200SC.  The  thermal  conductivity 
results  for  solid  alloyed  germanium  agree  well  with  the  data  in  [21]. 

Small  changes  in  heat  capacity  when  antimony  and  germanium  are  fused 
agree  with  what  has  been  obtained  for  other, elements.  The  values  of  molar 
heat  capacity  in  the  molten  state,  as  for  most  liquid  metals,  are  close  to 
7-8  cal/mole.  When  we  analyzed  the  thermal  conductivity  results,  our  atten¬ 
tion  was  drawn  to  the  discontinuous  increase  in  germanium  thermal  conductivity 
upon  fusion.  The  order  of  magnitude  of  thermal  conductivity  of  liquid 
germanium  is  the  same  as  for  most  liquid  metals.  This  thermal  conductivity 
behavior  confirms  the  finding  of  the  metallic  nature  of  conductivity  in 
liquid  germanium,  which  was  arrived  at  based  on  analysis  of  electroconduct¬ 
ivity  and  the  Hall  effect  [22].  The  behavior  of  the  Lorantz  number  f0r 
liquid  antimony  and  germanium  confirmed  the  regularity  found  for  other  liquid 
metals  (cf.  above). 

Close  to  the  melting  point  the  Lorentz  number  is  close  to  the  theoretical 
value,  but  is  gradually  reduced  with  rise  in  temperature.  To  illustrate 
this.  Figure  4  presents  the  results  of  thermal  conductivity  calculations, 
indicated  by  a  dashed  line,  carried  out  on  the  assumption  of  the  validity  of 
the  Widemann-Franz  law. 


zSee  p.  634 
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Figure  1.  Results  of  investigating  the  thermal  conductivity  cf 
niobium 


liquid  germanium,  liquid  antimony  and  liquid  gallium 


Figure  4.  Heat  capacity  of  germanium  and  antimony 
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Footnote* 


1.  To  p.631  These  data  have  been  obtained  with  participation  of 
A.  G.  Sorokin. 

2.  To  p.632  Yu.  L.  Stakheyev  participated  in  the  study. 
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THERMAL  CONDUCTIVITY  OF  STABILIZED  ZIRCONIUM  DIOXIDE  AT  HIGH  TEMPERATURES 
V.  Ya.  Cheykovskiy  and  A.  M.  Banayev 

Investigations  of  the  thermal  conductivity  of  refractory  oxides  in  the 
temperature  range  1400-2000*C  thus  far  are  very  few  in  number  [1,  2],  in 
spite  of  the  fact  that  data  at  these  temperatures  is  of  Bpeclal  practical 
interest.  The  authors  of  this  present  study  have  devised  an  experimental 
apparatus  for  determining  the  thermal  conductivity  coefficient  at  temper¬ 
atures  up  to  2100° C.  The  method  of  the  experiment  and  the  apparatus  design 
have  been  set  forth  quite  adequately  in  [3].  The  measurements  were  conducted 
by  the  tubing  method. 

Specimens  of  the  test  material,  cylindrical  in  configuration,  were  in 
the  form  of  a  set  of  disks  20  mm  high,  with  outer  diameter  of  50  mm  and  inner 
diameter  of  4-6  mm.  Dead-end  vertical  openings  2  mm  in  diameter  were  made 
with  ultrasonic  equipment  to  measure  temperature  in  the  specimens.  Tungsten- 
rhenium  (VR-5/20,  0.35  mm  in  diameter)  and  platinum-platinum-rhodium 
thermocouples  in  two-channel  alundum  tubings  were  placed  into  these  aper¬ 
tures. 

The  apparatus  was  supplied  from  a  220-v  alternating  current  network 
through  the  AOSK-25  single-phase  automatic  transformer  and  the  OSU-40  step- 
down  transformer.  Determination  of  the  thermal  flux  value  was  made 
electrically,  that  is,  by  measuring  the  current  passing  through  the  heater 
and  the  voltage  drop  over  its  working  section.  The  current  value  and  the 
voltage  drop  over  the  experimental  section  of  the  heater  were  measured  by  the 
R-56  type  alternating  current  potentiometer.  The  emf  of  the  thermocouples 
was  measured  by  the  PMC-48  type  direct  current  potentiometer.  The  VR-5/20 
type  thermocouples  employed  were  calibrated  with  the  OP-48  optical  pyrometer 
in  the  modernized  TW-2  high-temperature  vacuum  furnace  [4]. 
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Characteristics  of  Zr02  Specimens  and  Their  Chemical  Analysis 


j4*4*105  21,51  60V  Alka11  °»8  0.5  I.I8  0,275  5,25  0,122  0,18 

5Q*  0,068,  MS-  Honoclinic 


The  experimental  values  of  thermal  conductivity  obtained  relate  to  the 
mean  arithmetic  temperature. 

In  this  study  thermal  conductivity  of  refractories  based  on  zirconium 
dioxide  in  the  temperature  range  300-2000°  was  measured.  Table  1  lists  the 
characteristics  of  the  specimens  investigated  and  their  chemical  composition. 
As  we  can  see  from  Table  1,  the  test  materials  differed  by  different  ratio  of 
cubic  and  monoclinic  modifications  of  zirconium  dioxide,  the  amount  of 
impurities  in  the  form  of  various  oxides,  and  the  molding  pressure. 

Crystalline  structure  had  the  greatest  effect  on  the  results  of  the  study. 

We  first  of  all  must  take  note  of  the  variation  in  the  course  of  thermal 
conductivity  as  a  function  of  temperature  for  the  polymorphous  transformation 
(lOOO-^SO’C)  of  the  monoclinic  modification  of  zirconium  dioxide  into  the 
tetragonal  modification  [5].  Up  to  the  polymorphous  transformation,  the 
thermal  conductivity  of  refractories  rises  linearly  with  temperature  rise, 
but  thereafter  it  varies  so  slightly  that  to  the  first  approximation  it 
remains  constant  within  limits  of  experimental  error.  Further,  up  to  the 
phase  transition  point  the  thermal  conductivity  depends  on  phase  corap<  sitlon: 
the  higher  the  content  of  the  cubic  and  the  lower  the  content  of  the  mono- 
clinic  modification  of  zirconium  dioxide,  the  higher  the  thermal  conductivity 
of  the  material.  We  must  also  take  note  of  the  fact  that  at  temperatures 
above  the  polymorphous  transition  experimental  data  for  refractories  1  and  2 
are  practically  the  same. 

The  molding  pressure  did  not  affect  the  thermal  conductivity  value. 

Results  of  experimental  determination  of  thermal  conductivity  of  the 
specimens  investigated  are  shown  In  Figure  1,  and  Table  2  lists  the  smoothed 
values  of  the  thermal  conductivity  coefficient  for  rounded  temperature  values. 

The  maximum  scatter  of  experimental  points  relative  to  the  averaged 
value  does  not  exceed  the  limiting  random  error  of  the  experiment,  which  the 
authors  estimate  at  12  percent. 
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Table  2.  Smoothed-Over  Values  of  Thermal  Conductivity  Coefficients 

for  Zirconium  Dioxide 


w/nrdeg _ 

Material  No.  1  i  Material  Ho.  2 
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1.48 
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Figure  1.  Results  of  experimental  study  of  thermal  conductivity 
of  zirconium  dioxide  stabilized  with  calcium  oxide:  1,  2,  Materials 
No.  1  and  2  (cf.  Table  1);  the  averaged  valuee  of  the  experimental 
data  obtained  are  presented  on  the  curves. 
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THERMAL  CONDUCTIVITY  OF  POLYCAPROAMIDE  AS  A  FUNCTION  OF  TEMPERATURE 

AND  MOLECULAR  WEIGHT 


I.  I.  Chernobyl ’skiy,  V.  I.  Kirichenko  and  A.  N.  Piven 

Modern  views  of  the  nature  of  polymers  require  some  refinement  of  the 
heat  transfer  process  In  polymers  and  their  melts.  Studies  of  the  thermal 
conductivity  of  polymers  conducted  with  simultaneous  Investigation  of  vari¬ 
ations  or  structure  and  molecular  weight  allow  us  to  obtain  data  supple¬ 
menting  and  refining  modern  concepts  of  the  thermal  conductivity  method  and 
of  polymeric  bodies.  These  studies  are  necessary  both  for  advancing  the 
overall  theory  of  polymeric  properties  as  a  function  of  polymer  structure, 
and  also  for  substantiating  and  determining  the  main  parameters  of  processes 
of  obtaining  and  processing  polymers. 

The  aim  of  this  study  is  to  investigate  the  thermal  conductivity  of 
polycaproamide  as  a  function  of  temperature  and  molecular  weight  while 
simultaneously  studying  the  nature  of  supermolecular  structure. 

Determination  of  the  thermal  conductivity  was  conducted  by  the  method  of 
heating  two  plates  with  a  constant-power  flat  heat  source  [1],  The  thermal 

l 

conductivity  of  polycaproamide  was  studied  In  this  wor'k  as  a  function  of 
molecular  weight  in  the  temperature  range  20-250*C. 

The  test  specimens  were  made  in  I  the  form  of  flat  sheets  49.5  ram  in 
diameter  and  3-5  mm  thick.  By  special  processing  of  caprolactam  and  poly¬ 
caproamide  melts  —  degassifying  under  vacuum  —  monolithic  crystalline 
blocks  were  obtained  practically  free  of  moisture  (not  more  than  0.01  percent 
by  weight).  The  specimens  of  captolactsm-monomer  were  obtained  by  crystal¬ 
lizing  the  melt  in  a  metal  cell  with  no  access  to  air.  The  blocks  of 
polycaproamide  with  the  required  molecular  weight  were  obtained  by  polymer¬ 
ization  of  caprolactam  in  a  special  autoclave  under  specific  sets  of 
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conditions.  The  averaged  molecular  welg.  calculated  from  the  curves  of  the 
molecular-weight  distribution  was  adopted  as  the  molecular  characteristic  of 
polycaproamlda.  These  curves  were  obtained  by  the  method  of  diffusive 
salting  out  of  a  polycaproamlde  solution  developed  in  the  Institute  of 
Monomer  and  Polymer  Chemistry  of  the  Academy  of  Sciences  Ukrainian  SSR. 

The  melting  point  of  polycaproamlde  was  determined  with  the  aid  of  the  MIN-8 
microscope  by  observing  the  fusing  process  in  polarized  light.  To  exclude 
the  effect  of  crystallinity  on  thermal  conductivity  when  comparing  specimens 
of  polycaproamlde  of  different  structures  and  different  molecular  weights, 
specimens  with  identical  density,  determined  by  the  gradient  tube  method, 
were  selected. 

Physicochemical  characteristics  of  the  test  polycaproamlde  specimens 
are  listed  in  the  tabic.  Microphotographs  of  the  superaolecular  structure 
of  polycaproamlde  were  obtained  on  the  MRI-6  microscope  in  transmitted 
polarized  light. 

The  thermal  conductivity  of  caprolactam  and  polycaprolactam  as  a 
function  of  temperature  is  presented  in  Figures  1  and  2.  The  general  nature 
of  the  curves  describing  thermal  conductivity  as  a  function  of  temperature 
corresponds  to  the  concepts  of  heat  transfer  in  low-molecular  and  polymeric 
compounds.  It  follows  from  Figure  1  that  thermal  conductivity  of  capro¬ 
lactam  varies  slightly  with  temperature  rises  in  the  temperature  range  up  to 
the  melting  point  (about  69.5°C)  thermal  conductivity  rises,  but  at 
temperatures  higher  than  the  melting  point  it  drops  off  with  temperature 
rise.  For  polycap rolactasn,  thermal  conductivity  was  observed  as  a  function 
of  temperature  (Figure  2)  characteristic  for  most  polymers  [2,  3],  namely: 
thermal  conductivity  decreases  with  temperature  rise. 

In  this  present  work,  the  thermal  conductivity  of  caprolactam  and 
polycaprolactam  were  studied  in  temperature  ranges  within  which  the  first- 
order  structural  transition  takes  place  —  melting.  The  presence  of  this 
structural  transition  is  due  to  the  characteristic  maximum  on  the  thermal 
conductivity  versus  temperature  curve  for  caprolactam  (Figure  1)  and  the 
well  defined,  successive,  minimum  :  ad  maximum  on  the  very  same  curves  for 
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polycaproamide  (Figure  2) .  The  sharp  variation  of  many  properties  in  this 

case,  including  thermal  conductivity,  results  from  the  abrupt  variation  of 

the  energy  state  of  the  system.  The  effect  of  polymer  fusion  on  the  heat 

transfer  process  is  expressed,  first  of  all,  in  the  abrupt  drop  of  the 

3  5 

modulus  of  elasticity  of  the  macromolecules  (by  10  -10  times  [5]),  and  also 
in  the  considerable  rise  of  mobility  of  individual  segments  and  of  the  macro¬ 
molecules  as  a  whole. 

Table 

No.:  Mean  Molecu-:  Extent  :  Melting  :  Density  i  Size  of 

on:  Tar  Weight  j of  Poly-:  Point,  ;  at  20°C  ;  Spherolites, 

!  :  dispersej-'C  :  :  p 


1 

18000 

— H5SS — 
0,20 

210 

1,1590 

150+200 

2 

21000 

0,39 

218 

1,1610 

50*70 

3 

25000 

0,28 

222  * 

1,1593 

200+220 

4 

32000 

0,27 

215 

1,1590 

I 00+120 

Near  the  melting  point  in  polycaproamide  melt,  high  order  of  macro¬ 
molecules  is  possible,  as  a  result  of  which  the  thermal  conductivity  of  the 
polymer  must  rise  sharply  all  the  way  to  values  corresponding  to  the  crystal¬ 
line  state  (Figure  2) .  The  probability  of  the  highly  ordered  state  of 
polycaproamide  melt  drops  off  rapidly  with  temperature  rise,  and  thermal 
conductivity  is  already  governed  by  the  physical  regularity  corresponding  to 
the  state  of  the  single-phase  liquid  system,  that  is,  decreases  with  temper¬ 
ature  rise. 

Analysis  of  curves  a,  b  and  d  in  Figure  2  lead  to  the  conclusion 
that  in  spite  of  some  differences  in  structure  of  specimens  there  is  still  a 
marked  tendency  toward  rise  of  polycaproamide  thermal  conductivity  with 
increase  in  molecular  weight.  This  function  is  observed  even  more  distinctly 
for  molten  polycaproamide,  where  the  effect  of  the  supermole cular  structure 
(spherolite  type)  is  tcluded. 
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It  is  clear  from  the  microphotographs  we  have  taken  (Figure  3)  that  all 
the  polycaproaaide  specimens  studied  have  a  crystalline  structure  in  the 
form  of  spherolltes  of  fibrillary  structure  and  differ  in  spherolite  size. 

By  comparing  Figures  2  and  3  it  follows  that  the  spherolite  sizes  have  an 
effect  on  the  polymer  thermal  conductivity.  This  effect  evidently  is 
accounted  for  by  the  fact  that  the  thermal  conductivity  of  polymerB, 
including  polycaproamide ,  are  determined  by  the  number  and  mobility  of 
structure  elements  (spherolltes).  All  possible  changes  in  the  mobility  of 
structural  elements  can  be  accounted  for  by  change  in  their  size.  When 
spherolite  dimensions  are  reduced,  their  mobility  and  numbers  per  unit 
volume  rise.  This  can  account  for  the  riae  in  thermal  conductivity  for  the 
specimens  a  and  b  compared  with  the  thermal  conductivity  for  the  specimen 
(Figure  2).  Upon  analyzing  the  effect  on  thermal  conductivity  of  the 
supermolecular  structure  type  and  the  dimensions  of  structural  elements,  we 
must  also  take  into  account  the  homogeneity  of  structure. 


Figure  1.  Thermal  conductivity  of  caprolactam  as  a  function  of 
temperature:  •,  Values  obtained  in  the  present  study;  ♦  ,  Values 
found  in  [4] 


Figure  2.  Thermal  conductivity  of  polycaproamide  as  a  function 
of  temperature  for  various  molecular  weights:  a,  M  ■  18,000; 
b,  M  -  21,000;  c,  M  -  25,000;  d,  M  -  32,000 


Figure  3.  Microphotographs  of  the  structure  of  test  polycapro¬ 
amide  specimens:  a,  M  -  18,000;  b,  M  ■  21,000;  c,  M  ■  25,000; 
d»  M  -  32,000 
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